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RESUMEN 



Se presentan los resultados de los ensayos de laboratorio realizados en dos aisladores 
Elastomericos. Se ha determinado el factor de amortiguamiento, la rigidez lateral y el modulo 
de code para diferentes deformaciones de la goma y diferentes frecuencias de excitacion. 

La mejor opcion para que una estructura tenga un gran comportamiento sfsmico es que 
tenga aisladores de base o disipadores de energia. Por este motivo es que en el Centro de 
Investigaciones Cientificas, CEINCI-ESPE, se inicio en el 2006 las investigaciones en este 
campo y ahora con el concurso de la empresa privada se ha podido construir y ensayar los 
primeros aisladores elastomericos. 



ABSTRACT 



We present the results of laboratory tests in Elastomeric Isolator. We have determined 
the damping factor, the lateral stiffness and the shear modulus for different deformations of the 
rubber and different excitation frequencies. 



The best option to have a great seismic behavior in a structure is to have base isolators 
or energy sinks. For this reason, in the Scientific Research Center, CEINCI-ESPE, in 2006 
research were pioneered in this field and now with support of private enterprise, it has been 
able to build and test the first isolator made in Ecuador. 

1. INTRODUCCION 

En el mega sismo de Chile, del 27 de febrero de 2010 cuya magnitud fue M w =8.8 

varios edificios de hormigon armado, ubicados en Santiago de Chile, con muros de code 
sufrieron gran dano. En la figura 1 solo se presenta el caso de dos edificios. 
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Figura 1 Fallas en edificios con muros de corte durante el sismo de Chile de 201 0. 



En la parte superior izquierda, de la figura 1, se muestra un edificio muy esbelto, con 
muros de corte en forma de “ L” , “ T “.En la parte superior derecha se observa uno de los 
muros de corte de la planta baja, con una falla muy pronunciada a flexo compresion. Los 
edificios esbeltos no son adecuados en zonas sfsmicas y mas aun si es irregular en planta 
como en el presente caso. Aguiar (2010) 

En la parte inferior izquierda de la figura 1 se presenta otro edificio con muro de corte 
en el cual se presenta la falla por deslizamiento-cizalla. En el grafico de la derecha se aprecia 
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mejor esta falla debido a que el muro al moverse en los dos sentidos tritura al hormigon 
quedando las varillas expuestas y con el peso propio llegan a fallar el hierro. Probablemente 
este edificio fue disenado con un valor de reduccion de las fuerzas sismicas R muy alto que 
condujo a tener armadura horizontal y vertical baja. 




Figura 2 Hospital Militar con aisladores de base elastomericos. 



En contraste con el mal comportamiento de algunos edificios, hay otros con aisladores 
de base o disipadores de energfa que tuvieron un gran comportamiento durante el sismo de 
2007, en Santiago de Chile. Entre ellos se tiene el Hospital Militar cuya fachada se muestra en 
la parte superior de la figura 2, a la izquierda se observa que es una estructura irregular en 
elevacion ya que tiene un voladizo muy pronunciado en su fachada; a la derecha en cambio se 
aprecia una columna circular a la cual no le llegan las vigas en uno de los pisos por motivos 
arquitectonicos, esto no es aconsejable hacer en zonas sismicas. 

En el subsuelo del hospital existe un parqueadero, fotograffa inferior de la figura 2, con 
lo que se hace mas vulnerable la estructura ya que se tendrfa la falla denominada piso blando. 
Pero este edificio fue construido con aisladores de base elastomericos y tuvo un excelente 
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comportamiento durante el sismo a tal punto que sus ocupantes la noche del 27 de febrero de 
2010 pensaban que se trataba de un temblor al cual estan muy acostumbrados. 

Sobre la columna B9 de la figura 2 se ha colocado el aislador de base, en el 
parqueadero del Hospital Militar. Notese que hay dos ejes de vigas, una debajo de los 
aisladores y otra encima de los aisladores; todo esto para que trabajen como un todo. 

2. AISLADOR ELASTOMERICO CON PERNO PASADO 

En Aguiar et al. (2008) se presenta la forma como se realiza el analisis sfsmico de 
estructuras sobre aisladores de base elastomericos. Tambien se indica la existencia de 
aisladores sin nucleo de plomo y con nucleo de plomo y las ventajas que conlleva su empleo. 
Los aisladores elastomericos, son una especie de sanduche vulcanizado de goma y acero y en 
sus partes externas se tiene una placa de acero de mayor dimension con sus respectivos 
orificios por donde pasan los pernos de anclaje superior e inferior que conectan a la estructura. 

En la figura 3 se muestra el prototipo del aislador elastomerico ensayado, el mismo 
que esta compuesto por 6 capas de goma de 9 mm. y las dos exteriores son de 11 mm, de 
espesor e intercaladas con 7 placas de acero de 2 mm., de espesor. De tal manera que la 
altura solo de las gomas es de 76 mm. En la parte superior e inferior se tiene una placa de 
acero de 1 .5 cm., de espesor. La figura 3 corresponde a un code del aislador para mostrar que 
en la parte central existe un agujero para el perno pasante que tiene 19 mm., de diametro, el 
mismo que en la parte exterior esta soldado a la placa de acero. 




Figura 3 Geometrla de Aislador ensayado. 



Para la fabricacion del prototipo se utilizo un compuesto de caucho sintetico de primera 
calidad, el cual a su vez permite soportar altas cargas dinamicas, traccion, friccion, desgaste 
con Shore A 65 ± 5. La fabricacion qufmica para llegar a este compuesto se resumen en: 
Caucho sintetico en un 80 %; reciclado de llantas usadas (carga) reemplazando al caolfn en 
10%; Qufmicos y otros en un 10 %. Este es el primer aislador elastomerico fabricado en el 
Ecuador, que como se vera mas adelante no tuvo un comportamiento adecuado en el 
laboratorio razon por la cual se procedio a construir otro aislador en el cual se mejoro la calidad 
de la goma, se cambio el perno pasado por un perno macizo y con este aislador se obtuvo 
mejores resultados. 

A la izquierda de la figura 4 se aprecia el prototipo, la forma de las placas externas 
responde a la forma como se va a adaptar el prototipo a la maquina en la cual se realiza el 
ensayo. A la derecha se observa el perno pasante que se lo coloco como una especie de 
nucleo pero el objetivo fundamental era darle mayor fijeza al aislador. Los ensayos que a 
continuacion se indican se realizaron en el aislador sin y con el perno pasante. 
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Figura 4 Aislador elastomerico prototipo y perno utilizado como nucleo. 

3. AISLADOR ELASTOMERICO CON PERNO MACIZO 

En la figura 5 se indican los dos prototipos ensayados. Como se indico el primero es 
con perno pasado y el segundo con perno macizo. Precisamente en la figura superior se 
aprecia a la izquierda el aislador con el perno pasado el mismo que fue soldado a la placa de 
acero; el aislador funciono bastante bien hasta que se rompio la suelda; por ese motivo se 
volvio a construir otro aislador con las mismas propiedades y caracterfsticas del primero pero 
ahora se dejo un perno macizo pasado el mismo que se vulcanizo, es por eso que en la figura 
de la derecha se aprecia la goma que esta vulcanizada, ya no se utilizo suelda. 




Figura 5 Prototipos de aisladores de base ensayados. 
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Por cierto, las fotograffas fueron tomadas despues de haber realizado los ensayos por 
ese motivo es que estan un poco despintadas. En la parte inferior de la figura 5 se muestra una 
vista lateral de los dos aisladores, aparentemente parece que no tuvieron dafio pero si lo 
tuvieron, durante los ensayos la goma empezo a figurarse pero claro esta que una vez 
terminado el ensayo se cierran las fisuras y aparentemente da la impresion de que estan muy 
bien. 

4. FACTOR DE AMORTIGUAMIENTO Y RIGIDEZ 

El factor de amortiguamiento % y la rigidez horizontal del aislador k H se obtienen de la 
curva de histeresis, de la siguiente manera. A la izquierda de la figura 6 se muestra la maxima 
energfa elastica que absorbe el aislador, E s ; y a la derecha se muestra la energfa disipada por 

el aislador en un ciclo de histeresis, sea E D esta energfa. En Aguiar (2002) se presenta la 
deduccion de la ecuacion con la cual se obtiene 



e _ 1 ^_p 

4 7t E s 



(1 ) 



Donde % es el factor de amortiguamiento del aislador; E s es la energfa de 
deformacion elastica; E D es la energfa de deformacion disipada inelasticamente. 



La rigidez secante k H se obtiene como el cuociente entre la fuerza maxima en un ciclo 
de histeresis dividido para la deformacion maxima. Con la nomenclatuda indicada a la derecha 
de la figura 6 se tiene: 



k 



+ 

H 



F.L 



d 



+ 

max 



F~ 

i - max 

ii 



d~ 



2 



( 2 ) 



El significado de las variables que intervienen en la ecuacion ( 2 ) estan indicados en la 
figura 6, a la derecha. 




Figura 6 Diagrama de histeresis y energfas elastica e inelastica. 



Revista CIENCIA Vol. 13 N.- 2 de 2010 



111 



5. ENSAYOS REALIZADOS EN AISLADOR CON PERNO PASADO 

El aislador de base fue sometido a cargas cfclicas armonicas en la que se varfa la 
frecuencia de vibracion £2 ; para deformaciones de la goma previamente establecidas. Como 
se indico la altura solo de las gomas es 76 mm. Ahora bien, con respecto a este valor se 
determina la deformacion maxima lateral que se desea alcanzar en cada ensayo. Si se desea 
por ejemplo que la deformacion lateral sea del 10% se tiene que el desplazamiento lateral 
maximo es 7.6 mm., en cada direccion. 




Deformaci6n (mm) 



Figura 7 Curvas de histeresis para £2 - 0.1 Hz ■ Aislador sin pasador pasante. 



Para cada frecuencia de la excitacion £2 y para cada deformacion lateral 
preestablecida se realizaron cinco ciclos de carga de tal manera que cada una de las graficas 
de las figuras 7 y 8 hay 5 curvas. 

En la figura 7 se muestran las curvas de histeresis halladas sin el pasador pasante 
para una frecuencia £2 = 0.1 Hzy para deformaciones laterales del 10, 20, 30 y 40%. La 

frecuencia de excitacion de 0.1 /fzsignifica que en un minuto se realizan 6 ciclos. Mientras 
mas baja es la frecuencia menos ciclos se dan por minuto. Los aisladores de base pueden 
deformarse hasta un 200% de la altura de la goma pero en la practica es recomendable que el 
diseno se realice para un 50% de deformacion. Priestley et al. (2007). 

En la figura 8 se tienen las curvas de histeresis halladas para la misma frecuencia de 
excitacion y para las mismas deformaciones laterales, de la figura 7 pero para el aislador con 
pasador pasante. Las curvas son muy parecidas pero en realidad se requiere una mayor 
fuerza en el aislador con pasador pasante para lograr igual deformacion que en el aislador sin 
pasador pasante. 

Para apreciar la diferencia en las curvas de histeresis de las figuras 7 y 8, se presenta 
en la tabla 1, los valores de las fuerzas requeridas para obtener determinados desplazamientos 
de las curvas de histeresis; todo esto con la nomenclatura indicada en la figura 9 en que se 
han identificado cuatro segmentos de un ciclo de histeresis. Las curvas 1 y 3 son de carga y la 
2 y 4 de descarga. 
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En la tabla 1 se aprecia que para obtener un mismo desplazamiento lateral, en el 
aislador con perno se debe proporcionar una mayor fuerza, con respecto al aislador sin perno 
pasante pero no es muy significativa esta diferencia. Esto se debe a que el perno pasante no 
estaba ejerciendo una presion, por decirlo de alguna manera, con las laminas de acero y goma 
pero estaba soldado a las placas exteriores. 




Figura 9 Nomenclatura utilizada para identificar un ciclo de histeresis. 



6. RESULTADOS EN AISLADOR CON PERNO PASADO 

En la figura 10 se presenta, en la parte superior, el factor de amortiguamiento y en la 
parte inferior la rigidez horizontal del aislador sin perno pasador, a la izquierda y del aislador 
con perno pasador a la derecha. 
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Tabla 1 Valores de fuerzas y deformaciones promedio. 
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Con relacion al factor de amortiguamiento se puede indicar que a mayor deformacion 
lateral el factor de amortiguamiento disminuye y esto se justifica al observar las curvas de 
histeresis indicadas en las figuras 7 y 8, en que los lazos de histeresis disminuyen en su ancho 
conforme se incrementa la deformacion lateral. 

Para una deformacion lateral del 50% el factor de amortiguamiento medio en el aislador 
sin perno pasante es 3.7 % y en el aislador con perno pasante es 4.7 %. Si bien es cierto los 
valores hallados son bajos no es menos cierto que para una primera experiencia son 
aceptables. 

Con respecto a la rigidez horizontal se puede indicar que en el aislador sin perno 
pasante, la rigidez tiende a bajar conforme se incrementa la deformacion, que era el 
comportamiento esperado. En cambio el comportamiento es lo contrario en el aislador con 
perno pasante cuyos resultados se indican en la parte inferior derecha de la figura 10. 

A partir de la rigidez horizontal k H se determino el modulo de code G a con la 
siguiente ecuacion: 



Donde A es el area de la goma que trabaja al code, las dimensiones con las cuales se 
calculo el area son 13.5 por 14.5 cm.; H es la altura solo de las gomas, para el prototipo 
analizado este valor es 7.6 cm. Con los valores de k H indicados en la grafica inferior izquierda 

de la figura 8 y al aplicar la ecuacion ( 3 ) se hallo los valores de G a que se indican en la figura 

10. Para una deformacion del 50 % se tiene que: G a ~ 5 kg / cm 2 . Este valor se incremento 
notablemente con el segundo aislador, como se vera posteriormente. 
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Figura 10 Variacion de amortiguamiento y rigidez en aislador con y sin perno pasante. 




Figura 1 1 Rigidez al code G a de la goma, inferidos de los ensayos realizados. 

Los ensayos se iniciaron con una frecuencia de vibracion de 0.05 Hz., y se detuvieron 
en 0.15 Hz., debido a que la goma se desprendio del acero, tema que tambien se mejoro en el 
segundo aislador. 
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7. RESULTADOS EN AISLADOR CON PERNO MACIZO 

Con las experiencias adquiridas en el ensayo del primer aislador elastomerico 
construido en el Ecuador se construyo uno nuevo cambiando el perno pasado con un perno 
macizo pero como se indico esta vulcanizado, ya no soldado. Se mejoro la calidad de la goma 
y se volvio a ensayar, vale la pena mostrar como se obtiene las curvas de histeresis en el 
laboratorio. 

Las dos primeras fotograffas de la figura 12 muestran a uno de los aisladores durante 
el ensayo a cargas cfclicas y al final se aprecia la curva de histeresis que dibuja la maquina. 
Del estudio de la curva de histeresis se obtiene el factor de amortiguamiento la rigidez 

horizontal del aislador k h ; y a partir de la rigidez horizontal se encuentra el modulo de code de 

la goma G a . 




Figura 1 2 Ensayo de Aislador y curvas de histeresis. 



Se realizaron ensayos para frecuencias de vibracion que varfan desde Cl - 0.05 Hz 
hasta Cl = 0.45 Hz . A esta ultima frecuencia empezo a aparecer fisuras en la goma por lo que 
se suspendio el ensayo. Se destaca que en el primer prototipo (aislador con perno pasado) los 
ensayos llegaron a Cl - 0.15 Hz por lo que ahora se obtuvo una mejora considerable en el 
vulcanizado. Se ha mejorado pero todavia no es lo deseado, se continuara trabajando. 

En las figuras 13 y 14 se muestran las curvas de histeresis, promedio, para 
Cl = 0.05 Hz y para Cl = 0.40 Hz., respectivamente. Para bajas frecuencias los ciclos de 

histeresis son simetricos con respecto al origen, eso se observa en la figura 13 pero para altas 
frecuencias ya no son simetricas como ve en la figura 14. Esta falta de simetrfa en las altas 
frecuencias conduce a tener resultados desconcertantes en la rigidez del prototipo. 
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Figura 15 Variation del amortiguamiento con la deformation. Aislador con Perno Macizo. 




Figura 16 Variation de la rigidez horizontal con la deformation. Aislador con Perno Macizo. 



En la figura 15 se muestra la variation del amortiguamiento con la deformation. Es 
importante notar que para una deformation del 10%, el amortiguamiento varfa desde un 7% 
hasta un 10%, cantidades nada despreciables ya que van a aportar un buen amortiguamiento a 
la estructura. 



Para una deformation del 50% el factor de amortiguamiento varfa entre el 5% y el 
6.5%, que no esta tan mal sin embargo se va seguir investigando para elevar ese 
amortiguamiento. 

La tendencia del factor de amortiguamiento es muy clara, a medida que se aumenta la 
deformation lateral de la goma; el factor de amortiguamiento decrece. 
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Figura 17 Variation del modulo de corte con la deformation. Aislador con Perno Macizo. 

En la figura 16 la variation de la rigidez horizontal del aislador con la deformacion 
lateral de la goma. La tendencia general es que la rigidez va disminuyendo con la deformacion 
pero como se indico para altas frecuencias existen puntos que tienen un comportamiento no 
esperado. 

En la figura 17 se muestra la variation del modulo de corte G a con la deformacion 
lateral de la goma, destacando una vez mas que el modulo de corte de la goma se obtuvo a 
partir de la rigidez horizontal k h ; de la altura efectiva de goma H y del area que trabaja al 

corte A. Para una deformacion lateral del 50% el modulo de corte G a varia entre 8 y 10 
kg / cm 2 pero existen un par de puntos que estan muy por arriba de los valores indicados. Los 
valores obtenidos de G a son muy buenos pero para tener una mejor certeza se debe hallar a 
partir de ensayos de corte solo de la goma. 

8. COMENTARIOS Y CONCLUSIONES 

Se han presentado los resultados de los ensayos de laboratorio realizados en dos 
aisladores Elastomericos, construidos en el Ecuador, el uno con un perno pasante y el otro con 
un perno macizo. Se ha determinado el factor de amortiguamiento, la rigidez lateral y el modulo 
de corte para diferentes deformaciones de la goma y diferentes frecuencias de excitation. Del 
estudio realizado se desprende lo siguiente. 

• El primer aislador construido fue con perno pasado, en esta experiencia la cabeza del 
perno se soldo a la placa de acero de tal manera que se tuvieron dos tipos de 
resultados, los primeros que corresponden al perno soldado y despues cuando se 
desprendio la suelda. De igual manera no se tuvo una vulcanization buena, esto se 
detecto cuando se empezo a realizar los ensayos a altas frecuencias. De todas 
maneras fue una experiencia positiva toda vez que se gano experiencia en el estudio. 

• El segundo aislador construido fue con perno macizo, este trabajo de mejor forma que 
el primero, para una deformacion lateral del 50%, el factor de amortiguamiento 
promedio de los ensayos realizados esta en el orden del 6%, cantidad aceptable pero 
todavia no es lo que esperamos. Fallo el vulcanizado pero esta vez a mas altas 
frecuencias, se llego a 0.45 Hz. 
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Los resultados alcanzados son bastante promisorios y nos llenan de optimismo, ya que 
se vislumbra de que en un tiempo cercano podemos tener un aislador elastomerico de muy 
buenas caracterfsticas a un costo bajo. 
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RESUMEN 

Se obtiene la matriz de rigidez de un elemento finito rectangular mixto de ocho nodos, 
para modelar el comportamiento elastico de un suelo saturado en dos dimensiones. Se 
considera dieciseis grados de libertad para definir los desplazamientos y cuatro grados de 
libertad para el campo de las presiones pero luego se expande la matriz de rigidez del 
elemento para considerar tres grados de libertad por nodo. 

Posteriormente, se obtiene la matriz de rigidez de un elemento finito utilizando el 
programa KelemQ8C0Sat y se indican todos los programas son necesarios para utilizar este 
programa. 



ABSTRACT 

The stiffness matrix of a quadratic (eight nodes) mixed quadrilateral finite element was 
obtained in order to describe the constitutive elastic behavior of saturated soils subjected to the 
hypothesis of plain strain. Sixteen degree of freedom were consider for displacement unknown 
and four degree of freedom for the pressure field. 

Later on, a Matlab code for this approach is presented and a simple academic example 
was added. 

1. INTRODUCCION 

Es muy amplio el campo de aplicacion del elemento finito de suelo saturado, una de 
ellas es en el analisis sfsmico de un puente con pilotes, en donde una parte del pilote se 
encuentra en el agua y otra parte se introduce en un suelo completamente saturado hasta 
llegar a un estrato resistente. En la figura 1 se indica uno de los apoyos de los apoyos del 
puente sobre el rfo Carrizal, el mismo que esta conformado por 9 pilotes de 2 m., de diametro y 
tienen una longitud de 54 m., aproximadamente. Una parte del pilote se encuentra en contacto 
con la atmosfera, otra con el agua y otra con un suelo saturado. 

En Aguiar y Garcia (2010) se modelo el suelo que rodea a los pilotes con resortes de 
comportamiento no lineal, la misma que se obtuvo del diagrams carga deformacion del suelo. 
Evidentemente modelar el suelo con resortes es aproximacion; lo mejor es hacerlo con 
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elementos finitos y por ahora se presenta la formulacion elastica dejando para un proximo 
arttculo la formulacion inelastica (Mroginski, 2008). 




Figura 1 Uno de los apoyos del Puente sobre el rfo Carrizal. 



2. MARCO TEORICO 

El modelado de suelos parcialmente saturados ha sido abordado a partir de diferentes 
enfoques teoricos (Mroginski et al. 2010). Sin embargo, las teorfas mas difundidas en la 
comunidad cienttfica se agrupan en dos grandes categorfas, la primera de ellas esta basada en 
un proceso de homogenizacion y teorfa de mezclas (Lewis & Schefler 1998) conduciendo a 
sistemas de ecuaciones no simetricas, mientras que la segunda emplea tecnicas de 
aproximacion para obtener las relaciones volumetricas, como ser una combinacion de estados 
(Di Rado et al. 2009) que permite alternar entre un sistema de ecuaciones matriciales 
simetricas o no simetricas dependiendo de la influencia de la relacion entre el grado de 
saturacion y la succion. 

En este arttculo se trabaja con el modelo simplificado aplicado a un suelo saturado. 
Este modelo se basa en el cumplimiento de las ecuaciones de equilibrio mecanico y de un 
modelo del flujo, a partir de los cuales se obtiene el siguiente sistema de ecuaciones 
diferenciales Beneyto et al. (2005). 



£ 7 .. -a lP S tj 


(1) 


div(p w S w V vs )=a x £ u —a n p w 


(2) 



Pw 

Sujeto a las siguientes condiciones de contorno para el desplazamiento y el flujo 

kk ™ t \t 

p w —(-Vp'+p'g) n = q- enr» (3) 

N T cr = t en r» (4) 

Donde o\, es la tasa del tensor de tensiones efectivo, a\, es la tasa del tensor de 
tensiones totales, aprovechando la simetrta de el tensor de tensiones puede expresarse 
como cr = {cr x cr y r xy Y para problemas bidimensionales; a x y a u son constantes que 
dependen de las propiedades materiales del suelo; p w es el vector de velocidad de la presion 
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del agua; S l es la funcion delta de Dirac que vale la unidad para i = j , y vale cero para i * / ; 
/\, es la densidad del agua; S w es el grado de saturacion del agua que varfa de cero a uno. 
Para un suelo saturado S w = I ; v v s es la velocidad relativa del agua respecto del suelo; es 
un elemento del vector de velocidades de las deformaciones, para el caso piano 
£ = {£ x e y y xy \ ■ Por otro lado q w es el vector de flujo del agua; g la aceleracion de la 
gravedad; div(*) la divergencia de la funcion (•). Los parametros a x y a n para un suelo 
saturado valen (Beneyto et al. 2005): 



a x - a 

n a-n 

a u — 1 

K w K s 

Donde a es la constante de Biot que se obtiene con la ecuacion (7); n es la 
porosidad que se halla con la ecuacion (8); K w es el modulo de masa de la fase Ifquida y K s 
es el modulo de masa de la fase solida. 



( 5 ) 

( 6 ) 



a = 1 - — 



n = 



e o +1 



( 7 ) 

( 8 ) 



Donde K t es el modulo de masa que se halla con la ecuacion (9); e 0 es la relacion de 
vacios inicial del material. 

( 9 ) 



A t = m C t m 

Donde para estado piano de deformaciones 

m T ={ 1 1 0 1 } 

Mientras que para el estado piano de tensiones, 

m T ={ 1 1 0 } 



( 10 ) 



( 11 ) 



Por otra parte, C t es la matriz constitutiva que relaciona el vector de deformaciones de 
la fase solida con el vector que engloba las tensiones efectivas definido en (1) 

3. FORMULACION MATRICIAL 

Las ecuaciones diferenciales (1) y (2) son discretizadas en el dominio espacial 
empleando el metodo de Galerkin e interpolando las variables primarias uy p con (Zienkiewicz 

& Taylor, 1989) 

( 12 ) 

“ = iV -" (13) 

P = N r p ’ 



Donde N u , N son vectores que contienen a las funciones de forma para la fase 
solida y para la fase Ifquida; u es el vector que contiene a los desplazamientos nodales; 
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p w e s el vector que contiene a las presiones de la fase liquida. Mas adelante se vuelve a 
escribir (12) y (13); con estas ecuaciones se obtienen las siguientes ecuaciones de gobierno 
del problema de consolidacion de suelos saturados. 



A partir de las ecuaciones (1) y (2) se discretiza el dominio con lo que se encuentra la 
nueva formulacion para el caso de un suelo saturado, para las fases solida y liquida. Estas son: 



[\B T C T BdQ)u-[\B T am T S w N p dQ.)p w N U T pg dCl+\ N U T t dT (14) 



{\N T u am T BN U da)u+ JiV; 



n a—n 

1 

K. K„ 



I N p dQ p w + 



jfvivj 1 



kk' 



P 



VN p dQ 






+ 



S (VNj^P’gdCl+lN* ^-gjr = 0 



P 



p 



(15) 



Donde: VN P es la divergencia de las funciones de forma de la fase liquida; k es el 
tensor de permeabilidad, para un problema piano es de (2x2) y para un problema en tres 
dimensiones de (3x3); k rw es la permeabilidad relativa del agua; p w es la viscosidad 
cinematica. 



Expresando matricialmente las ecuaciones de gobierno (14) y (15) se obtiene la matriz 
de rigidez del elemento y el vector de cargas nodales. 



K Qsw 


u 


0 P 

L^WS A wwJ 


pW 



o 

o 


u 


L° H wwJ 


pW 



-f u 

f W 



(16) 



Luego la matriz de rigidez de un elemento finito de suelo saturado es: 



K = 



K 

~ Qws 



-Q 



SW 






WW 



(17) 



Donde K es la matriz de rigidez del elemento finito solido; Q sw es la matriz de 
acoplamiento solido-lfquido; Q ws es la matriz de acoplamiento Ifquido-solido. Para el metodo 
simplificado: Qs« =Qm.Pww es la matriz de rigidez de la fase liquida. 

Por otra parte, las matrices que estan inmersas en la ecuacion (17) se obtienen de la 
siguiente manera: 



K=\\b t C t BJ dsdt 


(18) 


Qsw - a n mT N p J ds dt 


(19) 


Qws = a n m T B J ds dt 


(20) 


II 


(21) 


H.,=\\(vN r ) Tkk . VN,Jdsd, 


(22) 
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Donde B es la matriz de compatibilidad que relaciona el vector de deformaciones £ 
con el vector de desplazamientos d ; C T es la matriz de elasticidad del material; J es el 
determinante de la matriz Jacobiana; a n esta definido en la ecuacion (6); m definido en la 
ecuacion (10) y (11), en el apartado 6 se amplfa la explicacion sobre este vector; N P es el 
vector que contiene a las funciones de forma de la fase Ifquida. 

4. DESCRIPCION DE ELEMENTO FINITO 

Se va a trabajar con un elemento finito de 8 nodos para definir el campo de 
desplazamientos y con un elemento finito de 4 nodos para las presiones, de esta manera se 
cumple con la condicion de convergencia de Babuska-Brezzi (Lewis & Schefler 1998), que 
recomienda que el orden de los grados de libertad de los desplazamientos sea mayor al orden 
de los grados de libertad de las presiones, para garantizar convergencia en la solucion. En la 
figura 2 se muestran los dos elementos finitos por separado por didactica pero mas adelante se 
los une. 





,V7 


V6 


VS 


U7 


Ud 




US 


j 


i 


.. 


L 




Vs 




V4 


ua 








Ui ^ 


UZ _ 


U3 






l ' 


V=3 






pi 






Figura 2 Elementos finitos para los desplazamientos y las presiones. 

Estrictamente el modelo numerico de la figura 2, tiene 20 grados de libertad, debido a 
que las presiones son cero en los nodos intermedios pero por facilidad de solucion se trabaja 
con 24 grados de libertad (gl), como se ilustra en la figura 3. 




Figura 3 Grados de libertad de elemento finito mixto. 
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5. FUNCIONES DE FORMA 

En la figura 4 se presentan las deformadas elementales, en coordenadas naturales, 
con las que se obtienen las funciones de forma en los nodos centrales, para la fase solida. Los 
de la parte superior corresponden a las funciones 8 y 4; los de la parte inferior a las funciones 2 
y 6. En cada una de ellas se presenta la ecuacion de la parabola y de la recta, las mismas que 
se encuentran con el interpolador de Lagrange, para tres puntos y dos puntos respectivamente. 
Al multiplicar estas dos ecuaciones se halla las funciones de forma, que se indican a 
continuacion. 



</> 4 ( S ,t)=^a+ S )(i-t 2 ) 



zL 




Figura 4 Funciones de forma para los nodos centrales. 

En Onate (1992) y Awruch & Di Rado (1998) esta bastante bien detallado el calculo de 
las funciones de forma para los nodos esquineros, sin embargo de ello, se describe el calculo 
de la funcion de forma (f) x {s,t) en tres etapas. En la primera existe la deformada unitaria en el 

nodo con lo cual se halla (f> x {s,t) = 0.25(l-,s)(l-/); en esta deformada elemental los nodos 

intermedios adyacentes se desplazan 1/2 pero por definicion de la deformada elemental con la 
que se halla la funcion de forma deben ser cero. En la segunda etapa, para cumplir con la 
condicion de la deformada elemental a la funcion de forma (j> 2 {s,t) se multiplica por (-1/2), 

con lo que al sumar a (j) x {s,t ) se obtiene cero en el nodo 2. En la tercera etapa se tiene que 
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hacer cero en el nodo 8 para esto se multiplica la funcion de forma 0 8 (s,?) por (-1/2) y se 
suma, de tal manera que: 



0x U r ) = 0\ ( 5 , t)~0. s U A (24) 

Al reemplazar las funciones de forma indicadas a la derecha de la ecuacion anterior se 
encuentra, luego de factorar la expresion final de 0 X (s,t) . AL proceder en forma similar con los 
otros nodos esquineros se encuentra. 



^i(M) = ^-(l-s)(l-f)(-s-f-l) 

4 (25) 

05 fa ^ “ (l + S ) (l + t ) (s + t — l) 

<z> 7 {s, t) = ^ (l - 5 ) (l + 1) (- 5 + 1 - 1) 

En la figura 5 se indican las deformadas elementales de las funciones de forma de los 
nodos esquineros. En la parte superior se tiene para los nodos 1 y 3; en la parte inferior para 
los nodos 5 y 7. 




Figura 5 Funciones de forma de los nodos esquineros. 



En la ecuacion (12) se agrupo las funciones de forma en el vector N u . Luego: 

N u = {0 (s, t) 02 (s,t) 00S, t) 00s, t) 00s, t) 00s, t) 00s, t) 0 8 Cs,O} ( 26 ) 

Las ordenadas de la elastica de un punto interior del elemento finito de ocho nodos, 
son las siguientes: 
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u(s,t)= (j. \(s,t)+u 2 fa(s,t)+U 3 fa(s,t)+U 4 fa(s,t) + 

+ u 5 fa {s, t ) + u 6 fa {s, t) + U-j fa {s, t) + u s fa (5, t) ^ 
v{s, t) = Vj fa (.S-, t) + v 2 fa {s, t) + v 3 fa {s, t) + v 4 fa (5, t) + 

+ v 5 0 b(s,f) + v 6 ^ 6 (j,r)+v 7 ^ 7 (s,f) + v 8 ^ 8 (s,f) 
Donde m, es la componente de desplazamiento horizontal del nodo i ; v,. es la 
componente de desplazamiento vertical del nodo i. En la ecuacion (12) se denomino u al 
vector que contiene a u(s,t), v(s,t ) . 

Las funciones de forma para el elemento finito de cuatro nodos con el cual se trabaja la 
fase Ifquida, son: 



^M = ^(l + s)(l-r) 

fa 3 (s,t) = — (l + 5')(l + r) 



(28) 



En la figura 6 se presentan las deformadas elementales para encontrar las funciones 
de forma de la fase Ifquida, para los nodos 1 y 2 en la parte superior y para los nodos 3 y 4 en 
la parte inferior. Notese que en todos los casos se tiene dos ecuaciones de dos rectas ya que 
se tienen dos condiciones; el producto de estas ecuaciones reporta las funciones de forma para 
la fase Ifquida. 

t 

Z.. 






Figura 6 Funciones de forma para la fase Ifquida. 



La presion p(s,t) en un punto interior del elemento finito, esta definida por: 



p{s,t)= p 1 fa' (s,t)+ p 2 fa v (s,t)+ p 3 fa v (s, t) + p 4 fa 4 {s,t) 



(29) 
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Donde p i es la presion en el nodo esquinero i . Por otra parte, el vector N p vale: 

Np = V"' (s,t) 0 2 (s,t) 0" (, S,t ) 01 Cm)} 

(30) 

Con relacion a la ecuacion (13) se debe indicar que el vector p w - {p l p 2 p 3 p 4 }. 

6. MATRIZ DE COMPATIBILIDAD B 

Por definicion de la matriz B se tiene: 

e = Bd (31) 

Donde £ es el vector de deformaciones y d es el vector que contiene los 
desplazamientos nodales. Para el caso piano el vector £ vale: 

du 
dx 
dv 

dy 

du dv 
dy dx 

El vector de desplazamientos, vale: 

d = {«, V! u 2 v 2 ••• u s vj (33) 

Luego la matriz B tiene 3 filas y 16 columnas. Se halla midiendo las deformaciones £ 
correspondientes a las deformadas elementales. Asf por ejemplo, la primera columna de la 
matriz B se hallara midiendo las deformaciones correspondientes a la deformada elemental 
Mj — 1 y las demas nulas. La segunda columna de B se halla de la deformada elemental 

Vj = 1 y las demas nulas. Sea B x la matriz que contiene a los valores de las deformaciones 
de las dos primeras columnas de B , que corresponden al nodo 1. 





d0 1 {s,t) 


0 


dx 


0 


d0 l (s,t) 


dy 


d0 1 {s,t) 


d0, (s,t) 


dy 


dx 



(34) 



Por lo tanto, la matriz B vale: 

B = [#1 B 2 B 3 B 4 B 5 B 6 B 1 B S ] 



( 35 ) 
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De tal manera la matriz B del nodo i vale: 



B,= 



dx 



0 

d0iO>J) 

dy 

dfiisj) d^(s,t) 



0 



dy dx 

7. MATRIZ JACOBIANA Y SU INVERSA 

Se define la matriz Jacobiana J de la siguiente manera: 



J = 



dx 

dx 

dt 



ds 

dy_ 

dt. 



( 36 ) 



(37) 



El determinante de la matriz Jacobiana es el area de influencia del punto de 
integracion, en el Metodo de la Cuadratura de Gauss. Para entender bien el significado de J , a 
continuacion se definen las coordenadas reales y las coordenadas naturales, tambien llamadas 
coordenadas computacionales. 

A la izquierda de la figura 7 se presentan las coordenadas reales de un elemento finito. 
Notese que se ha numerado al primer nodo en la parte inferior izquierda y se ha continuado la 
numeracion en forma antihoraria; las coordenadas del nodo 1, son (jc, , y, ) y ast 
sucesivamente. A la derecha de la figura 7 se muestran las coordenadas naturales, para el 
nodo 1, ahora le corresponde (-1,-1). Un punto cualquiera dentro del elemento finito se halla 
con las siguientes ecuaciones: 



X = ( 5 , t) + X 2 </> 2 ( S , t ) + X 3 </> 3 (5, t ) + x 4 0 4 ( S , t ) + 

+ x 3 (j) 5 (s, t ) + x 6 </> 6 (s, t ) + x 7 </> 7 (s, t) + X S 0 S (s, t ) 
y = yA (s, t ) + y 2 (/) 2 (s, t ) + y 3 (j) 3 (s, t ) + y 4 </> 4 (s, t ) + 

+ y 5 <f> 5 (s, t ) + y 6 (j) 6 (s, t ) + y 7 (j) 7 (5, t) + y 8 </> s (s, t ) 



(38) 




(- 1 , 1 ) 




M’-i) (0,-1) 

Figura 7 Coordenadas reales y naturales de un elemento finito Q8. 
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Para encontrar la matriz de compatibilidad B se debe tener en cuenta, lo siguiente: 



d(j) i _ d(f) i dx d(j) i dy 

ds dx ds dy ds 

Mi _ d(j), dx | d(j) ] dy 

dt dx dt dy dt 



( 39 ) 



Escrito en forma matricial, se tiene: 



'Mi ’ 




dx 


dy 




ds 






ds 


dx 


Mi 




dx 


dy 


Mi 


_ dt _ 




.dt 


dt_ 


_ dx _ 



~Mi~ 




ds 


— T 


Mi 


— J 


. dt _ 





Mi 

dx 

M,- 

dx 



(40) 



De donde la derivada de las funciones de forma respecto de las coordenadas reales 
puede obtenerse de: 



Mi 




Mi 




dx 


7-1 


ds 


(41) 


Mi 


= J 


Mi 


_ M _ 




. dt _ 





De tal manera que para hallar la matriz de compatibilidad B y la matriz de 
permeabilidad (22) se debe encontrar las funciones de forma con respecto a las coordenadas 
naturales, luego la matriz Jacobiana y finalmente se obtienen las derivadas de las funciones de 
forma en coordenadas cartesianas, mediante la ecuacion (41). 

Para terminar con el calculo de la matriz de rigidez de la fase solida, se debe indicar 
que la matriz de elasticidad para un problema de deformacion plana es (Oliver, 2000): 



C T 



e [ l-v) 

(l + v)(l-2v) 



1-v 

V 



0 



V 

1-v 

0 



0 

0 

l-2v 

2 



(42) 



Donde V es el modulo de poisson; E es el modulo de elasticidad del suelo. 



8. MATRICES DE ACOPLAMIENTO Y DE RIGIDEZ DEL FLUIDO 

La matriz de acoplamiento Q sw se puede hallar utilizando la ecuacion (19) o la (20). 
Los terminos que intervienen en esta ecuacion son: la matriz de compatibilidad B que tiene 3 
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filas y 16 columnas; a u que es un escalar; m es un vector de 4 columnas y 1 fila; N P es el 
vector que contiene las funciones de forma para la fase Ifquida y tiene 1 fila y 4 columnas. A 
continuacion se vuelve a copiar nuevamente la ecuacion (19) con la cual se halla Q sw . 



Qsw = {{^ r a n m T N P J ds dt ( 43 ) 

Para poder realizar el producto matricial indicado la matriz de compatibilidad B debe 
tener 4 filas en lugar de 3 filas, por lo que se aumenta una fila mas con ceros. Esta ultima fila 
corresponde a las deformaciones en sentido perpendicular al piano £ z que son nulas ya que el 
problema es en dos dimensiones. 

Con estas acotaciones, para el calculo de la matriz de acoplamiento Q sw se realiza el 
producto matricial indicado en la ecuacion (19) y se integra aplicando la cuadratura de Gauss. 

Para un suelo saturado se encuentra la matriz de rigidez del fluido aplicando la 

ecuacion (21), previamente se debe realizar el producto matricial indicado en dicha ecuacion e 
integrar aplicando la cuadratura de Gauss. Para tener presente su calculo se copia de nuevo la 
ecuacion (21). 

9. EXPANSION DE LA MATRIZ DE RIGIDEZ 

La matriz de rigidez del elemento finito, K e para suelo saturado, definido en la 

ecuacion (17) es de 20 por 20 , ya que la matriz de rigidez K es de 16 por 16; la matriz de 
acoplamiento Q sw de 16 por 4 y la matriz de rigidez del fluido P ww de 4 por 4. Ahora bien, 

esta matriz de rigidez se debe expander a los 24 grados de libertad del elemento indicado en la 
figura 3, para ello se procede en la forma indicada en la figura 8. 

En la parte superior izquierda de la figura 8, se muestra la matriz de rigidez de 20 por 
20 que inicialmente se obtiene, ahf se pueden ver las diferentes matrices que la conforman. 
Luego en la parte superior derecha se ve que se a anadido una matriz identidad cuadrada de 4 
por 4 y dos matrices cero de 20 por 4 y su transpuesta. La matriz identidad incrementada es 
para los grados de libertad de las presiones que se consideran en los nodos intermedios del 
elemento finito, que son cero, de tal manera que el vector de colocacion de las presiones para 
los nodos intermedios es cero, por lo que no influye el haber colocado la unidad. 

Ahora se debe obtener la matriz de rigidez del elemento finito de acuerdo a la 
numeracion de los grados de libertad indicados en la figura 3, para ello se procede a 
intercambiar filas y columnas; en la parte inferior izquierda de la figura 3 se muestra el 
intercambio de columnas y en la parte inferior derecha el intercambio de filas, se ha utilizado 
diferentes colores para ilustrar el cambio, sin embargo a continuacion se explica brevemente la 
forma como se cambian las columnas. 

Las dos primeras columnas corresponden a la matriz K luego quedan como esta, 
luego se puede ver en la figura 3 que el grado de libertad 3 corresponde a la presion por lo que 
se toma la primera columna de la matrices Q sw , P sw y se coloca a continuacion; las dos 

siguientes columnas corresponden a la matriz K por lo que queda igual; luego para el nodo 
intermedio la presion es cero por lo que se toma la primera columna de la matriz ceros y la 
identidad que se anadio y se intercambia; asi se continua con el intercambio de las columnas. 
Luego en el mismo orden se intercambian las filas. 
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La justification del intercambio de filas y columnas esta en la matriz de transformacion 
de coordenadas, descrito con detalle en Aguiar (2004). No solo se debe intercambiar las filas y 
columnas de la matriz de rigidez sino tambien las filas del vector de cargas. 
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Figura 8 Expansion de la matriz de rigidez del elemento finito. 



10.EJEMPLO DE APLICACION 

Se desea encontrar la matriz de rigidez del elemento finito indicado en la figura 9, 
considerando 9 nodos, el elemento es de 1 por 1 y se han numerado los nodos en la forma 
indicada, para que la matriz de rigidez tenga el menor ancho de banda. El modulo de 
elasticidad del material E = 1000 ; el modulo de Poisson v = 0.3 ; el espesor del elemento es 
e = 1.0; la relacion de vacios del material e 0 = 2; la rigidez del suelo k s =1000000; la 

rigidez del agua k w =1000000; la permeabilidad del suelo en sentido X es igual a la 
permeabilidad en sentido Y, k x - k y =0.0000864. Se omiten las unidades por facilidad de 
explicacion. 

En el Anexo 1 se indica el programa KelemQ8C0Sat debidamente documentado, 
razon por la cual se omite la explicacion de los datos de entrada, unicamente se indica como se 
opera en la modalidad consola. 

» ielem=1 ; 

» matpr=[1000 0.3 1.0 2 1000000 1000000 0.0000864 0.0000864 0 0]; 

»timer=[0.10 0.10 50 0.01 0.50]; % Aplicacion de carga y datos de integracion. 
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» conec=[1 3 5 7 8 6 4 2 1]; % Numeracion de nodos en forma antihoraria. 
» coord. =[0.0 0.0; 0.0 0.5; 0.5 0.0; 0.0 1.0; 1.0 0.0; 0.5 1.0; 1.0 0.5; 1.0 1.0]; 
l» inc=1 ; 
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Figura 9 Descripcion de ejemplo numerico. 

La matriz de rigidez del elemento finito que reporta el programa es: 

Kelem = 

1.0e+003 * 



Columns 1 through 6 



0.9615 

0.4541 

0.0002 

-1.1538 

-0.1496 

0.0000 

0.5128 

-0.0160 

-0.0000 

-0.3205 

-0.1068 

0.0000 

0.4808 

0.1870 

- 0.0001 

-0.6410 

-0.1068 

0.0000 

0.3526 

0.0160 

0.0001 

-0.1923 

-0.2778 

0.0000 



0.4541 

0.9615 

0.0002 

-0.2778 

-0.1923 

0.0000 

0.0160 

0.3526 

0.0001 

-0.1068 

-0.6410 

0.0000 

0.1870 

0.4808 

- 0.0001 

-0.1068 

-0.3205 

0.0000 

-0.0160 

0.5128 

-0.0000 

-0.1496 

-1.1538 

0.0000 



0.0002 

0.0002 

-0.0000 

- 0.0002 

0.0002 

0.0000 

0.0000 

0.0001 

-0.0000 

- 0.0002 

- 0.0001 

0.0000 

0.0001 

0.0001 

- 0.0000 

- 0.0001 

- 0.0002 

0.0000 

0.0001 

0.0000 

- 0.0000 

0.0002 

- 0.0002 

0.0000 



-1.1538 
-0.2778 
- 0.0002 
2.5641 
0.0000 
0.0000 
-1.1538 
0.2778 
0.0002 
0.0000 
-0.4274 
0.0000 
-0.6410 
-0.1068 
0.0001 
1 .0256 
-0.0000 
0.0000 
-0.6410 
0.1068 
- 0.0001 
0.0000 
0.4274 
0.0000 



-0.1496 

-0.1923 

0.0002 

0.0000 

1.2821 

0.0000 

0.1496 

-0.1923 

0.0002 

-0.4274 

-0.0000 

0.0000 

-0.1068 

-0.3205 

0.0002 

-0.0000 

-0.2564 

0.0000 

0.1068 

-0.3205 

0.0002 

0.4274 

-0.0000 

0.0000 



0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 



Columns 7 through 12 



0.5128 

0.0160 

0.0000 

-1.1538 



-0.0160 

0.3526 

0.0001 

0.2778 



-0.0000 

0.0001 

-0.0000 

0.0002 



-0.3205 

-0.1068 

- 0.0002 

0.0000 



-0.1068 

-0.6410 

- 0.0001 

-0.4274 



0.0000 

0.0000 

0.0000 

0.0000 
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0.1496 


- 0.1923 


0.0002 


- 0.4274 


- 0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.9615 


- 0.4541 - 0.0002 


- 0.1923 


0.2778 


0.0000 


- 0.4541 


0.9615 


0.0002 


0.1496 


- 1.1538 


0.0000 


- 0.0002 


0.0002 - 0.0000 


- 0.0002 


- 0.0002 


0.0000 


- 0.1923 


0.1496 - 0.0002 


1.2821 


0.0000 


0.0000 


0.2778 


- 1.1538 - 0.0002 


0.0000 


2.5641 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0010 


0.3526 


0.0160 - 0.0001 


- 0.1923 


- 0.2778 


0.0000 


- 0.0160 


0.5128 


0.0000 


- 0.1496 


- 1.1538 


0.0000 


- 0.0001 


- 0.0000 - 0.0000 


- 0.0002 


0.0002 


0.0000 


- 0.6410 


0.1068 


0.0001 


0.0000 


0.4274 


0.0000 


0.1068 


- 0.3205 - 0.0002 


0.4274 


- 0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.4808 


- 0.1870 - 0.0001 


- 0.3205 


0.1068 


0.0000 


- 0.1870 


0.4808 


0.0001 


0.1068 


- 0.6410 


0.0000 


0.0001 


- 0.0001 - 0.0000 


- 0.0002 


0.0001 


0.0000 


- 0.3205 


0.1068 


0.0002 


- 0.2564 


0.0000 


0.0000 


0.1068 


- 0.6410 - 0.0001 


0.0000 


1 .0256 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 



Columns 13 through 18 



0.4808 


0.1870 - 0.0001 


- 0.6410 - 0.1068 


0.0000 


0.1870 


0.4808 - 0.0001 


- 0.1068 - 0.3205 


0.0000 


0.0001 


0.0001 - 0.0000 


- 0.0001 - 0.0002 


0.0000 


- 0.6410 


- 0.1068 


0.0001 


1 .0256 


0.0000 


0.0000 


- 0.1068 


- 0.3205 


0.0002 


0.0000 - 0.2564 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.3526 


- 0.0160 - 0.0001 


- 0.6410 


0.1068 


0.0000 


0.0160 


0.5128 - 0.0000 


0.1068 - 0.3205 


0.0000 


- 0.0001 


0.0000 - 0.0000 


0.0001 - 0.0002 


0.0000 


- 0.1923 


- 0.1496 - 0.0002 


0.0000 


0.4274 


0.0000 


- 0.2778 


- 1.1538 


0.0002 


0.4274 - 0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.9615 


0.4541 - 0.0002 


- 1.1538 - 0.1496 


0.0000 


0.4541 


0.9615 - 0.0002 


- 0.2778 - 0.1923 


0.0000 


- 0.0002 


- 0.0002 - 0.0000 


0.0002 - 0.0002 


0.0000 


- 1.1538 


- 0.2778 


0.0002 


2.5641 - 0.0000 


0.0000 


- 0.1496 


- 0.1923 - 0.0002 


- 0.0000 


1.2821 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0010 


0.5128 


0.0160 - 0.0000 


- 1.1538 


0.1496 


0.0000 


- 0.0160 


0.3526 - 0.0001 


0.2778 - 0.1923 


0.0000 


0.0000 


- 0.0001 - 0.0000 


- 0.0002 - 0.0002 


0.0000 


- 0.3205 


- 0.1068 


0.0002 


- 0.0000 - 0.4274 


0.0000 


- 0.1068 


- 0.6410 


0.0001 


- 0.4274 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 



Columns 19 through 24 



0.3526 


0.0160 


0.0001 


- 0.1923 - 0.2778 


0.0000 


- 0.0160 


0.5128 


- 0.0000 


- 0.1496 - 1.1538 


0.0000 


0.0001 


0.0000 


- 0.0000 


0.0002 - 0.0002 


0.0000 


- 0.6410 


0.1068 


- 0.0001 


0.0000 


0.4274 


0.0000 


0.1068 


- 0.3205 


0.0002 


0.4274 - 0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.4808 


- 0.1870 


0.0001 


- 0.3205 


0.1068 


0.0000 


- 0.1870 


0.4808 


- 0.0001 


0.1068 - 0.6410 


0.0000 
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-0.0001 


0.0001 -0.0000 


0.0002 


-0.0001 


0.0000 


-0.3205 


0.1068 -0.0002 


-0.2564 


0.0000 


0.0000 


0.1068 


-0.6410 


0.0001 


0.0000 


1 .0256 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.5128 


-0.0160 


0.0000 


-0.3205 


-0.1068 


0.0000 


0.0160 


0.3526 -0.0001 


-0.1068 


-0.6410 


0.0000 


-0.0000 


-0.0001 -0.0000 


0.0002 


0.0001 


0.0000 


-1.1538 


0.2778 -0.0002 


0.0000 


-0.4274 


0.0000 


0.1496 


-0.1923 -0.0002 


-0.4274 


-0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.9615 


-0.4541 


0.0002 


-0.1923 


0.2778 


0.0000 


-0.4541 


0.9615 -0.0002 


0.1496 


-1.1538 


0.0000 


0.0002 


-0.0002 -0.0000 


0.0002 


0.0002 


0.0000 


-0.1923 


0.1496 


0.0002 


1.2821 


-0.0000 


0.0000 


0.2778 


-1.1538 


0.0002 


0.0000 


2.5641 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0010 



En su ejecucion el programs llama a los programas que se indican en el Anexo 2 y son 
los siguientes: 

• gauss_points Que define los puntos y pesos de la cuadratura de Gauss. 

• FluidCoef Que determina los coeficientes a p a n 

• ShapeQCO Que obtiene las funciones de forma para la fase solida y la fase liquids, 

encuentra ademas las derivadas respecto de las coordenadas naturales. 

• Jacobian_Matrix Determina la Matriz Jacobiana, su determinants y las derivadas de 

las funciones de forma con respecto a sus coordenadas reales. 

(Ecuacion 41). 

• MatrizB Encuentra la matriz de compatibilidad para la fase solida. 

• MatrizD Determina la matriz de elasticidad para el caso de deformacion plana. 

• Expand Expande la matriz de rigidez de 20 por 20 a 24 por 24. 
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ANEXO 1 

function [Kelem,Felem] = KelemQ8C0Sat(ielem,matpr, timer, conec, coord, inc); 

% 

% | Kaa Csw | |Au| |0 0 | |u| |fu| Kaa = B'*D*B ; fa = -B'*Sigma 
% | ||| = | | | | + | | Kal = B'*D*m*H ; fe = N'*Fext 

% j Cws Sww+Hww | |Ap| | 0 Hww | |p| |fp| Kla = H'*n*D*B ; fl = F(Sigma,Lambda J Grad2Lambda)*H 
% Kll = H'*(n'*D*m + h)*H - g*H'*P ; 

% nota: todo integrado en el dV 

% 

% Programado por: Javier Mroginski 

% Universidad Nacional del Nordeste (UNNE) 

% 

%— -Datos de entrada— 

% ielem es el numero de elemento que se analiza. 

% matpr es las propiedades de los materiales son 8 numeros mas 2 ceros al final 
% timer son los 5 datos de como se aplica el escalon de carga, convergencia, integration, AT y theta 
% conec es la numeration de los nudos en forma antihorario, empezando en vertice inferior izquierdo. Termina con el 
% numero de material. 

% coord, son las coordenadas de los nudos, dos columnas, la X y la Y. 

% inc es el incremento de carga, en este trabajo se toma inc=1 por que se requiere de un solo escalon para ver la 
% matriz de rigidez del elemento. 

% nelem Cantidad de elementos (Para ejemplo 1) 

% nnodo Numero de nudos por elemento (8) 

% nnod4 Numero de nodos de esquina (4) 

% npoin Numero de nudos totales (8) en la malla 
% nload Cantidad de cargas por nudo (para este ejemplo nload= 3) 

% ngaus Orden de interaction (2) 

% ndime Dimension del problema 2 dimensiones o 3 dimensiones (2) 

% ndofn Grados de libertad por nudo (3) 

% nfixe Numero de nudos restringuidos (7) 

% numat Cantidad de materiales (1) 

% nine Cantidad de incrementos (1) 

% nstre Codigo para Tension plana es 3 para deformation plana es 4 
% ntype Codigo que identifica tipo de problema. Solido 1 Suelo Saturado 2. 

% 

nelem=1 ;nnodo=8; nnod4=4; npoin=8; nload=3; ngaus=2; ndime=2; ndofn=3; 
nfixe=7; numat=1; ninc=1; nstre=4; ntype=2; 

%— La siguiente sentencia se debe habilitar cuando se ejecute program. mat— 

% y borrar los datos que acabo de dar. 

%global nelem nnodo nnod4 npoin nload ngaus ndime ndofn nfixe numat nine nstre ntype 

Kaa = zeros(nnodo*ndime , nnodo*ndime) ; Csw = zeros(nnodo*ndime , nnod4) ; 

Hww = zeros(nnod4 , nnod4) ; Kpp = zeros(nnod4 , nnod4) ; 

fu = zeros(nnodo*ndime ,1) ; fw = zeros(nnod4 , 1) ; 

Nu = zeros(ndime,nnodo*ndime); 

[gausp,weigp] = gauss_points(ngaus); 

corlo = coord(conec(ielem,1 :1 :nnodo),:) ; 
corl4 = corlo([1 3 5 7],:) ; 
kgaus = 0 ; imat=conec(ielem,9); 

D = MatrizD(matpr(imat,:), nstre); 

DTime = timer(inc,4) ; Theta = timer(inc,5) ; 

Kx = matpr(imat,7)./10 ; Ky = matpr(imat,8)./10 ; %se divide por 10 por las unidades 
Kperm = [Kx 0.0 ; 0.0 Ky ] ; 
m = [1 1 0 1]'; 

dens = matpr(imat,9) ; denw = matpr(imat,10) ; 

Ks = matpr(imat,5) ; Kw = matpr(imat,6) ; 
eO = matpr(imat,4) ; % Relacion de vacios inicial 
n = e0/(e0+1 ) ; % Porosidad 
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Satgr = 1.0; % Grado de saturacion, Satgr=1 -> 1 00% saturado 

g = [0 9.81]'; % Vector gravedad (en m/s2) 

[alpha, beta] = FluidCoef(D,m,Ks,Kw,n, Satgr); 

for igaus=1 :ngaus 
psi=gausp(igaus) ; 
for jgaus=1 :ngaus 
kgaus = kgaus + 1 ; 
eta=gausp(jgaus); 

[shap8,deri8] = ShapeQC0(psi,eta,nnodo) ; % Funcion de forma y derivada para u 
Nu(1 ,1 :2:nnodo*ndime)=shap8 ; 

Nu(2,2:2:nnodo*ndime)=shap8 ; 

[shap4,deri4] = ShapeQC0(psi,eta,nnod4) ; % Funcion de forma y derivada para p 
[detJB,DerCa,MatJB] = Jacobian_Matrix(corlo,shap8,deri8,ndime,nnodo) ; 

difV8 = detJB*weigp(igaus)*weigp(jgaus)*matpr(imat,3) ; 

MatrB = MatrizB(ndime,nnodo,nstre,corlo,shap8,DerCa) ; 

Kaa = Kaa + MatrB'*D*MatrB*difV8 ; % Matriz de rigidez de la fase solida 
[detJB,DerCa,MatJB] = Jacobian_Matrix(corl4,shap4,deri4,ndime,nnod4) ; 
difV4 = detJB*weigp(igaus)*weigp(jgaus)*matpr(imat,3) ; 

Csw = Csw + alpha.*MatrB'*m*shap4.*difV4 ; % Matriz de acoplamiento solido-agua 

Smat= beta.*shap4'*shap4.*difV4 ; 

Hww = DerCa*Kperm*DerCa'.*difV4 ; 

Kpp = Kpp + (Smat + Theta*DTime*Hww) ; 

fu = fu + Nu'*g.*(dens*(n-1 )+denw*n).*difV8 ; 
fw = fw + DerCa*Kperm*g.*denw.*difV4 ; 

end %jgaus 
end %igaus 

Ke = [+Kaa -Csw zeros(16,4); 

-Csw' -Kpp zeros(4,4) ; 
zeros(4,20) eye(4,4)] ; 

Fe = [-fu ; fw ; ones(4,1)] ; 

[Kelem,Felem] = Expand(Ke,Fe,nnodo,ndofn,ntype,ndime) ; 

ANEXO 2 



function [gausp,weigp] = gauss_points(ngaus) 



if ngaus==1 
gausp(1)=0.0; 
weigp(1)=2.0; 
elseif ngaus==2 

gausp(1)= -1/(sqrt(3)); %-0.5773502691 89626 

gausp(2)= 1/(sqrt(3)); %0.5773502691 89626 

weigp(1)=1 .0; 

weigp(2)=1 .0; 

end 

return 



function [alpha, beta] = FluidCoef(D,m,Ks,Kw,n, Satgr); 

% Ks = 1000000.0 ; % Rigidez del grano 
% Kw = 1000000.0 ; % Rigidez del grano 
% n = 2/(2+ 1 ) ; % Porosidad 

% Satgr = 1.0; % Grado de saturacion, Satgr=1 -> 1 00% saturado 

DRCom = m'*D*m ; 

DRCom = DRCom/9 ; 

Lvalu = 1.0- DRCom/Ks ; 

alpha = Satgr*Lvalu ; 

beta = (Lvalu-n)*Satgr/Ks + n*Satgr/Kw ; 

function [shape, deriv] = ShapeQC0(psi,eta,nnodo) 

S = psi ; T=eta ; 

S2 = S*2.0 ; T2 = T*2.0 ; 

SS = S*S ; TT=T*T ; 
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ST=S*T ; SST=S*S*T ; STT=S*T*T ; ST2=S*T*2.0 ; 
if nnodo==4 

shape(l) = 0.25*(1 - S)*(1 -T); 
shape(2) = 0.25*(1 + S)*(1 -T); 
shape(3) = 0.25*(1 + S)*(1 +T); 
shape(4) = 0.25*(1 - S)*(1 + T); 
deriv(1 ,1) = -0.25*(1-T); 
deriv(2,1) = 0.25*(1-T); 
deriv(3,1) = 0.25*(1+T); 
deriv(4,1) = -0.25*(1+T); 
deriv(1,2) = -0.25*(1-S); 
deriv(2,2) = -0.25*(1+S); 
deriv(3,2) = 0.25*(1+S); 
deriv(4,2) = 0.25*(1-S); 

elseif nnodo==8 

shape(1 )=(-1 .0+ST +SS+TT-SST-STT)/4.0 ; 
shape(2)=(1 .0-T-SS+SST)/2.0 ; 

shape(3)=(-1 .0-ST+SS+TT-SST+STT)/4.0 ; 
shape(4)=(1 .0+S-TT-STT)/2.0 ; 

shape(5)=(-1 .0+ST+SS+TT+SST+STT)/4.0 ; 
shape(6)=(1 .0+T-SS-SST)/2.0 ; 

shape(7)=(-1 .0-ST+SS+TT+SST-STT)/4.0 ; 
shape(8)=(1 .0-S-TT+STT)/2.0 ; 

deriv(1 ,1 )=(T+S2-ST2-TT)/4.0 ; 
deriv(2,1)=-S+ST ; 

deriv(3,1 )=(-T+S2-ST2+TT)/4.0 ; 
deriv(4,1 )=(1 .0 -TT)/2.0 ; 

deriv(5,1 )=(T+S2+ST2+TT)/4.0 ; 
deriv(6,1)=-S-ST ; 

deriv(7,1 )=(-T+S2+ST2-TT)/4.0 ; 
deriv(8,1 )=(-1 .0 +TT)/2.0 ; 

deriv(1 ,2)=(S+T2-SS-ST2)/4.0 ; 
deriv(2,2)=(-1 .0 +SS)/2.0 ; 

deriv(3,2)=(-S+T2-SS+ST2)/4.0 ; 
deriv(4,2)=-T-ST ; 

deriv(5,2)=(S+T2+SS+ST2)/4.0 ; 
deriv(6,2)=(1 .0 -SS)/2.0 ; 

deriv(7,2)=(-S+T2+SS-ST2)/4.0 ; 
deriv(8,2)=-T+ST ; 

end 

return 



function [detJB,DerCa,MatJB] = Jacobian_Matrix(corlo, shape, deriv,ndime,nnodo) 

detJB=0 ; 

MatJB=zeros(ndime,ndime); 

DerCa=zeros(nnodo,ndime); 

%%%% Matriz Jacobiana %%% 
for idime=1 :ndime 
for jdime=1 :ndime 
for inodo=1 :nnodo 

MatJB(idime,jdinne)=MatJB(idimeJdinne)+deriv(inodo,idinne)... 

*corlo(inodo,jdime); 

end 

end 

end 

%%% Jacobiano %%% 
detJB = det(MatJB) ; 
if detJB<=0.0 
stop ; 
end 

%detJB = MatJB(1 ,1)*MatJB(2,2)-MatJB(2,1)*MatJB(1 ,2) ; 
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%%% Matriz Jacobiana Inversa %%% 

MatJBi = pinv(MatJB) ; 

%%% Derivada de las funciones de forma en coord, globales %% 
for idime=1 :ndime 
for jdime=1 :ndime 
for inodo=1 :nnodo 

DerCa(inodoJdime) = DerCa(inodoJdime) + ... 

MatJBi(idimeJdime)*deriv(inodoJdime) ; 
end 
end 
end 
return 



function [MatrB]=MatrizB(ndime,nnodo,nstre,corlo, Shape, DerCa) 

MatrB=zeros(3,ndime*nnodo); 

% gpcod = corlo'*Shape' ; 
igaus=0; 

for inodo=1 :nnodo 
jgaus = igaus + 1 ; 
igaus = jgaus + 1 ; 

MatrB(1 Jgaus) = DerCa(inodo,1); 

MatrB(1 Jgaus) = 0; 

MatrB(2Jgaus) = 0; 

MatrB(2Jgaus) = DerCa(inodo,2); 

MatrB(3, jgaus) = DerCa(inodo,2); 

MatrB(3Jgaus) = DerCa(inodo,1); 

MatrB(4, jgaus) = 0.0 ; 

MatrB(4Jgaus) = 0.0 ; 

End 



function [D] = MatrizD(matpr,nstre); 

E=matpr(1); 

nu=matpr(2); 



if nstre==3 

D = E/(1-nu A 2)*[1 nu 0 nu ; nu 1 0 nu ; 0 0 (1-nu)/2 0 ; nu nu 0 1]; 
elseif nstre==4 

D = E*(1-nu)/(1+nu)/(1-2*nu)*[1 nu/(1-nu)0 0 ; nu/(1-nu) 1 0 0 ; ... 

0 0(1 -2*nu)/(2*(1-nu))0 ; 00 0 1]; 

D(4,4)=1; 

end 



function [Kelem,Felem] = Expand(Ke,Fe,nnodo,ndofn,ntype,ndime) 

% Ordena y expande las matrices de Rigidez elementales y ordena los vectores de carga 

% global nelem nnodo nnod4 npoin nload ngaus ndime ndofn nfixe numat nine nstre ntype 

Kelem = zeros(nnodo*ndofn , nnodo*ndofn); 

Kcol = zeros(nnodo*ndofn , nnodo*ndofn); 

if ntype==1&ndofn==6 % Solidos Cl 
%%% Intercambio de COLUMNAS 
indl =1 ; ind2=9 ; ind3=25 ; ind4=33 ; 
for inodo=1 :nnodo 
for idofn=1 :ndofn 

igl = (inodo-1)*ndofn + idofn ; 
if rem(inodo,2)==1 % nodo impar (de esquina) 
if idofn<=ndime 

Kcol(:Jgl) = Ke(:,ind1 ) ; 
indl =ind1 +1 ; 
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else 

Kcol(:,igl) = Ke(:,ind2) ; 
ind2=ind2+1 ; 
end 

elseif rem(inodo,2)==0 % nodo par (intermedio) 
if idofn<=ndime 

Kcol(:,igl) = Ke(:,ind3) ; 
ind3=ind3+1 ; 
else 

Kcol(:,igl) = Ke(:,ind4) ; 
ind4=ind4+1 ; 
end 
end 
end 
end 

%%% Intercambio de FI LAS 
indl =1 ; ind2=9 ; ind3=25 ; ind4=33 ; 
for inodo=1 :nnodo 
for idofn=1 :ndofn 

igl = (inodo-1)*ndofn + idofn ; 
if rem(inodo,2)==1 % nodo impar (de esquina) 
if idofn<=ndime 

Kelem(igl,:) = Kcol(ind1 , :) ; 
indl =ind1 +1 ; 
else 

Kelem(igl,:) = Kcol(ind2,:) ; 
ind2=ind2+1 ; 
end 

elseif rem(inodo,2)==0 % nodo par (intermedio) 
if idofn<=ndime 

Kelem(igl,:) = Kcol(ind3,:) ; 
ind3=ind3+1 ; 
else 

Kelem(igl,:) = Kcol(ind4,:) ; 
ind4=ind4+1 ; 
end 
end 
end 
end 

elseif ntype==2&ndofn==3 % Suelos Saturados 

%%% Intercambio de COLUMNAS 
indl =1 ; ind2=17 ; ind3=21 ; 
for inodo=1 :nnodo 
for idofn=1 :ndofn 

igl = (inodo-1)*ndofn + idofn ; 
if rem(inodo,2)==1 % nodo impar (de esquina) 
if idofn<=ndime 

Kcol(:,igl) = Ke(:,ind1) ; 
indl =ind1 +1 ; 
else 

Kcol(:,igl) = Ke(:,ind2) ; 
ind2=ind2+1 ; 
end 

elseif rem(inodo,2)==0 % nodo par (intermedio) 
if idofn<=ndime 

Kcol(:,igl) = Ke(:,ind1) ; 
indl =ind1 +1 ; 
else 

Kcol(:,igl) = Ke(:,ind3) ; 
ind3=ind3+1 ; 
end 
end 
end 
end 

%%% Intercambio de FI LAS 
indl =1 ; ind2=17 ; ind3=21 ; 
for inodo=1 :nnodo 
for idofn=1 :ndofn 

igl = (inodo-1)*ndofn + idofn ; 
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if rem(inodo,2)==1 % nodo impar (de esquina) 
if idofn<=ndime 

Kelem(igl,:) = Kcol(ind1,:) ; 

Felem(igl, :) = Fe(ind1 ,:) ; 
ind1=ind1+1 ; 
else 

Kelem(igl,:) = Kcol(ind2,:) ; 

Felem(igl, :) = Fe(ind2,:) ; 
ind2=ind2+1 ; 
end 

elseif rem(inodo,2)==0 % nodo par (intermedio) 
if idofn<=ndime 

Kelem(igl,:) = Kcol(ind1 , :) ; 

Felem(igl, :) = Fe(ind1 ,:) ; 
ind1=ind1+1 ; 
else 

Kelem(igl,:) = Kcol(ind3,:) ; 

Felem(igl, :) = Fe(ind3,:) ; 
ind3=ind3+1 ; 
end 
end 
end 
end 




Revista CIENCIA Vol. 13 N.- 2 de 2010 



143 



IMPLEMENTACION DE UN ENSAYO PCR MULTIPLEX PARA 
DETECCION DE GENES VANA, VANB Y VANC RELACIONADOS 
CON RESISTENCIA A GLUCOPEPTIDOS EN ENTEROCOCCUS 
FAECIUM Y ENTEROCOCCUS FAECALIS 



Maria Elena Terreros (1) , Marcelo Grijalva (1) , Patricia Jimenez (2) 

(1) Laboratorio de Biotecnologia Humana 
Escuela Politecnica del Ejercito, Sangolqui, Ecuador 

^Laboratorio de Biotecnologia Animal 
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RESUMEN 

Se ha implementado un ensayo genetico molecular basado en PCR multiplex para la 
deteccion temprana, confiable y precisa de 3 genes relacionados con resistencia a vancomicina 
(vanA, vanB y vanC), en dos especies de enterococo clmicamente importantes, Enterococcus 
faecium y Enterococcus faecalis. 

El sistema se aplico a muestras clmicas previamente aisladas mediante tecnicas 
tradicionales de cultivo microbiologico y deteccion de resistencia a antibioticos por 
antibiograma. La caracterizacion molecular determino que 2 de las 10 muestras clmicas 
recibidas eran Enterococcus faecium de genotipo vanA. Las 8 restantes fueron caracterizadas 
como Enterococcus faecalis sensibles a glucopeptidos. Existio concordancia entre los datos 
obtenidos por tecnicas convencionales de microbiologfa clfnica y el ensayo molecular 
implementado. La aplicacion de esta tecnica brindara a los pacientes rapidez y seguridad en la 
determinacion de infecciones por enterococos resistentes a vancomicina (ERV). 



ABSTRACT 

We implemented a molecular genetic test based on multiplex PCR for the early reliable 
and accurate detection of three genes associated with resistance to vancomycin (vanA, vanB 
and vanC) in two clinically important species of enterococci, Enterococcus faecium and 
Enterococcus faecalis. 

The system was applied to clinical samples which were previously isolated using 
traditional microbiological culture along with detection of antibiotic resistance by susceptibility 
testing. The molecular characterization assay found two out of 10 clinical samples to be a 
Enterococcus faecium vanA genotype. The remaining eight were characterized as 
Enterococcus faecalis susceptible to glycopeptides. There was consistency between data 
obtained by conventional clinical microbiology techniques and molecular testing. The application 
of this technique will help provide patients with a fast and reliable method to identify infectious 
events involving vancomycin resistant enterococci (VRE). 

1. INTRODUCCION 

La Vancomicina es un antibiotico glucopeptfdico utilizado en el tratamiento de 
infecciones por bacterias Gram positivas en caso de resistencia o alergia a otros antibioticos. 
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Su uso se ve limitado por los efectos secundarios que causa en el paciente, ademas de 
la reciente aparicion de brotes de cepas resistentes que ha encendido la alerta en muchos 
hospitales alrededor del mundo (Perez, 2002). 

En el Ecuador los dos primeros casos de ERV se detectaron en el Hospital de las 
Fuerzas Armadas de Quito, en agosto y septiembre del 2001. El tercer caso fue reportado en el 
Hospital Carlos Andrade Marin en el mismo ano. En todos ellos, la caracterizacion 
microbiologica tradicional determino la presencia de Enterococcus faecium resistente a 
vancomicina (Zurita et al, 2002). 

Datos del National Nosocomial Infection Surveillance (NNIS) de los Estados Unidos de 
America, revelan un aumento en los porcentajes de ERV desde 1989, con tasas que en el ano 
2000 se acercaron al 20% de todos los enterococos aislados (Mundy et al, 2000). 

Se han descrito seis genotipos de resistencia a vancomicina: vanA, vanB, vanC , vanD, 
vanE y vanG. Todos ellos comparten un mecanismo de resistencia muy similar basado 
fundamentalmente en la sustitucion del dipeptido terminal D-alanina-D- alanina, bien por D- 
alanina-D-lactato o por D-alanina-D-serina en la cadena de peptidoglucano de la pared celular 
bacteriana (Cetinkaya et al, 2000). 

Los genotipos vanA y vanB, adquiridos y transferibles, son los unicos que han estado 
implicados en brotes hospitalarios, por lo que su presencia implica la toma rapida de medidas 
de control. vanC se ha reportado en menor medida pero existe una prevalencia relativamente 
alta de casos a nivel mundial. Los genotipos de resistencia vanD, vanE y vanG se han 
encontrado solo en unos pocos casos y su importancia clfnica y epidemiologica es menos 
clara, ya que hasta el momento no han estado implicados en brotes (Perez, 2002). 

Clark et al (1993), Hayden et al (1993) y Quintiliani, Evers & Courvalin (1993) 
demostraron que los genotipos de resistencia a los glucopeptidos no siempre corresponden al 
fenotipo de resistencia determinado por las pruebas de sensibilidad a los antimicrobianos. Es 
por ello que el estudio de la epidemiologfa de la resistencia a los glucopeptidos es quizas el 
mas beneficiado mediante una combinacion de metodos fenotfpicos y genotfpicos de 
caracterizacion. 

El desarrollo de metodos moleculares rapidos, de alta sensibilidad y especificidad como 
la tecnologfa de Reaccion en Cadena de Polimerasa (PCR), que permite la amplificacion 
enzimatica in vitro de un segmento de ADN especffico (Cartwright, 1994), ha mejorado el 
tiempo y la precision del diagnostico de determinantes de resistencia a antibioticos, ofreciendo 
resultados en apenas unas horas luego de recibida la muestra (Perez, 2002). 

La PCR multiple o multiplex es una variante de la PCR convencional en la que dos o 
mas blancos geneticos son simultaneamente amplificados en la misma reaccion utilizando 
varios pares de cebadores o primers (Henegariu et al, 1997). 

El primer ensayo basado en PCR multiplex para la deteccion de resistencia a 
glucopeptidos (Dutka Malen et al, 1995) fue mejorado en investigaciones posteriores aunque, 
debido a la reciente caracterizacion de los genotipos vanD, vanE y vanG, dichos genes no han 
sido incorporados de manera rutinaria en los ensayos moleculares estandar (Depardieu et al, 
2004). 



El estudio de Depardieu F., Perichon B., & Courvalin P. (2004), utiliza una PCR 
multiplex que puede ser ajustada a los equipos, materiales y condiciones de trabajo de un 
laboratorio en particular, permitiendo la identificacion simultanea de genero y genotipo de 
resistencia. Esta puede ser una herramienta de gran utilidad diagnostica para el estudio de la 
epidemiologfa molecular de la resistencia a glucopeptidos en enterococo en nuestro pafs. 
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2. METODOS 

Este proyecto fue llevado a cabo en el Laboratorio de Biotecnologia Humana del 
Departamento de Ciencias de la Vida - Carrera de Ingenierfa en Biotecnologia de la Escuela 
Politecnica del Ejercito, Campus Sangolquf. Se ha contado ademas, con la colaboracion del 
Hospital Carlos Andrade Marin (HCAM) del Instituto Ecuatoriano de Seguridad Social (IESS) a 
traves del servicio de Laboratorio Clfnico - Laboratorio de Microbiologfa. El HCAM es parte de 
la Red Nacional de Control y Monitoreo de Resistencia a Antibioticos. 

El ensayo molecular incluyo las siguientes fases metodologicas: 

2.1 . Obtencion de muestras biologicas 

El ADN de las cepas tipo resistentes a glucopeptidos fue ordenado a partir de la 
Coleccion CCM (Czech-Slovak Collection of Microorganisms) y provisto por los Laboratories 
GeneProof de Brno, Republica Checa. Las cepas fueron: Enterococcus faecium ATCC 51559 
( vanA ), Enterococcus faecalis CCM 4647 ( vanB ) y Enterococcus gallinarum CCM 4054 ( vanC ). 

Durante el mes de Julio del 2009, el Laboratorio de Microbiologfa del HCAM 
proporciono cepas bacterianas, aisladas de muestras clfnicas, provenientes de pacientes con 
infeccion por enterococo. Estas cepas bacterianas fueron identificadas mediante tecnicas 
tradicionales de caracterizacion microbiologica y determinacion de susceptibilidad a 
glucopeptidos a traves de cultivo en agar nutriente y antibiograma por el metodo de difusion 
con disco de vancomicina (30 pg). Se consideraron tres criterios de resistencia: las muestras 
con diametro de inhibicion < 14 mm se consideraron resistentes; entre 15 y 16 mm como 
resistente intermedio y un diametro de inhibicion > 17 se determino como sensible (Chrystal, 
2002). Por este metodo se recolectaron 2 muestras resistentes, 3 cepas en el rango de 
resistencia intermedia y 5 cepas sensibles. 

2.2. Extraccion y cuantificacion de ADN bacteriano 

Se opto por el metodo de extraccion mediante Kit comercial tanto para la obtencion de 
ADN genomico como plasmfdico, pues asegura resultados mas confiables al extraer ADN de 
alta calidad, integridad, concentracion y pureza, ademas de permitirnos la optimizacion del 
tiempo empleado en el proceso. Para la obtencion de ADN genomico bacteriano se utilizo un 
protocolo modificado de extraccion para QIAamp DNA Blood Mini Kit (Qiagen, Alemania) y para 
la extraccion de ADN plasmfdico se utilizo el kit comercial QuickLyse Miniprep (Qiagen, 
Alemania) bajo el protocolo adjunto con el kit. 

El ADN extrafdo, tanto el genomico como el plasmfdico fue sometido a electroforesis 
por 45 min a 120 V en geles de agarosa al 0.8% y visualizado mediante luz UV. 

La cuantificacion de ADN se llevo a cabo mediante fluorometrfa, en el equipo Qubit™ 
(Invitrogen, Brasil) con el kit Quant-iT™ dsDNA HS Assay Kit. 

2.3. PCR multiplex 

En base al estudio de Depardieu et at (2004) se seleccionaron y ordenaron 5 pares de 
primers especfficos para la amplificacion de 3 genes de resistencia a glucopeptidos vanA, vanB 
y vanC y para la identificacion de 2 especies clfnicamente importantes de enterococo, 
Enterococcus faecium y E. faecalis. 

Utilizando el ADN de las cepas tipo a una concentracion estandar de 20 ng/ul se 
realizaron las pruebas y optimizacion del sistema PCR multiplex. La puesta a punto del sistema 
se realizo mediante tecnicas de optimizacion multidimensionales, evaluando variables como la 
temperatura de annealing, la concentracion de cloruro de magnesio, variacion del numero de 
ciclos y la utilizacion de adyuvantes. 
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El sistema se optimizo bajo las siguientes condiciones: 1.5 pi de ADN en 50 pi de 
reaccion que comprende 20 pmol de cada uno de los cinco pares de primers, buffer de PCR lx 
(200 mM Tris-HCI pH 8.4, 500 mM KCI), 1.5 mM de MgCI 2 , 0.1% de Triton X-100, Ipl de 
albumina de suero bovino (BSA lOmg/ml), 100 pM de cada dNTP y 1U de Taq polimerasa 
(Invitrogen, Brasil). 

La amplificacion se llevo a cabo en un termociclador Swift maxi modelo MX-BLC-4 
(ESCO, Singapur) bajo el siguiente programa: denaturacion inicial por 3 minutos a 94 2 C y 30 
ciclos de amplificacion a 1 min a 94 2 C, 1 min a 57 2 C y 1 min 72 2 C, extension final a 72 2 C por 7 
min y enfriamiento a 4 2 C por 10 min. 

La visualizacion de los productos de PCR se realizo mediante corrida electroforetica en 
geles de agarosa al 1%, a 120 V por 90 minutos. 

2.4. Sensibilidad anahtica 

A fin de determinar el comportamiento del sistema en respuesta a variaciones de la 
concentracion de ADN y establecer el Ifmite de deteccion de nuestra PCR multiplex, se llevo a 
cabo un ensayo de sensibilidad analftica mediante la evaluacion del sistema en 10 diluciones 
seriadas del ADN, partiendo de una concentracion de 20 ng/ml y un factor de dilucion 1 :2. 

2.5. Analisis y reporte de resultados 

Se analizo la presencia o ausencia de productos de PCR multiplex de pesos 
determinados, permitiendo la identificacion del patogeno presente en la muestra y su 
caracterizacion genetica de susceptibilidad ante glucopeptidos. 



3. RESULTADOS Y DISCUSION 

3.1 . Extraccion y Cuantificacion de ADN 

La calidad y cantidad de ADN es una de las condiciones indispensables para garantizar 
la correcta amplificacion de secuencias de genes mediante la tecnica de PCR (Surzycki, 2000). 
Varios metodos de extraccion de ADN de bacterias Gram positivas se han descrito, desde 
tecnicas muy sencillas y rapidas como la suspension directa en el master mix de PCR (Patel et 
al, 1997) o el metodo de ebullicion (Bell et al, 1998) ambas con resultados poco confiables. 
Otros metodos como el de lisis alcalina ofrecen una extraccion mas eficaz para PCR multiplex 
(Dutka-Malen et al, 1995). En nuestro trabajo, la extraccion de ADN de las muestras clfnicas 
mediante kit comercial nos permitio obtener material genetico de adecuada integridad y 
concentracion tanto para ADN genomico como plasmfdico. Los rangos de concentracion 
obtenidos oscilaron entre 18.67 ng/ul y 76.33 ng/ul, para el ADN genomico; y 3.6 ng/ul a 17.3 
ng/ul, para el plasmfdico. 

3.2. PCR multiplex 

En la Figura 1 se pueden observar los productos de PCR para cada uno de los 
genotipos y especies bacterianas estudiados. Para Enterococcus faecium vanA (Carril 1) se 
observan dos productos, uno de 732 pb que corresponde a vanA y otro de 1091 pb que es 
especffico para E. faecium. En el Carril 2, de Enterococcus faecalis vanB, de igual forma 
podemos visualizar dos bandas, una de 647 pb, correspondiente al genotipo vanB y un 
producto de 475 pb, que nos indica la especie bacteriana E. faecalis. En el carril 3 tenemos a 
E. galiinarum vanC, en donde el unico producto de PCR de 815/827 pb senala su genotipo de 
resistencia vanC. 
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E. jhecim 1091 pb 



S15/S27pb ^ 

vanA 732 pb 
vanE 647 pb - 



E. faecalis 475 pb 




Figura 1 Electroforesis en gel de agarosa al 1% P/V tenido con bromuro de etidio de los 
productos de PCR multiplex optimizado para deteccion de genes de resistencia vanA, vanB y vanC e 
identificacion de Enterococcus faecalisy E. faecium. M = marcadorde 100 pb (Invitrogen, Brasil); Carriles: 
1 = Enterococcus faecium vanA , 2 = Enterococcus faecalis vanB, 3 = E. gallinarum vanC, C- = control 
negativo de la reaction. 

3.3. Sensibilidad analitica 

El ensayo de sensibilidad analitica determino que la concentracion minima de ADN que 
permite la deteccion de la especie faecium es de 1.25 ng/pl, por otra parte la deteccion del 
genotipo de resistencia vanA, es posible hasta una concentracion de 0.156 ng/pl (Figura 2). La 
deteccion de la especie faecalis y del genotipo de resistencia vanB es posible hasta una 
concentracion de ADN de 0.078 ng/pl (Figura 3). 

Ml 23 456 789 10 C- 

E. faecium 
1091 pb 
vanA - 
732 pb 




Figura 2 Ensayo de sensibilidad analitica para el genotipo de resistencia vanA y para 
identificacion de Enterococcus faecium. M = marcador de 100 pb (Invitrogen, Brasil); Carriles: 1 = ADN a 
20 ng/pl, 2 = 10 ng/pl, 3 = 5 ng/pl, 4 = 2.5 ng/pl, 5 = 1.25 ng/pl, 6 = 0.625 ng/pl, 7 = 0.312 ng/pl, 8 = 0.156 
ng/pl, 9 = 0.078 ng/pl, 10 = 0.04 ng/pl, C- = control negativo del ensayo. 






148 



Revista CIENCIA Vol. 13 N.- 2 de 2010 



M 1 2 3 45 67 89 10 C- 



vanB — 
647 pb 

E. faecalis 
475 pb 




Figura 3 Ensayo de sensibilidad analftica para el gen de resistencia vanB e identification de 
Enterococcus faecalis. M = marcador de 100 pb (Invitrogen, Brasil); Carriles: 1 = ADN 20 ng/pl, 2 = 10 
ng/pl, 3 = 5 ng/pl, 4 = 2.5 ng/pl, 5 = 1.25 ng/pl, 6 = 0.625 ng/pl, 7 = 0.312 ng/pl, 8 = 0.156 ng/pl, 9 = 0.078 
ng/pl, 10 = 0.04 ng/pl, C- = control negativo del ensayo. 



3.4. Analisis de resultados 

Se analizo el ADN extrafdo de las muestras clfnicas en ensayos por separado para 
ADN genomico y ADN plasrmdico. 

En la Figura 4 se muestra la PCR multiplex de ADN genomico, en la cual se observa la 
presencia de productos de PCR especfficos para Enterococcus faecium vanA en las muestras 
1 y 2. En las 8 restantes (3-10) se observa la presencia de un unico producto de PCR que las 
identifica como Enterococcus faecalis y no se presentan otros productos de PCR, lo que es 
indicativo de susceptibilidad a antibioticos. En la Figura 5 podemos observar la PCR multiplex 
de ADN plasrmdico extrafdo de las muestras clfnicas, donde se obtuvieron resultados similares. 

Ml 23456789 10 C- 



E. faecium 
1091 pb 



vanA 
732 pb 



E. faecalis 
475 pb 




Figura 4 PCR multiplex de ADN genomico de las muestras provenientes del HCAM. M = 
marcador de 100 pb (Invitrogen, Brasil); Carriles 1-10 = muestras clfnicas, C- = control negativo. 
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M 1 2345 6789 10 C- 



E. faecium 
1091 pb 



vanA . 
732 pb 



E. faecalis 
475 pb 




Figura 5 PCR multiplex con ADN plasmfdico de muestras provenientes del HCAM. M = 
marcador de 100 pb (Invitrogen, Brasil); Carriles: 1-10 = muestras clfnicas; C- = control negativo. 



Las dos muestras clfnicas que resultaron positivas a la caracterizacion de resistencia a 
vancomicina provinieron de muestras de orina. Las infecciones urinarias son la forma mas 
comun de enfermedad clfnica producida por enterococo y entre los pocos factores de virulencia 
descritos en esta especie se encuentra su capacidad de adherirse a las celulas epiteliales 
renales y aunque la naturaleza de las estructuras moleculares responsables de esta capacidad 
no se conoce bien, existen evidencias de que la adherencia se llevarfa a cabo a traves de una 
sustancia de agregacion producida en respuesta a feromonas y que podrfa contribuir a la 
habilidad de las cepas productoras para causar infecciones urinarias (Murray, 2000). 

El origen mas evidente de una infeccion urinaria es la transmision ano-rectal, por ello 
se aconseja la realizacion de un estudio de colonizacion gastrointestinal por ERV en nuestro 
pafs, pues como afirma Smith et al (2005) esta puede permanecer asintomatica por anos. 
Ademas, es importante tener en consideracion que la transferencia genetica bacteriana ocurre 
muy a menudo en el tracto digestivo y la contaminacion de la piel en pacientes reingresados a 
hospitales se produce facilmente (Mundy et al, 2000). 

Las cepas resistentes obtenidas en esta investigacion fueron identificados como 
Enterococcus faecium de genotipo vanA. No se encontraron genotipos de resistencia vanB ni 
vanC, pese a ello no se puede descartar que en un estudio de prevalencia que incluya un 
mayor campo muestral, en un periodo no menor a un ano, se detecten muestras con estos 
genotipos. 

Asimismo, es importante determinar si la aparicion de ERV en un hospital ha sido 
causada por la propagacion de una sola cepa o por la aparicion simultanea de diferentes 
clones, porque esto puede afectar la prioridad de las medidas de control a tomar. Una variedad 
de metodos de caracterizacion molecular se ha utilizado para establecer el grado de relacion 
clonal de ERV, incluyendo la ribotipificacion, el analisis de plasmidos, PFGE, arbitrarily-primed 
PCR, entre otras (Cetinkaya et al, 2000). Estas tecnicas podrfan ser incluidas en posteriores 
investigaciones relacionadas a este tema. 

La evidencia de la circulacion en el medio de especies de enterococo resistentes a la 
vancomicina constituye un problema a enfrentar en forma multidisciplinaria por todo el equipo 
de salud, tanto medicos clfnicos, microbiologos, infectologos, administradores de salud y 
personal de apoyo (Bazet et al, 2005). 
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En aquellos hospitales en los cuales los ERV aim no se han detectado o continuan 
siendo raros, como es nuestro caso, las medidas de control incluyen el realizar pruebas de 
sensibilidad a vancomicina en todos los aislados de enterococo procedentes de muestras 
clinicas, asi como la vigilancia periodica en heces o en hisopados rectales de pacientes de 
riesgo como unidades de cuidado intensivo y hemodialisis (Perez, 2002). 

En presencia de una cepa resistente a vancomicina es imprescindible su confirmacion 
molecular, por lo cual seria deseable contar con un centro de referencia nacional donde se 
realice este tipo de estudios moleculares (Bazet et al, 2005). 

La cooperacion y la comunicacion entre los laboratories y el programa de control de 
infecciones facilitan los esfuerzos de manera sustancial (Cetinkaya et al, 2000). 

Esta investigacion constituye el punto de partida para nuevas y mayores 
investigaciones que incluyan no solo el analisis de colonizacion ano-rectal de ERV sino la 
puesta en marcha de programas de vigilancia que incluyan la deteccion de S. aureus resistente 
a vancomicina mediante tecnicas moleculares. 



4. CONCLUSIONES 

• Hemos optimizado un metodo de extraccion y purificacion de DNA genomico bacteriano 
a partir de cultivos puros de enterococo, utilizando un protocolo modificado para el 
QIAmp DNA Blood Mini Kit (Qiagen), el cual nos ofrece resultados confiables, ADN de 
alta calidad, adecuada integridad, concentracion y pureza, en menor tiempo. 

• El protocolo estandar de extraccion del QuickLyse Miniprep (Qiagen) permite la 
adecuada extraccion de ADN plasmidico bacteriano para E. faecalis y E. faecium. 

• La concentracion minima de ADN que puede ser detectada con fiabilidad por nuestra 
PCR multiplex es: 1.25 ng/pl para la especie E. faeciurrr, 0.156 ng/pl para el genotipo 
vanA, y 0.078 ng/pl para la especie E. faecalis y el genotipo vanB. 

• En este trabajo fue posible estandarizar y optimizar un sistema molecular PCR 
multiplex a partir de cepas puras y aislados clinicos, que permite la deteccion de los 
genes vanA, vanB y vanC de resistencia a glucopeptidos en Enterococcus faecium y 
Enterococcus faecalis, evitando la generacion de productos inespecificos y de falsos 
positivos. Los productos de de PCR obtenidos corresponden a los siguientes pesos 
moleculares: 475 pb para E. faecalis, 647 pb para vanB, 732 pb para vanA, 815/827 pb 
para vanC y 1 091 pb para E. faecium. 

• De un total de 10 muestras clinicas, se obtuvieron 2 cepas resistentes a vancomicina. 
La caracterizacion molecular determino que se trataba de Enterococcus faecium vanA. 
Las 8 muestras restantes fueron caracterizadas como Enterococcus faecalis sensibles 
a glucopeptidos resultados que confirman los datos reportados por tecnicas 
tradicionales de cultivo microbiologico y antibiograma del HCAM. 

• Este trabajo servira de base para proyectos futuros que abarquen temas relacionados a 
ERV en nuestro pais. 
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RESUMEN 

El agua potable tiene dos orfgenes desde los cuales se suministran a las poblaciones: 
las aguas superficiales (nos, lagos y embalses) y las aguas subterraneas, a traves de pozos y 
acufferos. En el Ecuador, en la Region Sierra, las aguas provienen de acufferos o rfos. Algunos 
acufferos y aguas superficiales estan contaminados con arsenico como es el caso de la laguna 
de Papallacta y el no Tambo que proveen de agua potable a sectores de Quito. Tienen 
concentraciones por sobre los 10 pg/L que es la norma permitida en Ecuador y por La 
Organizacion Mundial de la Salud. El presente trabajo propone una alternativa para remover el 
arsenico de las aguas contaminadas mediante la utilizacion de un consorcio y de una cepa 
axenica de cianobacterias a nivel de laboratorio. Se utilizaron tres pretratamientos a la biomasa: 
calor a 92C, autoclave a 121 C y sin pretratamient o. Se encontro que las cianobacterias 
pretratadas con el autoclave presentaron un 61.1% de remocion al partir de una concentracion 
de arsenico de 100 pg/L. El Qmax mas elevado lo tuvieron las cianobacterias tratadas con 
autoclave con 171.59 pgAs/gCianobacteria adsorbidos. La prueba de toxicidad aguda en 
Daphnia pulex no mato a ninguna Daphnia por lo que se concluyo que el agua de las 
cianoabcterias no es toxica para ellas. 



ABSTRACT 

Drinking water comes from two sources, which are delivered to human populations: one, 
is surface water such as rivers, lakes and reservoirs. The other is groundwater, through wells and 
aquifers. In Ecuador, in the Sierra Region, the water comes from rivers or aquifers. Some ground 
water and surface water are contaminated with arsenic such as Papallacta lagoon and Tambo 
River wich provide drinking water to certain urban areas of Quito and they has concentrations 
above 10 pg/L (the standard allowed in Ecuador and the World Health Organization). 

This project proposes an alternative to remove arsenic from contaminated water by using 
a consortium and an axenic strain of cyanobacteria in the laboratory by following these steps: i) 
cultivation of biomass in 5 L reactors, ii) crop and pretreatment of biomass using three methods: 
autoclave at 121 °C and 15 psi of pressure, heat to 92 °C and without pretreatment, iii) contact 
of the biomass pretreated in water contaminated with arsenic at concentrations of: lOOpg/L, 300 
pg/L, 500 pg/L, 700 pg/L for 30 minutes. The water samples were collected by filtration and water 
samples analysis was made through hydride generation coupled with an atomic absorption 
equipment. The cyanobacteria pretreated with autoclave showed a 61 .1% removal of the arsenic 
concentration of 100 mg / L, the Qmax was higher when the cyanobacteria were treated with 
autoclave, this pretreatment showed a 171.59 pgAs/gcyanobacteria adsorbed. Cyanobacteria 
were tested for acute toxicity to Daphnia pulex in this thesis project, the acute toxicity test did not 
kill Daphnia pulex. 
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1. INTRODUCTION 

El arsenico (As) es un elemento quimico que esta en la naturaleza en distintas formas; 
muy pocas veces se lo encuentra en forma solida, y se lo puede hallar en plantas y animales 
en su forma organica que es la menos toxica. La mayorfa de las veces, el arsenico se combina 
con el oxfgeno, cloro y azufre para formar compuestos inorganicos de arsenico que son mas 
toxicos (Lillo, 2003). 

El arsenico, dependiendo de la forma en la que se encuentre, puede tardar varios anos 
en degradarse y su presencia en el agua se da en forma de arsenico pentavalente [As (V)] y 
arsenico trivalente [As (III)] principalmente (Lillo, 2003). 

En America del Sur, la presencia de volcanes y de una topograffa en la cual abundan 
acufferos hace que muchas de las fuentes de agua que son destinadas para el consumo 
humano se encuentren contaminadas con arsenico (Curto etal., 2002). 

En las aguas destinadas para el consumo humano, el arsenico generalmente proviene de 
pozos profundos donde existe pirita o arsenopirita (combinacion de hierro, azufre y arsenico), 
aunque muchas veces tambien proviene de fuentes de agua geotermal como sucede en la 
laguna de Papallacta, Ecuador. La norma de la Organizacion Mundial de la Salud permite una 
concentracion maxima de arsenico en aguas destinadas para el consumo humano de 10 pg/L, 
lo cual confirma que existe un grave problema de contaminacion en el agua que se toma de la 
laguna de Papallacta para la potabilizacion y posterior distribucion en la ciudad de Quito (INEN, 
2009). 

En Ecuador existe una norma sobre la calidad de agua de bebida adoptada en el ano de 
1992, en la cual se describe que la concentracion permisible de arsenico es de 0.05 mg/L 
(BVSA, 2004). La concentracion establecida en Ecuador por el INEN es de 10 pg/L. 

Segun Alvarez, 2003, en su investigacion de adsorcion de arsenico en suelos de cultivo 
de Castilla y Leon, Espana, el suelo que mejor adsorcion de arsenico tuvo fue el arcilloso, con 
9,23 pg/g Suelo, lo cual indica que las cianobacterias adsorben gran cantidad de arsenico en 
comparacion con estos suelos. Es importante tomar en cuenta que lo que se pretendfa probar 
era si el arsenico puede quedarse retenido en aquellos suelos destinados para cultivo, y la 
retention fue minima. El contar con biosorbentes para arsenico, como lo son algunas 
cianobacterias, abarataria costos de remediation y, en caso de que se contase con nuevas 
cepas, que tengan una capacidad de adsorcion mayor que en estudios anteriores, seria un 
aporte cientifico para solucionar problemas ambientales. 

El uso de la adsorcion selectiva de arsenico, utilizando materiales biologicos, oxidos 
minerales, carbon activado o resinas polimericas ha tenido un crecimiento acelerado. La 
utilization de carbon se lo hace desde hace mucho tiempo atras en la historia (Mohan and 
Pittman, 2007). El carboncillo era utilizado para la filtration de agua de bebida por los ancianos 
hindues en la India, y la comida carbonizada era un absorbente medico y un agente purificante 
en Egipto 1500 A.C (Perez et al., 2006). 

El crecimiento de cianobacterias en reactores no demanda de mucha inversion dado que 
se pueden usar recipientes de gran tamano para su crecimiento masivo y, ademas, las 
cantidades de reactivos utilizadas para la elaboration de medio de cultivo son minimas (Fay, 
1992). La energia necesaria para la production de cianobacterias es pequenisima y 
equivalente a la de cinco aireadores para pecera que son utilizados para airear un reactor de 
50 L, la energia luminica puede ser aprovechada a partir del sol con un fotoperiodo 
luz/oscuridad de 12/12 (Koch, 2007). 

Segun Devercelli y Garcia, 2001, se debe considerar toxicas a todas las cianobacterias 
inicialmente aisladas de algun ecosistema, hasta que no se demuestre lo contrario. El 50% de 
las cianobacterias identificadas hasta la actualidad tienen toxinas que pueden causar danos a 
seres vivos (Roset , Aguayo, Munoz, 2001), de alii la importancia de realizar pruebas de 
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toxicidad con distintos organismos hasta llegar a pruebas con humanos, ya que, algunas 
toxinas pueden causar danos en algunos organismos mientras que en otros no. 

El analisis de datos de toxicidad de cianobacterias usadas en este proyecto fue 
cualitativo y cuantitativo, se observaron y contaron pulgas muertas sometidas a distintas 
concentraciones de agua tratada por las cianobacterias (Brena, 2002). Se utilizo el programa 
Probit de Biostat para analizar dichos datos (Nunez y Hurtado, 2005), el cual arroja valores de 
Concentracion Letal Media de las sustancias analizadas con organismos vivos, al ingresar 
datos de los organismos muertos y el total de organismos utilizados en las pruebas de 
toxicidad, ademas de las concentraciones de la sustancia que se sospecha que es toxica 
(Barrios et al, 2007). 

2. MATERIALES Y METODOS 

2.1 Obtencion de las cianobacterias 

Las cianobacterias fueron proporcionadas por el Laboratorio de Microbiologia 
Ambiental del CEINCI-ESPE. 

2.2 Cultivo de cianobacterias en medio liquido 



Las cianobacterias fueron entregadas en medio BGIIo solido en cajas petri (Boune y 
Castenholz, 2001). A partir del cultivo solido, se establecieron los cultivos Ifquidos en medio 
BGIIo liquido (Koch et al., 2007). 

Para inocular las cianobacterias en medio liquido se las tomaron del medio solido con 
asa bacteriologica y se inoculo en medio liquido mediante la inmersion del asa con una leve 
agitacion de acuerdo al protocolo de siembra de cianobacterias propuesto por Koch et al., 
2007. Los frascos inoculados se mantuvieron con aireacion constante a un fotoperiodo 
luz/oscuridad de 16/8 con temperatura de 22±2 C y humedad de 36±5% durante 30 dias 
hasta el momento de la cosecha (Rippka et a/., 1979). 

2.3 Preparacion de agua contaminada con arsenico 



La preparacion del agua contaminada con arsenico se realizo con arseniato acido de 
sodio heptahidratado (Na 2 HAs0 4 .7H 2 0) segun Maeda et al., 1987. 

2.4 Recuperation de la biomasa 



Para poder realizar los ensayos de adsorcion de arsenico es importante contar con la 
biomasa limpia de medio de cultivo, con la menor cantidad de humedad posible para evitar 
interferencias en los ensayos de adsorcion. 

Se recupero la biomasa mediante centrifugacion a 4000 revoluciones por minuto (rpm) 
en tubos falcon de 50 mL durante 15 min. Posteriormente se la enjuago con agua destilada 
esteril durante 5 min para otra vez centrifugarla a 4000 rpm durante 15 min (Arunakumara, 
Xuecheng, Song, 2007). Para estandarizar la biomasa de los reactores, fue necesario hacer 
las diluciones necesarias y ajustar al mismo numero de cianobacterias en las diferentes 
repeticiones y tratamientos del diseno experimental llevado a cabo en esta investigacion. Las 
diluciones fueron realizadas al momento de iniciar los ensayos con el agua contaminada con 
arsenico. 
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2.5 Tratamientos aplicados a la biomasa 

2.5.1 Tratamiento de la biomasa mediante autoclave 



La biomasa fue resuspendida en agua destilada esteril para poder autoclavarla a 
121 C y 15 Psi de presion. El objetivo de este tratamiento fue lisar las celulas y obtener una 
mayor exposition de grupos funcionales presentes en las paredes celulares, para que exista 
una mejor adsorcion del arsenico. Despues del autoclavado, se centrifugo la biomasa a 4000 
rpm y se elimino el sobrenadante, para enjuagar la biomasa sobrante con agua esteril durante 
5 min y luego centrifugarla nuevamente para obtenerla limpia y lista para las pruebas de 
adsorcion (Ferrari, Guzman, Silva, Alcaraz, Silva, Gonzalez, 2004). 

2.5.2 Tratamiento de la biomasa mediante calor a 92C durante una hora 



La biomasa fue resuspendida en agua destilada esteril para someterla a calor. Se 
sumergieron los tubos falcon de 50 mL conteniendo la biomasa en un vaso de precipitation con 
agua destilada y nucleos de ebullicion, sobre una plancha de calentamiento y se mantuvieron 
60 min mientras el agua estuvo en ebullicion. Se centrifugo a 4000 rpm para eliminar el 
sobrenadante y luego se la enjuago con agua destilada esteril (modificado a partir del 
protocolo de Ferrari et al, 2004). 

2.5.3 Biomasa sin tratamiento 



La biomasa en los tubos falcon fue centrifugada a 4000 rpm se la resuspendio en agua 
esteril. Se elimino el agua mediante centrifugation a 4000 rpm para los ensayos con arsenico. 

2.6 Resuspension de la biomasa en agua contaminada con arsenico 



Se vertio 100 mL de agua contaminada con arsenico en frascos Boeco de 500 mL, se 
resuspendio la biomasa de cianobacterias y se colocaron los frascos en un agitador orbital 
SHO-2D WiseShake a 100 rpm (modificado a partir de Ferrari, 2004). 

2.7 Toma de muestras y cuantificacion de arsenico 



Se tomo una muestra de 5 mL cada 3, 9, 15, 21 y 30 min. Las muestras fueron filtradas 
para contar con agua libre de cianobacterias, y preservadas con acido clorhfdrico al 37% hasta 
alcanzar un pH de 2. La cuantificacion se realizo mediante generation de hidruros acoplado a 
un espectrofotometro de absorcion atomica (Modificado a partir de Torralba, 1996 y Varian, 
2004). 



2.8 Ensayo de toxicidad con Daphnia pulex 



Las pulgas de agua fueron facilitadas por el senor Daniel Orrico, estudiante de la 
Carrera de Ingenieria en Biotecnologia de la Escuela Politecnica del Ejercito. 
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2.8.1 Preparacion de agua reconstituida para Daphnia pulex 



Su preparacion se la realizo en un garrafon de 19 L de capacidad en agua destilada, 
para ello se adiciono 2,4 g de MgS0 4 , 3,84 g de NaHC0 3 y 0,16 g de KCI. Se agito hasta 
disolver completamente las sales. Paralelamente se disolvio 2,4 g de CaS0 4 .2H 2 0 en un litro 
de agua destilada; se agito durante 24 h en un agitador orbital para que se disuelva 
completamente. Al terminar de solubilizar el CaS0 4 .2H 2 0, se lo incorporo a los 19 L de la 
solucion de sales para completar un total de 20 L de agua reconstituida (Diaz eta. I, 2006). 

2.8.2 Neonatos utilizados para pruebas de toxicidad 



El dia anterior a realizarse las pruebas de toxicidad, se extrajeron con una pipeta de 
2500 pL todos los neonatos presentes en los estanques de crecimiento, para asi, al dia 
siguiente, poder contar con neonatos de entre 0 y 24 h que son los utilizados para pruebas de 
toxicidad universal (Diaz et al., 2006). 

2.8.3 Prueba de sensibilidad de Daphnia pulex con dicromato de potasio 



Las pruebas de sensibilidad en estos organismos se hacen para constatar su utilidad 
en pruebas de toxicidad, dado que si no tienen sensibilidad a los toxicos de referencia 
dificilmente tendran sensibilidad a posibles toxicos en el agua tratada con las cianobacterias 
(Diaz etal, 2006). 

Se preparo una solucion madre de dicromato de potasio de 100 mg/L en agua 
destilada, a partir de la cual se realizaron diluciones en agua reconstituida para obtener 
concentraciones de 10 mg/L, 1 mg/L 0,1 mg/L, 0,05 mg/L y 0.01 mg/L (Diaz et al, 2006). 

Se pusieron 10 neonatos de menos de 24 h en las diluciones de dicromato de potasio 
con el fin de conocer la concentracion letal media para este toxico de referencia. Se efectuaron 
dos repeticiones por cada concentracion y se mantuvo a los neonatos por 48 h en dichas 
soluciones. Ademas se preparo un bianco con agua reconstituida. Pasadas 48 h, se contaron 
los neonatos muertos, a los que se identified por falta de movimiento y en casos de dificil 
observacion a simple vista se observo la existencia de ritmo cardiaco mediante el 
estereomicroscopio (Diaz etal, 2006). 

Para las pruebas de toxicidad tanto con el agua tratada como con el medio de cultivo 
de las dos cianobacterias se llevaron a cabo diluciones en balones de 100 mL. De los 100 mL 
de cada dilucion se utilizaron 25 mL en cajas petri para cada repeticion con un total de tres 
repeticiones por cada dilucion. Por cada 25 mL de muestra se anadieron 10 neonatos menores 
a 24 h y se los mantuvo por 48 h, y se observaron al estereomicroscopio para contar los 
neonatos vivos y muertos (Diaz C et al., 2006). 

3. RESULTADOS Y DISCUSION 



3.1 Remocion de Arsenico 

Se obtuvieron cultivos masivos de cianobacterias para los ensayos de remocion de 
arsenico, se cosecharon las cianobacterias cuando la D0 75 o=1 • En los ensayos de remocion 
de arsenico, se obtuvieron porcentajes de remocion para cada pretratamiento de la biomasa y 
para el tipo de cianobacteria utilizada en los ensayos de remocion. La interaccion con el mejor 
porcentaje de remocion fue la del autoclave-consorcio-100 pg/L, con un 61.10% de remocion, 
le siguio la interaccion calor 92'C-consorcio-IOO pg/L con un 56.73% de remocion como se 
puede observar en la Figura 1 . 
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Los porcentajes mas bajos de remocion se obtuvieron cuando las cianobacterias no 
fueron sometidas a pretratamiento con porcentajes de remocion de 3.07% para la cepa 15 y de 
4.39% para el consorcio como se aprecia en la Figura 1 . 







Figura 1 Porcentajes de remocion de arsenico partiendo de una concentracion inicial de arsenico de 100 
pg/L con distinto pretratamiento y tipo de cianobacteria. CN=Consorcio, C15=Cepa 15. 



La remocion del arsenico se fue dando paulatinamente durante el tiempo de exposicion 
de las cianobacterias al metaloide. A los 30 minutos de exposicion, se pudo apreciar la 
disminucion de la concentracion de arsenico en el agua como se observa en la Figura 2. 




Figura 2 Disminucion del arsenico presente en el agua en el tiempo para una concentracion de partida 
de Arsenico de 100 pg/L. 



Se analizaron los datos en el programa SPSS 15 con la prueba de Tukey al 5%. Se 
pudo constatar que ningun pretratamiento es similar a otro. Se realizo un ANOVA en el que se 
indica que todos los factores utilizados en esta investigacion influyen en el porcentaje de 
remocion, ademas de las interacciones entre dichos factores. 



La cantidad de arsenico adsorbido se analizo mediante isotermas de adsorcion 
(Langmuir y Freudlich) realizando las regresiones lineales de las siguientes ecuaciones: 



Lineal: q = K*Ce 

Q max* K * Ce 

Langmuir: q = 



1 + k*Ce 
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Freundlich: q-K*Ce n 

Donde Ce es la concentracion en el equilibrio 
Qmax: adsorcion maxima 

Las isotermas de adsorcion realizadas para cada pretratamiento permitieron conocer 
cual adsorbe mayor cantidad de arsenico. El pretratamiento que obtuvo el Qmax mas alto de 
adsorcion fue el del autoclave a 121 C realizado al consorcio con 171.59 pgAs/gCianobacteria. 
Las isotermas se presentan en las Figuras 3 y 4. 




Figura 3 Isoterma de adsorcion para el consorcio tratado con autoclave utilizado en la remocion de 

arsenico. 




Figura 4 Isoterma de adsorcion de Langmuir para el consorcio tratado con autoclave utilizado en la 
remocion de arsenico. 

Posiblemente, la gruesa capa de peptidoglicano con la que cuentan las cianobacterias 
permitio la adsorcion del arsenico en su pared celular, por la presencia de grupos aminos en su 
estructura. No solamente la pared celular pudo interactuar en dicha adsorcion, sino que 
tambien organelos intracelulares con grupos aminos cargados positivamente. 

El mejor pretratamiento fue del autoclave, de igual manera el calor fue un 
pretratamiento eficiente, pero no supera al autoclave como se aprecia en la Tabla 1 . 
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Tabla 1 : Cantidad de Arsenico adsorbido en las cianobacterias con distintos tratamientos. 





Calor 


Autoclave 


SinTratamiento 


Cepa 15 


Consorcio 


Cepa 15 


Consorcio 


Cepa 15 


Consorcio 


Qmax 


101.829717 


107.628005 


123.660692 


171.594446 


1.694911428 


2.66883554 



3.2 Analisis de toxicidad con Daphnia pulex 

El ensayo de sensibilidad para Daphnia pulex con dicromato de potasio resulto en una 
DL50 = 0.604, lo cual permitio continuar con los ensayos de toxicidad del agua tratada con 
cianobacterias. 

Los ensayos de toxicidad realizados al agua tratada con cianobacterias dieron 
resultados negativos para muerte de Daphnia pulex, es decir no es toxica para este organismo. 



4. CONCLUSIONES 

El mejor pretratamiento realizado a la biomasa fue el del autoclave a 121 C y 15 psi de 
presion por obtener la mejor remocion de arsenico con un porcentaje de remocion de 61.10% 
en un lapso de tiempo de 30 minutos. 

La cepa 15 y el consorcio pretratados mostraron menor capacidad de adsorcion de 
arsenico en comparacion con investigaciones similares en las que se utiliza el quitosano como 
material adsorbente (Cumbal et al, 2006). Habrfa que analizar otras variables que podrian 
influir en la posible utilizacion de estas cianobacterias en escalados mayores tales como 
costos, facilidad de produccion, rapidez de produccion y otros. 

La cepa 15 y el consorcio que no se sometieron a pretratamiento remueven cantidades 
insignificantes del arsenico presente en el agua sintetica a nivel de laboratorio. 

La sensibilidad al dicromato de potasio de Daphnia pulex utilizada en esta investigacion 
dio positiva en la prueba de toxicidad aguda y permitio continuar con las pruebas de toxicidad 
en el agua de las cianobacterias. 

Tanto las cianobacterias como el agua tratada durante los ensayos de remocion de 
arsenico en esta investigacion no fueron toxicos para Daphnia pulex, al utilizar el metodo de 
toxicidad aguda con Daphnia magna de Diaz et al 2006. 

Los resultados obtenidos podrian promover una posterior investigacion a mayor escala 
en lo que respecta a remocion de arsenico en un biofiltro. 
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RESUMEN 

Contemplando criterios de resistencia y de desempeno sismico, se planted una tecnica 
inedita de diseno estructural, aplicable a las edificaciones de albanileria confinada de mediana 
altura (hasta de cinco pisos). Esta tecnica se basa en los resultados de multiples experimentos 
realizados en el Peru y en otros paises, asi como en una serie de estudios teoricos y en las 
ensenanzas dejadas por los sismos pasados que afectaron a este tipo de edificaciones. En el 
ano 2006, la tecnica quedo materializada en la actual Norma E.070 “Albanileria” del 
Reglamento Nacional de Edificaciones. 



Como derivacion de la tecnica indicada, se planted una tecnica similar para el diseno 
estructural de las viviendas de adobe confinado de hasta dos pisos, considerandose el 
comportamiento sismico de esta albanileria de baja calidad. Se debe resaltar que en el ano 
1998 se reforzaron viviendas reales de adobe de un piso, utilizando mallas de alambre 
electrosoldado cubiertas de mortero, que simulaban vigas y columnas de confinamiento; estas 
viviendas se comportaron exitosamente ante los sismos ocurridos en el sur del pais en los anos 
2001 y 2007. 



SUMMARY 

A new technique for structural design applicable to confined masonry buildings of medium 
height (up to five stories) was developed, including strength and seismic performance criteria. 
This technique is based on the results of multiple experiments performed in Peru and other 
countries, as well as in a series of theoretical studies and the lessons left by past earthquakes 
that affected this type of buildings. In year 2006, the technique was materialized in the current 
Masonry code (Norma E.070 in Spanish) of the National Building Code. 



As a derivation of the mentioned technique, a similar technique was developed for the structural 
design of confined adobe houses up to two stories high, considering the seismic behavior of this 
low quality masonry. It must be mentioned that in 1998 actual one-story adobe houses were 
reinforced, using welded wire meshes covered with mortar, which simulate confining beams and 
columns; these houses behaved successfully under the earthquakes that affected the southern 
part of the country in years 2001 and 2007. 
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1. PROPOSITOS DE LA INVESTIGACION E HIPOTESIS PLANTEADAS 

La albanilerfa confinada por elementos de concreto armado fue creada por ingenieros 
italianos, despues que el sismo ocurrido en 1908 en Messina, Sicilia, arrasara con las 
edificaciones de albanilerfa no reforzada. En nuestro pais, la albanilerfa confinada empezo a 
utilizarse despues de la ocurrencia del terremoto de 1940 en Lima, pero sin criterios de 
ingenierfa. Recien en el ano 1982, basandose en los resultados de escasos experimentos, se 
promulgo la primera Norma E-070 “Albanilerfa” [1], donde se empleaba criterios de diseno 
estructural por esfuerzos admisibles. Sin embargo, la fuerza sfsmica que se utilizaba estaba 
asociada a aceleraciones maximas del orden de 0.1 g en suelo duro, mientras que los sismos 
severos, de acuerdo a la Norma E.030 [3], pueden generar aceleraciones cuatro veces 
mayores (0.4g) para ese tipo de suelo, cifra que duplica al factor de seguridad (del orden de 
dos) que por fuerza cortante exigfa la Ref.1 para los muros confinados. Por consiguiente, es 
necesario que los muros de estas edificaciones sean disenados contemplando su incursion en 
el rango inelastico ante los sismos severos, previendo su falla por fuerza cortante, pero de tal 
modo que queden en un estado que sea economicamente reparable. 



1.1 Objetivos 

Los objetivos de la tecnica de diseno (Fig.1) son dos: 1) ante la accion de sismos 
moderados (o frecuentes, de hasta 0.2g de aceleracion en suelo duro), buscar que los muros 
se comporten elasticamente, sin que se agrieten; y, 2) ante los sismos severos (o raros, con 
aceleracion de hasta 0.4g en suelo duro [3]) que podrfan causar la falla por code en los muros, 
tratar que queden en un estado economicamente reparable, sin que disminuya su resistencia. 
Para los casos en que la edificacion se ubique sobre suelos blandos, las aceleraciones 
indicadas deben amplificarse por el factor “S” que especifica la Norma E.030 [3]. Debe 
resaltarse que por primera vez en el Peru, una tecnica de diseno estructural contempla los 
criterios de desempeno sfsmico mencionados. 



Fig.1 

Objetivos e 
hipotesis de la 
tecnica de 
diseno. 




1 .2 Hipotesis 

Previamente debe indicarse que en una estructura de albanilerfa confinada, los muros de 
albanilerfa proveen la resistencia y rigidez lateral, mientras que los confinamientos 
proporcionan la ductilidad. Por tanto, para lograr los objetivos indicados en el acapite 1.1, se 
adoptaron las siguientes hipotesis: 1) los confinamientos deben ser capaces de absorber la 
carga que produce el agrietamiento diagonal de la albanilerfa (VR en la Fig.1), para que la 
resistencia se mantenga estable durante la etapa inelastica; y, 2) debe existir una adecuada 
densidad de muros, medida a traves de la suma de sus resistencias a fuerza cortante (I VR), 
para dotar a la edificacion de la suficiente rigidez y resistencia, de tal modo que los muros 
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puedan soportar el sismo moderado sin mostrar danos, asf como quedar en un estado 
reparable despues del sismo severo. 

Para lograr el cumplimiento de las hipotesis senaladas, la estructura debe satisfacer lo 
siguiente: 



• Integracion de los Elementos Estructurales. Este tipo de edificacion se caracteriza por 
construirse primero el muro de albanilena, en segundo lugar se vacfa el concreto de los 
elementos verticales de confinamiento (columnas); y finalmente, se vacfa el concreto del 
techo en conjunto con el de las vigas, segun se muestra en la Fig.2. Esta secuencia 
constructiva produce una adecuada adherencia entre los materiales involucrados y un 
comportamiento integral de la estructura. 



Fig.2 

Secuencia en la 
construccion de 
la albanilena 
confinada. 





• Uso de Ladrillos Solidos. Puede emplearse ladrillos solidos de arcilla, sflico-calcareos o 
de concreto, que presenten hasta un 30% de perforaciones en su cara de asentado; de otro 
modo, la unidad calificarfa como hueca y terminarfa triturandose ante los sismos severos, 
lo que resulta peligroso. Las unidades huecas, a diferencia de las solidas, se trituran por las 
continuas aberturas y cerramientos de las grietas diagonales, perdiendose parte del area 
del muro y la resistencia en la etapa inelastica, segun se ha observado en diversos 
experimentos (Fig. 3). Cabe indicar que podrfa emplearse unidades huecas, solo si se 
llegase a demostrar que los muros se comportaran elasticamente (sin grietas diagonales) 
ante los sismos severos. 



Fig.3 

Trituracion de 
unidades huecas 
con perdida de 
area del muro, a 
diferencia de las 
unidades solidas. 




• Uso de Diafragmas Rigidos. Las losas aligeradas o macizas se comportan como laminas 
axialmente rfgidas para cargas contenidas en su piano, por lo que al estar integradas a los 
muros, uniformizan sus desplazamientos laterales ante los sismos, ademas de proveerles 
arriostramiento en sus bordes horizontales. Solo en el ultimo techo podrfa emplearse 
diafragmas flexibles (techo metalico o de madera, Fig.4), con la condicion de que existan 
vigas soleras de concreto armado que arriostren horizontalmente a los muros. 





166 



Revista CIENCIA Vol. 13 N.- 2 de 2010 




• Praxis Constructiva Adecuada. Experimentalmente se ha observado que defectos 
importantes en la construccion reducen sustancialmente la resistencia de los muros, tanto a 
carga vertical como a carga sfsmica, por lo que esos errores deben evitarse. Entre estos 
errores se tiene: grosores de juntas de mortero mayores que 15mm, cangrejeras en las 
columnas, secuencia constructiva inadecuada (Fig. 16), etc. 



2. METODOLOGIA EMPLEADA 



Para plantear la tecnica de diseno estructural, buscando cumplir con los objetivos 

especificados en el acapite 1.1 y usando las hipotesis indicadas en el acapite 1.2, se utilizaron 

procedimientos teoricos y experimentales. Los pasos que se siguieron fueron: 

• Evaluacion experimental de la resistencia a fuerza cortante de los muros (VR), 
empleandose materiales y mano de obra local. En estas pruebas ademas se determinaron: 
la deriva maxima para la cual el sistema queda en un estado reparable, el factor de 
reduccion de las fuerzas sfsmicas elasticas (“R = 3” en la Fig.1), la manera de predecir 
rigideces laterales, etc. 

• Analisis experimental de la influencia de diversos parametros estructurales sobre la 
resistencia al code de los muros (VR), tales como: la esbeltez, la carga de gravedad, etc. 

• Analisis experimental de la influencia de diversos parametros de la construccion que 
influyen en la resistencia al code de la albanileda (grosor de juntas, tipos de mortero, tipos 
de arena, etc.), mediante ensayos de probetas pequenas de albanileda (pilas y muretes, 
Fig. 12). 

• Formulacion de la resistencia al code, correlacionando los resultados de los muros a escala 
natural con aquellos obtenidos en el ensayo de probetas pequenas de albanileda. 

• Deteccion y analisis de las fallas usuales en las edificaciones de albanileda, despues de 
haber estado sujetas a sismos reales y simulados (experimentales). 

• Determinacion teorica de la densidad minima de muros requerida para cumplir con los 
objetivos indicados en 1.1, utilizando programas de analisis estructural inelastico paso a 
paso. 

• Formulacion de la tecnica de diseno y comprobacion mediante ensayos sfsmicos 
simulados. 
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3. TRABAJOS REALIZADOS y ANALISIS DE RESULTADOS 

La tecnica de diseno que se plantea se encuentra basada en: 1) los multiples ensayos 
estaticos y dinamicos realizados en el Laboratorio de Estructuras Antisismicas de la Pontificia 
Universidad Catolica del Peru y en el extranjero; 2) una serie de analisis teoricos; y, 3) las 
ensenanzas dejadas por los terremotos ocurridos en el Peru y en otros pafses [4], 

Los trabajos experimentales se iniciaron en el afio 1979 y continuan hasta la fecha. 
Estos trabajos se encuentran reportados en una serie de articulos, ilustraciones y videos que 
pueden ser vistos por Internet en el blog de “Investigaciones en Albanileria” 
http://bloa.pucp.edu.pe/albanileria [5]. En este documento solo se citaran los resultados de los 
proyectos mas relevantes. 



3.1 Ensayos de Carga Lateral Ciclica en Muros a Escala Natural 

Mediante ensayos de carga lateral ciclica con desplazamiento horizontal controlado (V-D, 
Fig. 5), se analizo la influencia de una serie de parametros sobre el comportamiento sismico de 
los muros de albanileria confinada, llegandose a obtener las conclusiones que se indican a 
continuacion. 




• La calidad de los ladrillos influye sustancialmente sobre la resistencia a fuerza cortante de 
los muros, asi como sobre su rigidez lateral (Fig. 6). No debe emplearse unidades huecas, 
porque se trituran disminuyendo drasticamente el area de code, la resistencia y la rigidez 
(Fig-7). 




Fig-6. Calidad de los ladrillos. 



Fig-7. Ladrillos huecos. 
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• La rigidez lateral elastica de los muros confinados, puede calcularse teoricamente 
transformando la seccion de las columnas de concreto en area equivalente de albanileria. 

• La deriva para la cual los muros quedan en un estado economicamente reparable, es y = 
0.005. Para esta deriva, el factor de reduccion de las fuerzas sfsmicas elasticas, calculado 
experimentalmente mediante el criterio de igualacion de energfas, puede estimarse en R = 
3. 

• Cuando la carga vertical genera en el muro esfuerzos axiales (c) mayores que 0.05 f'm, 
donde f'm es la resistencia a compresion axial de las pilas de albanileria (Fig. 12), se 
produce una reduccion importante de la resistencia inelastica a fuerza cortante del muro 
(Fig. 8) para derivas menores que 0.005; salvo que se le refuerce horizontalmente con 
varillas continuas ancladas en las columnas de confinamiento, en una cuantia nominal de 

0.001 (Fig. 9). Aun con este refuerzo horizontal, el esfuerzo axial “a” no debe ser mayor que 
0.15 f'm. 





• Las vigas que acoplan a los muros confinados (dinteles), deben ser preferentemente 
peraltadas (Fig. 10), para elevar la rigidez lateral y la resistencia a fuerza cortante. Estas 
vigas deben ser disenadas para soportar las cargas del sismo moderado, menores a las 
que provocan la falla por code del muro (sismo severo), tratando de que fallen por flexion 
en sus extremos, para que disipen energfa sfsmica antes de que se fracturen 
diagonalmente los muros. Es decir, ante los sismos, las vigas de acople deben funcionar 
como la primera Ifnea resistente. 
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Fig. 10 

Efectos del peralte de 
las vigas de acople. 




• La esbeltez de los muros, expresada a traves de la relacion M / (V L), donde “M” es el 
momento flector, “V” la fuerza cortante y “L” la longitud del muro, reduce la resistencia al 

code en un factor a = V L / M, con 1/3 < a< 1. Este efecto fue estudiado 
experimentalmente (Fig.11) y teoricamente mediante el metodo de elementos finitos, 
llegandose a la misma conclusion. Este hecho contribuye a explicar el porque las 
edificaciones de albanileria fallan por code en el primer piso (Fig. 15), ya que alii el 
momento flector es maximo. 




Fig.11. Efectos de la esbeltez del muro sobre la resistencia a fuerza cortante. 
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3.2 Formulacion de la Resistencia al Corte (VR) 

Los resultados de los diversos experimentos indicados en 3.1, fueron correlacionados 
con los obtenidos en los ensayos de las pilas a compresion axial y de los muretes a 
compresion diagonal (Fig. 12). La mejor correlacion se obtuvo con los muretes, ya que la forma 
de falla por corte en los muretes se refleja tambien en los muros a escala natural. Esto permitio 
obtener la ecuacion 1 . 

Ladrillos de arcilla y de concreto vibrado : VR - 0.5 v'm a t L + 0.23 Pg (1 ) 

Ladrillos sflico-calcareos: VR - 0.35 v'm a t L + 0.23 Pg 

En la ecuacion 1 , el termino “a” esta limitado a: 1/3 < a = VL/M <1; “t” es el espesor 
efectivo del muro, descontando el tarrajeo (salvo que este se aplique sobre una malla anclada 
al muro); “L” es la longitud total incluyendo las columnas; “Pg” es la carga de gravedad; y, v'm 
es la resistencia caracterfstica unitaria a corte puro de los muretes, obtenida como el valor 
promedio de la carga de rotura entre el area bruta de la diagonal cargada menos una 
desviacion estandar. 




Fig. 12. Ensayos de pilas a compresion axial y de muretes a compresion diagonal. 



3.3 Ensayos de Simulacion Sismica 



Mediante ensayos de especfmenes a escala reducida de 2 y 3 pisos (Fig. 13, Ref. 7), 
hechos en el simulador de sismos de la Universidad Catolica, se verificaron los resultados de 
los ensayos de carga lateral cfclica (acapite 3.1). Adicionalmente se obtuvo las siguientes 
conclusiones. 



Fig. 13 

Ensayos de 

simulacion 

sfsmica. 
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• La ecuacion 1 (acapite 3.2) para determinar la resistencia a corte quedo corroborada. En 
los pisos superiores al primero, donde la fuerza cortante sfsmica fue menor que la 
resistencia teorica “VR” correspondiente, no se produjo el agrietamiento diagonal de los 
muros (Fig. 13). 

• En sismos severos, los extremos de las columnas estan sujetos a fuertes compresiones por 
flexion (Fig. 13), por lo que sus estribos deben ser disenados para evitar el aplastamiento y 
el pandeo del refuerzo vertical. Al evitarse la falla por aplastamiento en el extremo de la 
columna, se elimina tambien la falla por compresion de la albanilerfa circundante. 

• En el rango elastico (sin grietas), la distribucion de las fuerzas sfsmicas en la altura de los 
especfmenes fue del tipo triangular (similar a la de la Norma E.030 [3]), pero al formarse la 
grieta diagonal en el primer piso, cambio a una distribucion practicamente uniforme. 

• La rigidez y la resistencia de los muros fueron similares a las obtenidas en los ensayos 
estaticos; sin embargo, la ductilidad decrecio en la prueba dinamica, debido al mayor 
deterioro que el ensayo dinamico produce tanto en la albanilerfa como en los 
confinamientos. 



3.4 Deteccion y Analisis de las Fallas Usuales 

• Falla por Corte. Tanto en los experimentos (figuras 13 y 14) como en los sismos reales 
(Fig. 15), se ha detectado que las edificaciones de albanilerfa fallan por fuerza cortante 
inicialmente en el primer piso. Esto se explica porque la fuerza cortante en el primer piso es 
mayor que la de los pisos superiores; y porque los efectos de esbeltez (Fig. 11), que 
reducen la resistencia al corte, son mayores en el primer piso, ya que allf el momento 
flector es maximo. Este hecho fue contemplado en la tecnica de disefio. 




Fig. 14. Experimento sfsmico japones en 
un edificio de 5 pisos de albanilerfa 
armada a escala natural [11]. 





Fig. 15. Falla por corte en el primer piso 
de la edificacion ante sismos reales. 



• Volcamiento. Se han reportado fallas de la albanilerfa en los ultimos pisos por cargas 
sfsmicas perpendiculares al piano del muro (Fig. 16). Esto se debio a errores cometidos en 
el proceso constructivo, donde primero se construyeron las columnas para despues 
levantar la albanilerfa, perdiendose la accion de arriostramiento vertical que deberfan 
proporcionar las columnas. Esta secuencia constructiva no integra a la albanilerfa con las 
columnas. 
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Fig. 16 

Pisco, 2007. 
Falla por 
volcamiento. 




• Falla por Cizalle. Este tipo de falla (Fig. 17) se ha detectado en las juntas de construccion 
columna-solera y columna-cimentacion. Los estribos en las columnas no controlan este tipo 
de falla horizontal, que debe ser evitada mediante refuerzo vertical y area de concreto. 

• Efectos del Alfeizar. El alfeizar de ventanas reduce la altura efectiva del muro 
rigidizandolo; lo cual hace que el muro absorba fuerza cortante en exceso, por lo que es 
conveniente aislarlo del alfeizar (Fig. 18). Tambien se recomienda este aislamiento porque 
el muro principal es portante de carga vertical, mientras que en el alfeizar la carga vertical 
es nula; esa diferencia de cargas causa desplazamientos relativos que muchas veces 
producen la fisuracion de la albahilena en la zona de encuentro, formandose una junta 
vertical que da lugar a que el alfeizar quede sin arriostre en ese borde ante cargas sfsmicas 
perpendiculares a su piano. 




• Falta de Confinamiento y de Arriostramiento. Debido a que el sismo genera fuerzas en 
los dos sentidos de la direccion longitudinal del muro estructural y en la direccion 
perpendicular a su piano, es necesario: 1) que la albahilena este completamente bordeada 
por los confinamientos; 2) que las columnas no esten espaciadas mas de dos veces la 
distancia que existe entre las vigas soleras; y, 3) que la relacion altura libre de la albahilena 
(distancia entre las soleras) dividida entre su espesor, sea menor que 20. Por ejemplo, en 
muros que presentan una sola columna, no se puede controlar el grosor de las grietas 
diagonales (Fig. 19), por lo que estos muros calificaran como “Muros No Confinados”. Por la 
gran degradacion de resistencia que tienen los “Muros No Confinados” despues de su falla 

por code, su aporte a la resistencia sfsmica del edificio (LVFt) en el rango inelastico debe 
considerarse nula. 
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3.5 Determinacion de la Densidad de Muros Confinados (Resistencia Global) 

Para demostrar la segunda hipotesis indicada en 1.2 (Fig.1) y asf lograr los objetivos de la 
tecnica, se analizaron teoricamente una serie de edificios estructurados por muros de 
albanilerfa confinada (uno de ellos aparece en la Ref.4). Se utilizo un programa de analisis 
inelastico paso a paso (TODA, Fig. 20), empleandose senales sfsmicas peruanas que fueron 
escaladas en sus aceleraciones para simular sismos moderados (0.2g en suelo duro) y severos 
(0.4g en suelo duro). 

La resistencia a fuerza cortante “VR” de cada muro confinado fue calculada con la 
ecuacion 1 (acapite 3.2), y en cada uno de los edificios se vario la densidad de muros, 
expresada a traves de la suma de sus resistencias a code (E VR), hasta cumplir con los 
objetivos indicados en 1.1. Cuando la densidad de muros fue excesiva, el comportamiento fue 
elastico, incluso ante el sismo severo; mientras que cuando la densidad de muros era baja, la 
demanda de ductilidad resulto muy alta, sobrepasandose la deriva de 0.005 para la cual el 
sistema es economicamente reparable. 

El analisis teorico determino el cumplimiento de los objetivos cuando en cada direccion del 
edificio se tuvo: E VR > V, donde “V” es la fuerza cortante especificada por la Norma E.030 [3] 
para la condicion de sismo severo (Fig.1). Por tanto, la resistencia global minima que debe 
proporcionarse a la edificacion en cada direccion y en cada piso es: E VR > V. 



Fig.20 

Analisis Sismico Inelastico 
Paso a Paso, mediante el 
programa de computacion 
TODA. 
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3.6 Determinacion de Fuerzas Internas en los Confinamientos del Piso 
Agrietado 

Producida la falla por fuerza cortante en los muros de un pano (con una fuerza “VR”, 
Fig. 21), las fuerzas internas en los confinamientos pueden calcularse por equilibrio. Se asume 
que las columnas y la solera no presentan momento flector, debido a que su curvatura por 
flexion es nula, porque los confinamientos continuan conectados a la albanilerfa incluso 
despues del agrietamiento diagonal (Fig.21). De este modo, las columnas trabajaran a traccion 
(T), a fuerza cortante (Vc) y a compresion (C), mientras que la solera trabajara a una traccion 
igual a “Vc”. 

Para el caso de muros de varios pisos y de un pano (Fig.21), las fuerzas internas en los 
confinamientos tambien pueden obtenerse por equilibrio. Para ello se asume que los pisos 
superiores al primero no se agrietan y que el cortante de agrietamiento se distribuye en cada 
nivel adoptando un patron de fuerzas laterales uniforme (f = VR / N, donde N es el numero de 
pisos). 




Fig.21 

Calculo por equilibrio de las 
fuerzas internas en los 
confinamientos de los muros 
de un pano. 




C 



f 



f 



Para el caso de muros con varios tramos (Fig.22), o muros conectados por vigas de 
concreto, o muros con paredes transversales, el problema es hiperestatico. De esta manera, se 
opto por solucionarlo mediante programas de computo, recurriendo a modelos matematicos 
como el indicado en la Fig.22, donde f = VR / N. Como conclusion se lograron las formulas 
simplificadas que aparecen en la Tabla 1. Estas formulas prestan atencion especial a las 
columnas localizadas en los dos extremos del muro, para que estas funcionen como topes que 
eviten el deslizamiento de la albanilerfa; de esta manera, las franjas de albanilerfa 
comprendidas entre las grietas diagonales, aportaran resistencia a carga lateral en la etapa 
inelastica, como se ilustra en la Fig. 23, descargando el trabajo de los confinamientos. 





Revista CIENCIA Vol. 13 N.- 2 de 2010 



175 



Fig.22 

Muro de varios 
pisos y varios 
tramos. Modelo 
analitico para 
hallar las fuerzas 
internas en los 
confinamientos. 





Tabla 1. Formulas para hallar las fuerzas internas en las columnas del primer piso. 



COLUMNA 


Vc (fuerza cortante) 


T (traccion) 


C (compresion) 


Interior 


V7?l Lm 
L (Nc + l) 


h 

VRl --Pc 

L 


VRl h 

Pc- 

2L 


Extrema 


VRl Lm 

1.5 

L (Nc + l) 


F - Pc 


Pc + F 



Fig.24 

llustracion de algunos 
parametros que 
intervienen en las 
formulas de la Tabla 
1. 




Los parametros que intervienen en las ecuaciones de la Tabla 1 (algunos de los cuales 
aparecen en la Fig.24) son: 
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VR1 = cortante de agrietamiento diagonal en el primer piso (calculado con la ecuacion 1) 
Lm = longitud del tramo mayor (LI o L2) > 0.5L. En muros de un tramo: Lm = L 

L = longitud total del muro, incluyendo las columnas de confinamiento 

h = altura del primer piso (o del piso agrietado en analisis) 

Nc = numero de columnas del muro en analisis. En muros de un tramo: Nc = 2. 

Pc = carga vertical tributaria en cada columna (incluye la carga “Pt” existente en un 
ancho tributario “B”, proveniente de muros transversales a la columna en analisis) 
M = Mul - V 2 VRI h 

Mul = momento flector en la base del primer piso producido por el sismo severo 
F = M / L = carga axial en las columnas extremas producida por el momento flector 

“M” 



4. RESULTADO: TECNICA DE DISENO. 

Como resultado de los trabajos indicados en el acapite 3, se elaboro la tecnica de diseno 
para las edificaciones de albanileria confinada. En esta tecnica se acepta que las edificaciones 
inevitablemente fallaran por fuerza cortante en su primer piso, cuando son sometidas a sismos 
severos, salvo que tengan una densidad de muros excesiva. A pesar de que la falla es por 
code, los experimentos (acapite 3) demostraron que es posible dotar de cierta ductilidad a los 
muros confinados, siempre y cuando los elementos de confinamiento sean capaces de soportar 
la fuerza asociada al agrietamiento diagonal de la albanileria “VR”. Cabe destacar que en el 
experimento japones [11] mostrado en la Fig. 14, los muros del primer piso fallaron por code y 
la degradacion de resistencia se produjo despues de superarse la deriva de 0.005 (adoptada 
como limite maximo en nuestra tecnica), producto de la trituracion de los extremos de los 
muros (Fig.25). 



Fig. 25 

Experimento japones en un 
edificio de albanileria 
armada de 5 pisos [11]. 
Falla por code en el primer 
piso y envolvente de los 
lazos histereticos V-D. 




El proceso de diseno planteado consta de cinco pasos: 1) verificacion de la densidad 
minima de muros confinados en las direcciones principales del edificio; 2) diseno por carga 
vertical; 3) analisis elastico ante el sismo moderado, evitando el agrietamiento diagonal de los 
muros; 4) determinacion de la resistencia al code “VR”; y, 5) diseno para sismos severos. 
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4.1 Paso 1. Verificacion de la Densidad Minima de Muros Confinados. 

Con el objetivo de evitar el colapso de la edificacion por deficiencia en la resistencia a 
carga lateral (Fig.26), en la etapa de predimensionamiento se debe proporcionar en cada 
direccion principal del edificio, una densidad minima de muros confinados, tal que se cumpla la 
expresion 2. Los muros restantes pueden carecer de confinamientos, e incluso ser portantes de 
carga vertical, ya que al estar conectados todos los muros al diafragma (losa de techo), su 
desplazamiento lateral sera controlado por el conjunto de muros estructurales confinados. 

£L t ZUS^ N (2) 

Ap ~ 56 

En la expresion 2, Z, U y S son los factores empleados en el calculo de “V” y se 
encuentran especificados en la Norma E.030 [3]; “N” es el numero de pisos del edificio; “L” es 
la longitud horizontal total del muro; “t” es su espesor descontado el tarrajeo, salvo que este se 
aplique sobre una malla conectada al muro; y, “Ap” es el area en planta del piso tfpico. 

En caso que no se cumpla la expresion 2 en alguna de las direcciones del edificio, se 
debe agregar muros de concreto armado (“placas”) o incrementar el espesor “t” en algunos 
muros de albahileria confinada. De emplearse placas (Fig.27), su espesor “t” en la expresion 2 
debera afectarse por la relacion de los modulos de elasticidad del concreto y de la albahileria 
(Ec/Em). 




Fig .26. Colapso de una vivienda con baja densidad 
de muros en la direccion de la fachada. Notese que 
los muros largos no contribuyen a la resistencia en la 
direccion transversal. 




Fig.27. Edificio con placa central. 



4.2 Paso 2. Diseno por Carga Vertical. 

Por cualquier metodo racional (metrado de cargas con areas tributarias, por ejemplo), se 
debe calcular la fuerza axial producida por las cargas de gravedad (Pg) en cada muro. Luego, 
el esfuerzo axial (a = Pg /(L t)) actuante debera ser menor que 0.15 f'm, donde f'm es la 
resistencia a compresion axial de las pilas de albahileria (Fig. 12); en caso contrario, habra que 
mejorar la calidad de la albahileria (f'm), aumentar el espesor del muro, o transformarlo en una 
placa. 

Si el esfuerzo axial (a) excede a 0.05f'm, debe agregarse en el muro una cuantia nominal 
de refuerzo horizontal (0.001), consistente en varillas continuas de un diametro no mayor que 
14”, ancladas en las columnas de confinamiento (Fig. 9). 

Una manera de reducir los esfuerzos axiales en los muros (a), es mediante el empleo de 
losas de techo armadas en dos direcciones (macizas o aligeradas), ya que estas distribuyen la 
carga proveniente del techo sobre los muros ubicados en ambas direcciones del edificio 
(Fig. 28). 
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4.3 Paso 3. Analisis Elastico ante Sismo Moderado y Diseno de otros 
Elementos. 

En esta tecnica se define al “sismo moderado” como aquel que produce fuerzas sfsmicas 
iguales a la mitad de las causadas por el sismo severo. La fuerza cortante en la base del 
edificio generada por el sismo severo, se determina utilizando la expresion dada por la norma 
sfsmica peruana E.030 [3]: V = ZUSCP/R, con R = 3. Bajo la definicion dada, para evaluar 
la fuerza cortante basal producida por el sismo moderado, puede emplearse la misma 
expresion, pero con R = 6. 

Es preferible que el analisis estructural sea hecho utilizando las fuerzas del sismo 
moderado, porque las fuerzas del sismo severo podrfan producir cortantes en los muros que 
superen a su resistencia al agrietamiento diagonal (VR); esto no significa que el muro haya 
colapsado, sino que incurrio en el regimen inelastico, generandose una redistribucion de 
fuerzas cortantes entre los diversos muros, que no es contemplado por los programas de 
analisis estructural elastico. 

El analisis sfsmico se debe realizar empleando cualquier procedimiento reconocido. En el 
modelo estructural se debe contemplar: 1) la accion de diafragma rfgido dada por las losas de 
techo; 2) los alfeizares de ventanas y tabiques no aislados de la estructura principal (*); 3) la 
participacion de los muros transversales a la direccion en analisis; y, 4) el hecho que los muros 
estan compuestos por dos materiales integrados: la albanilerfa y el concreto de las columnas; 
en este caso, el concreto puede transformarse en area equivalente de albanilerfa aplicando el 
criterio de la seccion transformada. En el analisis puede emplearse las siguientes expresiones 
empfricas que permiten obtener el modulo de elasticidad “Em” de la albanilerfa y su modulo de 
code “Gm”. 

• Em (ladrillos de arcilla) = 500 f'm 

• Em (ladrillos sflico calcareos) = 600 f'm 

• Em (ladrillos de concreto vibrado) = 700 f'm 

• Gm (ladrillos de cualquier material) = 0.4 Em ... (modulo de Poisson = 0.25) 

(*) Para simplificar el modelo matematico, es recomendable aislar los tabiques que presenter! 
discontinuidad vertical, asf como los alfeizares de ventanas (Fig. 18). 

El diseno estructural en la etapa de sismo moderado comprende lo siguiente: 

• Determinadas las fuerzas internas producidas por el sismo moderado en cada muro 
(momento flector “Me” y fuerza cortante “Ve”), debe verificarse que la fuerza cortante “Ve” 
no exceda al 55% de su resistencia al agrietamiento diagonal “VR”. La finalidad de lo 
indicado es que los muros se comporten elasticamente (sin grietas diagonales) ante el 
sismo moderado. 
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• Debe verificarse que las derivas inelasticas en cualquier piso no sean mayores que 0.005, 
de lo contrario habra que rigidizar al edificio. Estas derivas se calculan [3] como: 0.75 Ae R 
/ h, con R = 6 y Ae = desplazamiento elastico relativo entre pisos consecutivos. 

• Con la finalidad de que las vigas que conectan a los muros disipen energfa sfsmica antes 
que se agrieten los muros (Fig. 29), se plantea disenarlas en esta etapa a la rotura, de tal 
modo que su falla sea por flexion, tal como se indica en la Norma E.060 [6]. Para ello, debe 
emplearse un factor de amplificacion por efectos sfsmicos igual a 1.25. Una vez calculado 
el refuerzo longitudinal (As), se podra evaluar los momentos plasticos en sus extremos 
(Mp) y, posteriormente, por equilibrio, podra calcularse la fuerza cortante maxima (Vm) 
asociada al mecanismo de falla por flexion, y con ella disenar los estribos de confinamiento. 

Debe senalarse que estas vigas no colapsaran cuando actue el sismo severo (con fuerzas 
iguales al doble del sismo moderado), por la sobre resistencia que tienen, dada por: el 

factor de amplificacion (1.25), el factor de reduccion de resistencia a flexion ((j)= 0.9) y 
porque el diseno se realiza utilizando el esfuerzo de fluencia del acero (fy), que 
aproximadamente es 1.5 veces menor que la resistencia maxima del acero en su zona de 
endurecimiento. 



• De existir placas de concreto armado, tambien deben disenarse a la rotura en la etapa de 
sismo moderado, amplificando los esfuerzos elasticos por un factor de 1 .25, de tal modo 
que fallen por flexion, tal como se especifica en la Norma E.060 [6]. 

• La cimentacion tambien debe disenarse en la etapa de sismo moderado. El suelo de 
cimentacion no colapsara cuando ocurra el sismo severo, puesto que su resistencia 
admisible presenta un factor de seguridad del orden de tres, mayor que la relacion de 
fuerzas proporcionadas por el sismo severo y el moderado (dos). 

4.4 Paso 4. Determination de la Resistencia al Agrietamiento Diagonal “VR”. 

De acuerdo a lo indicado en el acapite 3.2, empleando la ecuacion 1 y utilizando las 
fuerzas internas provenientes del analisis sfsmico elastico (Me, Ve, acapite 4.3) para hallar el 
factor “a”, podra determinarse en cada muro de la edificacion su resistencia al agrietamiento 
diagonal “VR”. 

Ladrillos de arcilla y de concreto vibrado: VR - 0.5 v'm a t L + 0.23 Pg 




Fig.29. Diseno de las vigas de acople en la etapa de sismo moderado. 



( 1 ) 

Ladrillos silico-calcareos: 



VR = 0.35 v'm at L + 0.23 Pg 
1/3 < a = Ve L/Me <1 
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Para facilitar la labor del ingeniero estructural, se proporciona la Tabla 2, donde aparecen 
valores de la resistencia caracterfstica unitaria a code puro (v'm) de muretes construidos con 
varios tipos de ladrillo, asentados con mortero en proportion volumetrica cemento-arena 1 :4. 
Para otros tipos de unidades o de mortero, deberan efectuarse los ensayos de compresion 
diagonal (Fig. 12). 



Tabla 2. Resistencias caracteristicas de la albanileria. 



RESISTENCIAS CARACTERISTICAS DE LA ALBANILERIA MPa (kg / cm 2 ) 



Materia 


Denomination 


UNIDAD 


PI LAS 


MURETES 


Prima 


Comercial 


fb 


f 

J m 


v m 




King Kong Artesanal 


5.4 (55) 


3.4 (35) 


0.50 (5.1) 


Arcilla 


King Kong Industrial 


14.2 (145) 


6.4 (65) 


0.80 (8.1) 




Rejilla Industrial - INFES 


21.1 (215) 


8.3 (85) 


0.90 (9.2) 


Sflice-cal 


King Kong Normal 


15.7(160) 


10.8 (110) 


0.95 (9.7) 


Dedalo 


14.2 (145) 


9.3 (95) 


0.95 (9.7) 


Concreto 

vibrado 


KingKoncreto 


14.2 (145) 


13.2 (135) 


1.00 (10.2) 



4.5 Paso 5. Diseno para Sismos Severos. 

Este paso se compone de cinco partes: 

Paso 5.1: Verificacion de la resistencia global en cada piso del edificio (E VR > V) 

Con el valor de “VR” calculado para cada muro confinado, se obtiene en cada piso la 
suma de estos valores en cada direction principal del edificio (27 VR = resistencia global); la 
cual debe ser mayor que la fuerza cortante proporcionada por el sismo severo (V), especificada 
por la Norma E.030 [3]. Esta verificacion determina si la cantidad de muros confinados es 
suficiente como para que la estructura cumpla con los objetivos de la tecnica propuesta 
(acapite 1.1). En caso exista deficiencia de muros, podra anadirse “placas” de concreto 
armado, o aumentar el espesor “t” de algunos muros confinados (cambiando el aparejo, o con 
un tarrajeo aplicado sobre una malla anclada al muro), o mejorar la calidad de la albanileria 
(v'm), o incrementar el peralte de las vigas de acople para disminuir “Me” y, en consecuencia, 

aumentar “a” y “VR”. 

Si la resistencia proporcionada por los muros confinados (E VR) en el primer piso, 
incluyendo las placas de concreto armado en caso hubiesen, excede en mas de tres veces 
(factor “R = 3” de reduction de la fuerza sfsmica elastica) al cortante basal del sismo severo 
“V”, la estructura se comportara elasticamente ante el sismo severo (sin agrietarse 
diagonalmente). En tal situation, los muros requeriran confinamientos minimos que arriostren 
lateralmente a la albanileria ante cargas sfsmicas transversales al piano. Este es el caso de 
muros muy largos que se deforman en su piano principalmente por fuerza cortante, allf el 
refuerzo vertical de las columnas es rnmirno al no existir practicamente desplazamientos 
verticals por flexion que deformen al refuerzo. 

Paso 5.2: Fuerzas internas ante el sismo severo (Vu, Mu) y verificacion del 

agrietamiento diagonal de los muros localizados en los pisos superiores al 
primero (i > 1) 

La tecnica supone que ante la action del sismo severo, los muros del primer piso se 
agrietan diagonalmente, alcanzando su resistencia VR1. De este modo, las fuerzas internas 
elasticas (Me, Ve) se veran amplificadas en la relation VR1 / Vel, pasando a adoptar los 
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valores “Mu” y “Vu” asociados al mecanismo de falla por code en el primer piso (Fig.30). Este 
factor de amplificacion se calcula en el primer piso de cada muro y debe estar comprendido 
entre 2 y 3. El Ifmite “2” es la relacion de fuerzas del sismo severo y moderado, mientras que el 
Ifmite “3” es el factor “R” de reduccion de la fuerza sfsmica. El hecho de que el factor de 
amplificacion sea distinto en cada muro, indica una falla progresiva en los muros, en funcion a 
su nivel de resistencia “VR”. 



VR1 v 
Vel 



Fig.30 

Amplificacion 
de esfuerzos en 
un muro del 
edificio, 
agrietado en el 
primer piso. 




Esfuerzos 
ante el 

SISMO 

SEVERO 

Vui 

Mui 




Si en un piso “i” superior al primero (i > 1) de algun muro, se obtuviese que la fuerza cortante 
“Vui” supera a la resistencia al code “VRi”, entonces el muro de ese piso tambien se agrietara 
diagonalmente y sus confinamientos deberan ser disenados en forma similar al primer piso. 

Paso 5.3: Determination de las fuerzas internas en los confinamientos del primer 

piso 



Las fuerzas internas en los confinamientos de los muros del primer piso (T, C y Vc, ver la 
Fig. 21), se obtienen en forma aproximada mediante las formulas indicadas en la Tabla 1 del 
acapite 3.6, o utilizando modelos matematicos como el mostrado en la Fig.22. 



Paso 5.4: Diseno de los confinamientos del primer piso agrietado 

El diseno de las columnas de confinamiento se realiza empleando las formulas 
especificadas en la Norma de Concreto Armado E.060 [6], sin amplificar los valores de rotura 
de la Tabla 1 : T, C y Vc, pero contemplando los factores de reduccion de resistencia (c» de la 
Norma E.060. La resistencia a compresion del concreto (f'c) debera ser mayor que 175 kg/cm 2 . 

El area del refuerzo vertical de las columnas se determina disenandolo por corte-friccion 
(Vc) y traccion (T) combinada, ya que ambos efectos actuan en simultaneo, como se muestra 
en la Fig.21. Por lo menos debera emplearse un area de acero Asm = 0.1 f'c Ac / fy, donde 
“Ac” es el area de la columna, o 4 varillas con diametro de 8mm, lo que sea mayor, formando 
una canastilla. El area minima de acero indicada (Asm), es la requerida para absorber la carga 
que produce la fisuracion del concreto de la columna en traccion por flexion del muro. 

El area del nucleo de concreto (confinado por los estribos, Fig.31) se determina 
disenandolo por aplastamiento con una compresion igual a “C”, como si se tratase de una 
columna sin esbeltez. Luego, agregando el recubrimiento, el area total de la columna no debe 
ser menor que la obtenida en su diseno por corte-friccion (Vc), o que 15 t (cm 2 ), donde “t” es el 
espesor del muro. 

Los estribos se disenan para evitar la falla por aplastamiento del nucleo de concreto 
(Fig.31), debiendose colocar por lo menos estribos de diametro Vi”: 1 @ 5, 4 @ 10, r @ 20cm. 

La zona a confinar con estribos a corto espaciamiento (extremos de las columnas, 
Fig.32), abarca una longitud de 1.5 veces el peralte de la columna o 45cm, lo que sea mayor. 
Debe indicarse que los estribos no requieren disenarse por code, ya que la albahileria aun 
fragmentada aporta resistencia lateral (Fig. 23), ademas del aporte dado por el refuerzo 
horizontal en caso existiese; asimismo, los estribos son ineficientes cuando la falla en las 
columnas es por corte-friccion. 
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Fig.31 

Falla por aplastamiento 
de una columna 
(izquierda) y 
desprendimiento del 
recubrimiento, 
manteniendo el nucleo 
confinado sin danos 
(derecha). 





10 cm 



[] adicionales 
@ 10 cm 



t 



n i nr-^ i i~ 


i i l l l n 




Min. 45 cm o 1.5 d 

[] de confinamiento , o 

Junta rayada 
^ humedecida 




Sobrecimiento 

1:8 + 30% de piedra mediana (3") 



2[]s 

Adicionales 

45 cm 
o 

15 d 



Cimiento 

[] Montaje Concreto Ciclopeo 

- 1:10 + 30% de piedra grande (10") 



[] <t> 1/4", 1 @5, 

4 @ 1 0, r @ 25 cm 



Fig.32. Disposicion minima de estribos en las columnas de confinamiento. 



El refuerzo longitudinal de las vigas soleras, encargadas de transmitir las cargas sismicas 
desde la losa de techo hacia la albanileria, se disena a traccion (Ts = Vfe VR1 Lm / L), 
agregando estribos minimos de diametro 14 ”: 1 @ 5, 4 @ 10, r @ 20cm. A diferencia de las 
vigas dinteles que se deforman por flexion, las soleras no necesitan disenarse por fuerza 
cortante ni por flexion, debido a la gran seccion transversal que existe encima del primer piso 
(Fig.33); asimismo, como las soleras siguen la deformacion del muro (de acuerdo a la hipotesis 
de Navier, la seccion transversal del muro no se curva por flexion), requeriran area minima en 
su seccion transversal. 




Paso 5.5: Diseno de los confinamientos de los pisos superiores no agrietados 

En los casos en que “Vui” sea menor que “VRi”, para i > 1, la albanileria absorbe la 
fuerza cortante sin agrietarse diagonalmente, por lo que el muro no requiere refuerzo horizontal 
y las columnas no necesitan disenarse por corte-friccion. En estos casos, es suficiente disenar 
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a las columnas extremas del muro por traccion (T) y aplastamiento por compresion (C), 
producidas por “Pc” y el momento flector Mui = Mei VR1 / Vel (ver en la Fig. 24: F = Mui / L). 

Las columnas internas de estos muros no necesitan disenarse ante acciones coplanares, 
puesto que siguen la deformada de la albanileria no agrietada. Sin embargo, los 
confinamientos deberan ser capaces de funcionar esta vez como arriostres de la albanileria 
ante acciones sismicas perpendiculares a su piano; para ello, el espaciamiento maximo entre 
las columnas no debera ser mayor que dos veces la distancia que existe entre las vigas 
soleras, y la relacion altura libre de la albanileria entre su espesor debera ser menor que 20. 



5. CONCLUSIONES 

• La tecnica de diseno sismico expuesta, ha sido comprobada con exito en experimentos 
estaticos y dinamicos realizados en muros de albanileria confinada a escala natural y 
reducida. Esta tecnica es inedita a nivel mundial y por primera vez en el Peru, aplica 
conceptos de desempeno de la estructura ante sismos moderados y severos. 

• La tecnica de diseno expuesta en este documento, ha sido incorporada en la Norma 
Tecnica de Edificaciones E.070 “Albanileria” del 2006 [2], para el diseno sismico de las 
edificaciones estructuradas con muros de albanileria confinada. 

• La falla por code que se desarrolla en las edificaciones de mediana altura estructuradas 
con muchos muros, debe tomarse con naturalidad. Aunque esta falla sea fragil, su efecto 
no es negativo, siempre y cuando se adopten las siguientes precauciones: 1) proporcionar 
una adecuada densidad de muros para limitar los desplazamientos laterales inelasticos; 2) 
disenar a las vigas de acople como una primera linea resistente, para que disipen energia 
sismica antes que se fracturen los muros; 3) adicionar refuerzo horizontal cuando el 
esfuerzo axial sea excesivo; y, 4) disenar a las columnas de confinamiento como ultima 
linea resistente. 

• Una manera de producir la falla por flexion en los muros es mediante la inclusion de 
refuerzo horizontal importante. Sin embargo, la forma que tienen los ladrillos, sin canales 
que permitan alojar a las varillas horizontales, hacen que como maximo se pueda utilizar 
varillas de diametro !A” embutidas en las juntas de mortero; de emplearse varillas de mayor 
diametro, el grosor de las juntas horizontales se incrementaria, lo que disminuiria 
significativamente la resistencia al agrietamiento diagonal. Cabe resaltar que los 
experimentos hechos en muros donde se uso el doble de la cuantia nominal de refuerzo 
horizontal especificada en esta tecnica (0.001), proporcionaron un incremento de la 
ductilidad al reducirse el deterioro de la albanileria, pero solo incrementaron la resistencia 
maxima en 10%. Asimismo, en muros no agrietados diagonalmente, este refuerzo no 
trabaja al trasladarse en conjunto con el muro. 

• Las vigas que conectan coplanarmente a los muros (dinteles) juegan un papel importante. 
Cuando estas vigas tienen un peralte razonable (menor que 60cm), rigidizan a los muros y 
retardan su falla por fuerza cortante. Asimismo, disminuyen el momento flector en la base 
del muro, lo que trae consigo: 1) un incremento de la resistencia al code del muro al 
disminuir los efectos de esbeltez; 2) una reduccion en el tamano de la cimentacion; 3) una 
disminucion del refuerzo vertical a colocar en las columnas de confinamiento; y, 4) una 
reduccion de la compresion en las columnas extremas (causa de la falla por 
aplastamiento). Ademas, las vigas peraltadas limitan el giro por flexion de los muros, con lo 
cual, decrece la posibilidad de una falla por punzonamiento de la losa de techo, causada 
por la rotacion de la seccion del muro. 
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6. POSIBILIDADES DE APLICACION 



6.1 Aplicacion en la Norma E.070 “Albanileria” y en Codigos Extranjeros 

Desde el ano 2006, la tecnica de diseno expuesta se aplica oficialmente en el Peru, al 
habersela incorporado en la Norma E.070 “Albanileria” [2] del Reglamento Nacional de 
Edificaciones. 

Tambien es posible que los criterios vertidos en este documento puedan aplicarse en 
reglamentos de otros paises, aunque los valores numericos deberan adecuarse a la realidad de 
cada pais. En este sentido, debe destacarse que el articulo donde se describe esta tecnica: 
“Design Proposal of Confined Masonry Building”, obtuvo una mencion honrosa en el evento 
Tenth North American Masonry Conference, realizado en junio del 2007 en St. Louis, Missouri, 
USA, bajo la organizacion de The Masonry Society y la University of Missouri-Ftolla [8]. 



6.2 Aplicacion en el Diseno del Adobe Confinado 

Pensandose en los casi 10 millones de peruanos que habitan en viviendas de adobe, 
vulnerables ante los sismos, los conceptos vertidos en la tecnica expuesta fueron aplicados al 
caso del adobe confinado para viviendas de hasta dos pisos. Para ello, se trato al adobe como 
una albanileria de baja calidad (la resistencia unitaria al code es del orden de 0.5 kg/cm 2 ) y se 
utilizaron confinamientos de concreto pobre (f'c = 100 kg/cm 2 ) para abaratar costos. En este 
caso, las columnas de confinamiento se disenan solo por corte-friccion, ya que las tracciones y 
compresiones son minimas, mientras que la viga solera debe disenarse para servir de arriostre 
horizontal a los muros, puesto que los diafragmas son flexibles (techos de madera, Fig. 34). 




Bajo los conceptos indicados se diseno un modulo de dos pisos (Fig. 35). En el primer 
piso se adiciono refuerzo horizontal de Vi” cada seis hiladas (para proveer integracion entre 
todos los adobes y las columnas), sobre una capa de mortero cemento-arena 1 :5 para evitar su 
corrosion. Este modulo se ensayo ante sismos simulados en mesa vibradora, obteniendose 
resultados exitosos [9] que se utilizaron para refinar la propuesta de diseno sismico que 
aparece en la Ref.5. 

Cabe destacar que el articulo “Seismic Behavior of a Two Storey Model of Confined 
Adobe Masonry” [9], sobre el comportamiento sismico del modulo mencionado, obtuvo el 
primer puesto en el simposio "Eleventh Canadian Masonry Symposium", realizado en Toronto, 
Canada, en junio del 2009, bajo la organizacion de Me Master University y Canada Masonry 
Design Centre. 
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Fig. 35. Confinamientos y refuerzo horizontal en el primer piso del modulo de adobe ensayado en 
mesa vibradora. Una de las conclusiones fue la necesidad de reforzar horizontalmente tambien al 
piso 2, para integrar todos los adobes con las columnas. 



6.3 Aplicacion en las Viviendas Existentes de Adobe No Reforzado de un 
Piso 



En las viviendas de adobe existentes, carentes de refuerzo, los confinamientos (vigas y 
columnas) se simulan a traves de franjas de mallas de alambre electrosoldado de 1mm de 
grosor, con cocadas de % de pulgada (Fig. 36). Estas mallas se recubren con una capa de 
mortero 1:5 (cemento-arena), a fin de protegerlas de la intemperie y proveer rigidez, 
especialmente en la zona central superior de los muros, ya que los diafragmas (usualmente 
compuestos por vigas de madera, cafia y estera con torta de barro) son flexibles. 



Estas mallas se colocan en las dos caras del muro y se interconectan entre sf mediante 
alambre #8. Previamente, se hacen perforaciones de 5x5cm cada 50cm en el muro, donde se 
instala el alambre #8, para despues taponarlas con mortero 1 :5. Luego la malla se clava a la 
pared usando chapas, o laminas de metal o de madera en la cabeza del clavo (Fig. 36 y Ref.5). 
La malla que Simula a la columna no necesita conectarse a la cimentacion ni al techo (no se 
requiere desmontar el techo existente), porque los muros se deforman principalmente por 
fuerza cortante, desarrollando en sus extremos tracciones y compresiones rmnimas por flexion. 




En el ano 1998 (Fig.36), veintitres viviendas existentes ubicadas en distintas zonas 
sfsmicas del Peru, fueron reforzadas bajo el sistema descrito, en el marco del proyecto GTZ- 
CERESIS-PUCP [12], Seis de estas viviendas, ubicadas en Tacna y Moquegua, soportaron sin 
ningun dano el terremoto ocurrido el ano 2001 (Fig. 37). Otras cinco viviendas, ubicadas en lea, 
soportaron el terremoto de Pisco del 2007 [10], sin ningun dano (Fig.38); debiendose resaltar 
que una de ellas es la estacion sfsmica “Parcona” del Instituto Geoffsico del Peru (IGP), donde 
se registro una aceleracion maxima de 0.5g en suelo blando. En ambos terremotos, las 




186 



Revista CIENCIA Vol. 13 N.- 2 de 2010 



viviendas de adobe no reforzadas, vecinas a las reforzadas, sufrieron fuertes danos o 
colapsaron (figuras 37 y 38). 

Este comportamiento exitoso trascendio internacionalmente, logrando el apoyo de 
Alemania e Italia para la construccion de mas de 500 viviendas de adobe, que se asignaron a 
las familias afectadas por el terremoto del 2001 en varios pueblos de la zona altoandina de 
Arequipa [12]. Debe mencionarse que esta vez, por tratarse de viviendas nuevas, se utilizaron 
vigas soleras de concreto pobre, adicionales a las mallas electrosoldadas que continuaron 
simulando a los confinamientos verticales (figuras 39 y 40). Este nuevo sistema de refuerzo fue 
desarrollado por los autores del presente documento, quienes ademas supervisaron la 
construccion respectiva. 




Fig .37. Viviendas reforzadas y no reforzadas en Moquegua despues del sismo del 2001 . 




Fig.38. Viviendas reforzadas y no reforzadas en lea despues del sismo del 2007. 




Fig.39. Vivienda en Construccion. 




Fig.40. Vivienda terminada 


Notese la solera de concreto con 




en la zona altoandina de 


dientes en los extremos y la malla. 




Arequipa. 
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RESUMEN 

Para poder resolver en forma manual y comparar los resultados con programas de 
computacion, se realiza el analisis estatico de una columna de suelo de un metro de ancho y 
cinco metros de alto, la misma que esta confinada por un lado y sometida al empuje del agua. 

Se consideran dos modelos numericos de analisis, en el primero se trabaja con un 
elemento finito rectangular Q4, que considera dos grados de libertad por nodo y en el segundo 
se trabaja con un elemento finito mixto rectangular, para suelo saturado, Q8, que tiene tres 
grados de libertad por nodo, dos desplazamientos y una presion. Se comparan los resultados 
obtenidos con estos dos modelos en terminos de desplazamientos. 

Finalmente en base al ejemplo desarrollado, se presents una manual de como utilizar 
el programa GID, en sus dos fases como un Pre y Postproceso. Se describe tambien un 
programs desarrollado en Matlab para resolver el problems de elementos finitos considerando 
suelo saturado. 



ABSTRACT 

In this article the static analysis of soil column is carried out in order to compare a 
manual solution against a numerical procedure based on the finite element method for soil 
mechanic. The soil column is about 1m wide and 5m height. 

Two different numerical analysis are consider, the first one use a bilineal quadrilateral 
finite element Q4 with two degree of freedom per node, while, the second method use a 
quadratic quadrilateral finite element Q8 with three degree of freedom per node (two 
displacement and the pore pressure). The results obtained are compared in terms of 
displacements freedom degree. 
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Finally based on the developed example, a short manual is presented in order to 
explain the software GID as Pre and Postprocessor, and also a Matlab subroutine is described 
in order to solve the coupled consolidation problem of saturated soil by finite elements method. 



1. INTRODUCCION 

Si bien es cierto en este articulo se realiza el analisis estatico de un problema 
elemental, como lo es una columna de suelo confinado sometida a cargas hidrostaticas, no es 
menos cierto que se puede extender la aplicacion a problemas tan complejos como es el 
analisis de la Presa la Esperanza, que esta ubicada en la Parroquia Quiroga, aproximadamente 
a 12 km., de Calceta, en la Provincia de Manabf, construida sobre el rfo Carrizal. En un futuro 
cercano, se realizara el analisis estatico de esta presa que se muestra en la figura 1 . 




Figura 1 Vista aerea de la Presa la Esperanza. 



La Presa la Esperanza, almacena 45 millones de m 3 de agua y sirve para atender a 
25000 Ha., fue disehada en la decada de los 70 del siglo pasado. Un code transversal en el 
centra de la misma, se indica en la figura 2. 

Dentro de un cfrculo, en la figura 2, se ha identificado los materiales que conforman la 
Presa, son 7 en total, los mismos que se describen en la tabla 1 . La presa es de materiales 
sueltos con nucleo central de arcilla y espaldones de arenisca; su altura es de 47 m. sobre el 
cauce y 696 m. la longitud de su coronacion. Los taludes son 3.5:1 aguas arribas y 2.75:1 
aguas abajo. (Intecsa y Geosisa 1982), (Dragados y Geosisa 1993). 

Aguas arriba del eje de la presa se construyo una Atagufa que se utilizo para el desvfo 
del rfo durante la construccion; posteriormente quedo integrada en el cuerpo de la propia presa. 
En el contacto del talud aguas abajo del nucleo de la presa y del espaldon de arenisca, se 
intercalo un dren de gravas graduadas que al llegar a la base de cimiento continua 
horizontalmente hasta el pie de aguas abajo. El talud de aguas arriba se protegio con escollera 
y el de aguas abajo con pasto para impedir la erosion producida por la lluvia. 
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SECCION TIPO 
PRESA LA ESPERANZA 

EJE REFERENCIA (R=2093.50m) EJE REFERENQA PRESA (R=2000.00m) 

DIQUE DE PROTECCION 




Sobre el antiguo cauce del no se construyo bajo la presa una pantalla de 
impermeabilizacion de Hormigon Plastico para evitar las fugas de agua a traves de los 
materiales que rellenan el cauce, y en las laderas se efectuaron cortinas de inyeccion de 
lechada de cemento. 

La pantalla de impermeabilizacion de hormigon plastico es de 210 m. de longitud y 1 m 
de espesor, atraviesa los 60 m de profundidad de relleno hasta empotrarse en el fondo rocoso. 
Las arenas que sirven de cimentacion al espaldon de aguas abajo se consolidaron mediante 
compactacion dinamica para evitar el hundimiento y ruptura en la presa por licuefaccion de las 
mismas. Los tipos de inyecciones que se aplicaron en las laderas fueron: 

• Cortina trilineal con tubo-mangito y convencional, y pantalla bilineal de Jet Grouting 
sobre la margen derecha. 

• Pantalla monolineal y bilineal de Jet Grouting en la margen izquierda. 



Se ha descrito en forma rapida, la futura aplicacion del metodo de los elementos finitos. 
Ahora, en este artfculo dos son los objetivos principales que se desean presentar. El primero 
de ellos es comparar la respuesta de una columna de suelo confinada ante cargas hidrostatica 
del agua, empleando dos modelos de elementos finitos, el Q4 que considera dos grados de 
libertad por nudo y el suelo completamente seco. Aguiar (2010), y el elemento finito Q8 para 
suelos saturados Mroginski (2008), Aguiar y Mroginski (2010), para los dos casos la carga 
lateral del agua sobre la columna de suelo se aplica en un solo escalon de carga. 

El segundo objetivo es mostrar como se utiliza el programa GID (Sistema Grafico 
Interactive), desarrollado en el Centro Internacional de Metodos Numericos, CINME, de la 
Universidad Politecnica de de Cataluna, debido a que es una herramienta muy poderosa para 
utilizarlo como Pre y Postproceso, en combinacion con otros programas de elementos finitos. 
Zarate etal. (2001). 

En este artfculo se ilustra la forma de uso del GID como Pre Proceso, para: definir la 
geometrfa de la columna de suelo confinado, generar la malla de elementos finitos a utilizar, las 
condiciones de borde de la malla, las cargas actuantes, el tipo de material de la malla, la forma 
de aplicacion de la carga hidrostatica. En esta modalidad se generan varios archivos los 
mismos que deben ser abiertos y utilizados por otros programas de computacion que resuelvan 
el problema de elementos finitos, en este caso los desarrollados por Mroginski (2008). 

Luego se describe el programa de Mroginski (2008), con el cual se leen los datos del 
GID, se obtiene la matriz de rigidez de la estructura y se encuentran los desplazamientos para 
cada escalon de carga y tambien se indica la forma como se graban los resultados para luego 
ser utilizados por el programa GID como Postproceso. 
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Tabla 1 Materials que conforman la Presa la Esperanza. Dragados y Construcciones (1994). 



Zona 


Funcion 


Caracteristicas 

materiales 


Especificaciones 

compactacion 


Yacimiento 

recomendado 


Volumen 
aprox. (m3) 


1 


Nucleo 


Arena arcillosa tamano 
maximo:10cm. Porcentaje 
que pasa: on.- 4: 70-10 
on.- 200: >30 
(LI <60%:; Ip: 7-35%) 


Rodillo pata de cabra o con 
compactadores de 10T 
esp.capa:en funcion terraplen 
de ensayo control de 
humedad +2 a +4 h.op. 
Densidad minima 98% d.m.p 


Ma-2{cobertura) 


1 '980.000 


2 


Filtro 


Arena fina y limpia porcentaje 
que pasa: sn.-1 0:100 

on. -40: 90-100 on.-100:5-25 
on. -200: 0-5 


Rodillo vibratorio liso 5T 
espesor capas: 25cm riego 
200 l/m3 densidad relativa 
minima:60% 


Arena de San 
Jacinto 


6.000 


3 


Dren 


Gravas y arenas g radii ad as 
porcentaje que pasa: 
on. -3/4": 75-1 00 on.4: 40-70 
on.-IO: 25-5 on. -40: 0-15 
on. -200: 0-5 


Rodillo vibratorio liso 5T 
espesor capas: 25cm riego 
200l/m3 densidad relativa 
minima:60% 


Cantera de 
basalto 


113.000 


4 


Espalddn 


Grava arenosa con finos 
tamano maximo:20cm 
porcentaje que pasa: 
on. 4: 0-70 on. -40: 0-60 
on. -200:0-30 


Rodillo pata de cabra o con 
compactadores de 10T 
esp.capa:en funtidn terraplen 
de ensayo control de 
humedad +2 a +4 KGp. 
Densidad minima 98% D.M.P 


Productos 

excavacion 


1715.000 


5 


Protection 
t a 1 lid aguas 
arriba 


Arenisca calcarea compacta 
porcentaje carbonato >20% 
densidad >1.5 tamano 
m axi mo :75c ni tamano 
medio::30cm 


Rodillo vibratorio liso 10T 
espesor colocation: 100cm 
riego no se precisa 


Ma-2 (arenisca 
calcarea) 


110.000 


5a 


Proteccidn 
talud aguas 
arriba 

seleccionado 


Roca basaltica dura desgaste 
los angeles < 35% ataque 
sulfatos (2.5)< 10% tamano 
maximo: 75cm tamano 
medio: 30cm 


Rodillo vibratorio liso 10T 
espesor colocation: 100cm 
riego no se precisa 


Cantera de 
basalto 


42.000 


6 


Rellenos 

estabilizacion 

quebradas 

ladera 

izquierda 


Producto de excavacionftodo 
uno) tamano maximo::60cm 
materia organica <5%(peso) 


Rodillo pata de cabra o con 
compactadores espesor 
capas:50cm numero de 
pasadas:2 


Productos 

excavacion 


152.000 


7 


Relleno pie 
drenante 


Bolones. grava y arena 
tamano maximo: 10cm 
porcentaje que pasa 3/4pulg.: 
30-80 n.4: <20 


Rodillo vibratorio liso 10T 
espesor colocacion: 50cm 
numero pasadasA riego: 200 
I/m 3 


Cantera de 
basalto 


27.000 



Para la etapa de Postproceso se indica la forma de uso para presentar en forma grafica 
los desplazamientos encontrados. Todo esto relacionado a la columna de suelo confinada 
sometida a carga hidrostatica. 
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2. ELEMENTOS FINITOS Q4 Y Q8 

En la figura 3 se presentan los dos elementos finitos que se van a utilizar en el analisis 
estatico, se trata de dos elementos de clase C 0 en que para un nodo determinado se garantiza 

la continuidad de la funcion; los elementos de clase C, para un nodo determinado garantizan 
la continuidad de la funcion y su primera derivada. En general se dice que un elemento es de 
clase C m si su campo de desplazamientos tiene continuas las m primeras derivadas. Ofiate 
(1992). 



k 

7 

^ 


k 8 


k 6 

> < 13. 


I 14 ft 1 


k 10 

^ o 


F* 

i fcl 




i 

15^ 

, 3 . U, 


p 7 P 6=0 

k 16 j 

1 P8=0 ' 

ft 


P 5 

k 8 

*— ► 7 

P 4=0 

, 5 fc . 


F 1 

i 


L k 

2 

Q4 


, * *1 

4 Pi 


1 ^ : 

2 P2=0 

Qs 


^ P3 



Figura 3 Descripcion de los elementos finitos a utilizar. 



Para el elemento finito Q4 se trabaja con funciones de forma lineales, en cambio para 
el elemento finito Q8 con funciones de forma parabolicas. En Aguiar y Mroginski (2010) se 
indican las funciones de forma en coordenadas naturales para los elementos finitos Q4 y Q8. 
De igual manera en dicha referenda se indica la forma como se obtiene la matriz de rigidez de 
un elemento finito saturado, asociado a los 24 grados de libertad indicados a la derecha de la 
figura 3; se destaca que en esta figura con la letra p se ha identificado a las presiones de poro. 

Lo importante es destacar que el elemento finito de la izquierda de la figura 3 es para 
un suelo completamente seco, en cambio el elemento finito de la derecha es para un suelo 
completamente saturado. Evidentemente no se debe utilizar el modelo de la izquierda para el 
analisis estatico de la Presa la Esperanza ya que el suelo esta saturado pero se lo emplea en 
este artfculo para ilustrar el calculo de la matriz de rigidez por ensamblaje directo, ya que son 
menos numeros y de paso para comparar los resultados que se obtienen con el modelo de 
suelo saturado. Ademas de ello en la primera parte del artfculo la carga se aplica en un solo 
escalon de carga. 



3. DESCRIPCION DEL PROBLEMA 

En la figura 4 se muestra la columna de suelo confinado, que se va a analizar, la 
misma que tiene 1.0 m., de base por 5.0 m., de altura y el espesor del elemento es de 1.0 m. 
Se discretiza el dominio en 5 elementos finitos, los mismos que se han numerado de abajo 
hacia arriba, como se ve en la figura 5. Estrictamente las fuerzas horizontales deben actuar en 
el centra de presiones de los elementos pero el objetivo del artfculo es ilustrar el calculo 
manual y posteriormente el uso del GID. 
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Se ha considerado que la columna de suelo confinado, puede ser del Nucleo, Ataguia o 
Espaldon de la Presa la Esperanza. De tal manera que se tienen tres tipos de materiales, los 
mismos que se indican en la tabla 2, en donde se indica el tipo de material; el peso especifico 
del suelo 7; el % de humedad del material % W (no se utiliza en el modelo); la variacion del 

modulo de elasticidad del suelo E ; la variacion del modulo de Poisson del suelo V ; la 
compresibilidad de los granos k s y la permeabilidad del suelo en los dos sentidos, k x - k y . 

De tal manera que se resolvieron tres problemas en todos ellos se trabajo con los valores 
rmnimos y con los valores maximos. 



Tabla 2 Propiedades del suelo. 



Material 


y 

(Suelo) 

[KN/m3] 


% de 
W 




E [KN/m2] 


V 


k s 

Compresibilidad 
De los granos 

[KN/m2] 


k x= k y 

Permeabi- 

lidad 

[m/dia] 


Min. 


Max. 


Min 


Max. 


Nucleo 


15.090 


26 


0.7296 


30600 


40800 


0.25 


0.48 


1142752.048 


4.23E-04 


Ataguia 


14.990 


27 


0.7412 


30600 


40800 


0.25 


0.48 


1157786.982 


2.51E-04 


Espaldon 


15.620 


23 


0.7286 


55100 


138700 


0.16 


0.36 


1141449.313 


3.28E-03 



La descripcion de los materiales del Nucleo y del Espaldon, estan indicados en la Tabla 
1. El material de la Ataguia es similar al Nucleo. Dragados y Construcciones (1994). 

La compresibilidad del agua se considero igual para los tres materiales y tiene un valor 
de k w = 1019700 KN / m 2 (Rigidez del agua). Mroginski (2008). 

4. ANALISIS ESTATICO 

La ecuacion basica de equilibrio para los problemas estaticos es la siguiente: 



Q = Kq (1) 

Donde Q es el vector de cargas generalizadas; K es la matriz de rigidez de la 
estructura, que se halla por ensamblaje directo y q es el vector de coordenadas generalizadas 
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que contendra los desplazamientos en los nodos para el modelo que no considera el agua y 
que contendra los desplazamientos y presiones para el modelo que considera suelo saturado. 

Para encontrar la matriz de rigidez de la estructura K se debe encontrar en primer 
lugar la matriz de rigidez de cada uno de los elementos finitos en coordenadas globales (Aguiar 
y Mroginski 2010) y los vectores de colocacion. Aguiar (2004). 

El vector de colocacion de un elemento finito, esta compuesto por los grados de 
libertad de los nodos del elemento, indicados en forma antihorario, empezando desde cualquier 
nodo pero para sistematizar el calculo es conveniente que se empiece por el nodo inferior 
izquierdo. Para el modelo numerico de calculo indicado a la izquierda de la figura 3 el vector de 
colocacion tiene 8 elementos; en cambio para el modelo numerico indicado a la derecha de la 
figura 3 se tienen 24 elementos. 

Cuando solo hay cargas en las juntas, como en el presente caso, el calculo del vector 
de cargas Q es directo y solo se tendra valores en los puntos de aplicacion de las cargas, los 
restantes valores son nulos. 

En este artfculo, el analisis estatico finaliza con la determinacion del vector de cargas 
q debido a que interesa ilustrar la obtencion de la matriz de rigidez por ensamblaje directo, en 

forma manual (con la ayuda de programas en Matlab) y utilizando el programa GID. 

5. GRADOS DE LIBERTAD DE MODELOS UTILIZADOS 

En la figura 5 se presentan los grados de libertad de los dos modelos utilizados, a la 
izquierda se tiene el primer modelo que no considera la presencia de agua en el suelo y se 
trabaja con un elemento finito Q4, por esta razon se tienen dos grados de libertad por nudo. Se 
ha numerado en primer lugar los desplazamientos horizontales (de 1 a 6) luego los 
desplazamientos verticales (de 7 a 16). 

A la derecha de la figura 5 se muestra el modelo dos para un suelo saturado, en el cual 
se trabaja con el elemento finito Q8. Unicamente para comparar los resultados se numeraron 
en primer lugar los desplazamientos horizontales y verticales de los nodos, en el mismo orden 
del modelo uno luego se ha numerado los restantes grados de libertad. En total se tiene 67 
grados de libertad con el modelo dos; notese que no existen presiones en la parte superior de 
la columna de suelo ya que son nulas por encontrarse a presion atmosferica. 

En las tablas 3 y 4 se indican los vectores de colocacion de los elementos finitos para 
los modelos uno y dos. Como se habfa indicado para el modelo uno, los vectores de colocacion 
tienen 8 elementos, en cambio para el modelo dos, se tienen 24 elementos. 



Tabla 3 Vectores de colocacion de modelo 1 . 



Elemento 


Vector de colocacion 


1 


[1 


0 


0 


0 


0 


8 2 


7] 


2 


[2 


7 


0 


8 


0 


10 3 


9] 


3 


[3 


9 


0 


10 


0 


12 4 


11] 


4 


[4 


11 


0 


12 


0 


14 5 


13] 


5 


[5 


13 


0 


14 


0 


16 6 


15] 
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Tabla 4 Vectores de colocacion de modelo 2 



Elemento 


Vector de colocacion 


1 


[1 


0 


17 


18 


0 


19 


0 


0 


20 


0 


21 


22 


0 


8 


23 


24 


25 


26 


2 


7 


27 


28 


29 


30] 


2 


[2 


7 


27 


24 


25 


26 


0 


8 


23 


0 


31 


32 


0 


10 


33 


34 


35 


36 


3 


9 


37 


38 


39 


40] 


3 


[3 


9 


37 


34 


35 


36 


0 


10 


33 


0 


41 


42 


0 


12 


43 


44 


45 


46 


4 


11 


47 


48 


49 


50] 


4 


[4 


11 


47 


44 


45 


46 


0 


12 


43 


0 


51 


52 


0 


14 


53 


54 


55 


56 


5 


13 


57 


58 


59 


60] 


5 


[5 


13 


57 


54 


55 


56 


0 


14 


53 


0 


61 


62 


0 


16 


0 


63 


64 


0 


6 


15 


0 


65 


66 


67] 



6. RESULTADOS DE SOLUCION MANUAL 

En las tabla 5 se presentan los desplazamientos encontrados, con los dos modelos de 
analisis y para el segundo caso se muestran los resultados hallados con los valores mmirnos y 
maximos, de los datos de suelo indicados en la tabla 2. En la parte superior izquierda se 
muestran los resultados para el caso de que el material de la columna de suelo sea el del 
Nucleo de la Presa la Esperanza; en la parte superior derecha lo propio pero para la Ataguia y 
en la parte inferior para Espalddn. 
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Tabla 5 Desplazamientos encontrados con los dos modelos. 



ATAGLUA 


(Va lores min.) 


(Va lores max.) 


Desplazamientos 


Desplazamientos 


gdl 


Modelo 1 


Modelo 2 


Modelo 1 


Modelo 2 


4 nodos 


8 nodos 


nodos 






-0.0003 


0.0002 


-0.0003 


0.0002 


2 


0.0015 


0.0015 


0.0011 


0.0014 


3 


0.0010 


0.0013 


0.0008 


0.0011 


4 


0.0005 


O.OOOS 


0.0005 


0.0007 


5 


0.0004 


0.0004 


0.0003 


0.0004 


6 


-0.0002 


0.0000 


-0.0001 


0.0000 


7 


0.0003 


0.0004 


0.0003 


0.0004 


3 


0.0000 


0.0004 


0.0001 


0.0004 


9 


0.0003 


0.0011 


0.0006 


0.0009 


10 


0.0006 


0.0011 


0.0007 


0.0010 


ii 


0.0006 


0.0016 


0.0009 


0.0014 


12 


O.OOOS 


0.0016 


0.001 


0.0015 


13 


0.0007 


0.0019 


0.0011 


0.0016 


14 


0.0009 


0.0019 


0.0012 


0.0017 


15 


0.0007 


0.0020 


0.0011 


0.0017 




0.0009 


0.0021 


0.0013 


0.0018 



NUCLEG 


(Va lores min.) 


(Va lores max.) 


Desplazamientos 


Desplazamientos 


gdl 


Modelo 1 


Modelo 2 


Modelo 1 


Modelo 2 


4 nodos 






8 nodos 




-0.0003 


0.0002 


-0.0003 


0.0002 


2 


0.0015 


0.0015 


0.0011 


0.0014 


3 


O.OOIO 


0.0013 


O.OOOS 


0.0011 


4 


0.0006 


0.0008 


0.0005 


0.0007 


5 


0.0004 


0.0004 


0.0003 


0.0004 


6 


-0.0002 


0.0000 


-0.0001 


0.0000 


7 


0.0003 


0.0004 


0.0003 


0.0004 


3 


0.0000 


0.0004 


0.0001 


0.0004 


9 


0.0003 


0.0011 


0.0006 


0.0009 


10 


0.0006 


0.0011 


0.0007 


0.0010 


11 


0.0006 


0.0016 


0.0009 


0.0014 


12 


O.OOOS 


0.0016 


0.001 


0.0015 


13 


0.0007 


0.0019 


0.0011 


0.0016 


14 


0.0009 


0.0019 


0.0012 


0.0017 


15 


0.0007 


0.0020 


0.0011 


0.0017 




0.0009 


0.0020 


0.0013 


0.0018 



ESPALDON 


(Valores min.) 


(Va lores max.) 


Desplazamientos 


Desplazamientos 


gdl 


Modelo 1 


Modelo 2 


Modelo 1 


Modelo 2 




8 nodos 






1 


-0.0002 


0.0000 


-0.0001 


0.0000 


2 


O.OOOS 


O.OOOS 


0.0003 


0.0004 


3 


0.0006 


0.0006 


0.0002 


0.0003 


4 


0.0004 


0.0004 


0.0001 


0.0002 


5 


0.0002 


0.0002 


0.0001 


0.0001 


6 


-0.0001 


0.0000 


0.0000 


0.0000 


7 


0.0002 


0.0002 


0.0001 


0.0001 


3 


-0.0001 


0.0002 


0.0000 


0.0001 


9 


0.0001 


0.0005 


0.0001 


0.0002 


10 


0.0002 


0.0006 


0.0002 


0.0003 


11 


0.0002 


O.OOOS 


0.0002 


0.0004 


12 


0.0003 


O.OOOS 


0.0002 


0.0004 


13 


0.0002 


0.0010 


0.0002 


0.0004 


14 


0.0003 


0.0010 


0.0003 


0.0004 


15 


0.0002 


0.0010 


0.0002 


0.0004 


16 


0.0004 


0.0010 


0.0003 


0.0005 



7. USO DE PROGRAMA GID COMO PRE PROCESADOR 

El programa GID, desarrollado en la Universidad Politecnica de Cataluna, sirve como 
un pre procesador, que genera la malla, los datos de la estructura y archivos de datos que 
seran utilizados por el programa de elementos finitos con el cual se trabaje. Posteriormente 
GID trabaja como un pos procesador para presentar en forma grafica los resultados. 
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A continuacion se indica un pequeno manual de uso del programa GID, con los datos 
de la columna de suelo que se ha venido analizando. 



i. Se ejecuta GID y se selecciona en la parte inferior izquierda el icono de malla con lo 
que aparece lo indicado en la figura 6. 

ii. Se crea linea para lo cual se selecciona el icono que esta en un circulo en la figura 7. 
Posteriormente con el Mouse se dibuja la malla, cuando esta en el ultimo punto 
aparece una pantalla que dice: Procedimiento de creacion de un punto. Pulse la 
imagen para obtener un punto existente o para crear un punto nuevo. Luego hay dos 
dibujos con los tftulos: NO JUNTAR y JUNTAR. Si esta bien se selecciona JUNTAR y 
para terminar se selecciona la tecla Esc. Luego de todo esto se tiene lo mostrado en la 
figura 7. 

iii. Se crea superficie para esto se debe seleccionar primero el icono indicado en un 
circulo en la figura 8. Posteriormente con el Mouse se selecciona toda el area a crear y 
finalmente se pulsa Esc., quedando de esta forma lo indicado en la figura 8. 



Proyecto: UNNAMED 



Archivo Ver Geometrfa Utilidades Datos Malla Calcular Ayuda 



- IDI x| 



00^1 Q QL, ” I 

> H 

p <3 

&zz 

/ 



~IM? 



a 

-u J 

i 

-•& 

£ 5S» 

* 

# 

& 

Qa 



L. 



£$ta es la version de EVALUAClON. 



x=6.853 

y=5.646 



Orden: j 



+L 



/inicio] Elementos Finitos 



JE 7 " 



GiD Proy... Art_03 . doc - Microsoft ... | 

Figura 6 Pantalla de inicio de programa GID. 






iv. Luego se indica el elemento finito que se va a utilizar, para ello se selecciona el icono 
MALLA luego el icono ELEMENTO CUADRATICO y finalmente CUADRATICO. 

v. Luego el tipo de elemento finito con el cual se va a trabajar para esto se selecciona el 
icono MALLA luego el icono TIPO DE ELEMENTO y finalmente CUADRILATERO. Con 
lo que aparece una ventana que dice: Ventana de Informacion. Seleccione Surface 
para asignar el tipo de elemento Quadrilateral y mas abajo esta el icono ACEPTAR que 
se debe pulsar. 

vi. Luego con el Mouse hago click en un punto de la malla (aparecen dos rectangulos). 

vii. En la figura 9 se indica en un circulo el icono, que se encuentra al final de la parte 
izquierda que se debe seleccionar para indicar el tamano de los elementos finitos, al 
seleccionar este icono aparece una ventana que dice: Ventana de Entrada de valor. 
Entrar el tamano de los elementos a generar. En este caso los elementos son de 1 .0 m. 
Se indica este valor con lo cual aparece otra ventana que dice: Ventana de dialogo. 
Malla generada pulsar aceptar para verla. Numero de elementos cuadrilateros = 5. 
Numero de nodos = 28. Al seleccionar ACEPTAR aparece el grafico de la figura 9. 
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Abandonando creacion de linea. 0 nuevas lineas, 0 nuevos puntos — 

Entrar los puntos para deiinir la linea (Escape para terminar) 


i 

■ z=C 


i 

10 


Orden: | §g 
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Figura 7 Creacion de Lineas. 
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Creado una nueva superficie NURBS plana. Puede continuar 
Entrar las lineas para definir la Superficie NURBS (Escape para terminar) 


1 


x=-11.64 

y=-T0.16 

z=0 


' Orden: j 






-f 
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| 0 Luke Filewalker 


J 111 


1- 
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Figura 8 Creacion de superficie. 



viii. 



Para ver la numeracion de los nodos y elementos se utiliza el icono que esta dentro del 
circulo, indicado en la figura 10, es un icono con numeros. Se selecciona el icono y se 
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da enter con lo que aparece la estructura, muy pequena razon por la cual se selecciona 
el icono de agrandar la imagen y se obtiene lo indicado en la figura 10. 
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Figura 9 Malla de elementos finitos con longitudes. 



La numeracion de los nodos es totalmente diferente a la realizada en forma manual, el 
programa GID numera los nodos de tal forma que el ancho de banda de la matriz de 
rigidez de la estructura es la menor posible. 



ix. Ahora se debe especificar que los archivos que se generan con el programa GID se 
almacenen en una carpeta dentro de la particion CONSAT.gid. Para esto se selecciona 
el icono DATOS, luego TIPO DEL PROBLEMA y finalmente CONSAT. En la particion 
CONSAT.gid se encuentran un conjunto de archivos que especifican el formato de 
salida de los datos del preproceso de GID, los cuales deben ser lefdos por el programa 
de calculo de Elementos Finitos. 



Se debe indicar las restricciones de los nudos, para ello se selecciona el icono DATOS 
luego CONDICIONES y finalmente se selecciona PRESCRIBED DISPLACEMENT. 
Luego de hacer todo esto aparece la pantalla indicada en la figura 11. Se seleccionan 
los casilleros de los grados de libertad que se restringen y luego se selecciona 
ASIGNAR. Ahora se deben indicar con un click el nudo que tiene esas restricciones, si 
son varios nudos con la misma restriccion como en el presente caso en que los nudos 
26 y 27 tienen los mismos apoyos, con el Mouse se hace un rectangulo y se 
seleccionan los dos nudos. De forma similar se hace con los nudos 6, 8, 10, 13, 15, 18, 
20, 23 y 25 que tienen los mismos apoyos. 

Si se tuviera un problema con condiciones iniciales, en el casillero para escribir 
cantidades se indicarfan estas condiciones. En la figura 12 se indica la forma como se 
dan las restricciones para los nudos 1 y 3. 

La restriccion del nudo 28 se realiza de forma individual ya que posee dos condiciones 
diferentes a las demas (restriccion en X y en Y). Si se restringe el nudo 28 en X y 
despues en Y, automaticamente el programa opta solo la ultima restriccion (Y). De la 
misma manera se hace con el nudo 5 que tiene restricciones en X y presion. 
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Figura 10 Numeracion de nodos y elementos. 



i «sb <>\ 13:00 




Initio I Elementos Finitos 



| Art_03.doc - Microsoft ... 1 1 ^ C 



16:35 



Figura 1 1 Pantalla para indicar las restricciones de los nudos. 



En la figura 12 falta por seleccionar TERMINAR, con lo que desaparece el rectangulo 
que selecciona a los nudos 1 y 3. 
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Figura 12 Restricciones para los nudos 1 y 3. 

x. Se indican las fuerzas laterales, para ello, bajo la pestana de PRESCRIBED 
DISPLACEMENT se selecciona PRESCRIBED FORCE, se indican las fuerzas 
actuantes en los nodos y se obtiene la figura 13. 
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Figura 13 Asignacion de fuerzas a los nodos. 



i a 



Se selecciona ASIGNAR, luego con el Mouse se indica el nodo que tiene esas fuerzas, 
se lo indica y se coloca TERMINAR. Asi se procede con todos los nodos cargados. 
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A la altura de PRESCRIBED FORCE, se encuentra a la derecha un icono de libreta, el 
mismo que se selecciona y luego se dan los datos de los materiales del perfil de suelo. 
En la figura 14 se indican los datos que se han dado para el presente problema. 
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Figura 14 Propiedades de uno de los suelos. 



Ahora se selecciona ASIGNAR, como todos los elementos tienen el mismo material 
con el Mouse se selecciona un rectangulo y se asigna a toda la malla, finalmente se 
coloca TERMINAR. 



Ahora se dan datos para la forma como se aplican las cargas en los nodos, para esto 
se selecciona CONDICIONES. 

Se graban los datos, para lo cual se da un nombre que sera el de la carpeta que se 
crea en CONSAT. En este caso se ha llamado MANTA1 . 

Finalmente se selecciona el icono CALCULAR. Con lo que se crean las siguientes 
carpetas dentro de MANTA1 . 

Manta1.dat contiene la conectividad de los nodos, de cada uno de los elementos 
finitos. Los datos creados son: 



26 27 28 25 23 22 21 24 1 

21 22 23 20 18 17 16 19 1 

16 17 18 15 13 12 11 14 1 

11 12 13 10 8 7 491 

4 7 8 6 5 3 1 2 1 

Los primeros 8 valores corresponden a la numeracion de los nodos, en forma 
antihorario, empezando por el nodo inferior izquierdo. 

Manta1-1.dat que contiene las coordenadas de los nudos. En total se tienen 28 filas 
de datos, en cada una de ellas esta la coordenada en X, la coordenada en Y., en m. 
Manta1-2.dat que contiene las restricciones de los nodos. Para el presente ejemplo se 
tiene. 
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1 0 


0 


1 


0.00000e+00 


0.00000e+00 


0.00000e+00 


30 


0 


1 


0.00000e+00 


0.00000e+00 


0.00000e+00 


5 1 


0 


1 


0.00000e+00 


0.00000e+00 


0.00000e+00 


6 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


8 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


10 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


13 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


15 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


18 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


20 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


23 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


25 1 


0 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


26 0 1 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


27 0 1 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 


28 1 


1 


0 


0.00000e+00 


0.00000e+00 


0.00000e+00 



El primer dato es el numero del nodo restringido, luego desplazamiento en X, en Y, y 
presion. Para los nodos 1 y 3 la presion no existe. Para los nodos 26, y 27 los 
desplazamientos en Y no existen, en cambio los desplazamientos en X en los nodos 6, 
8, 10, 13, 15, 18, 20, 23 y 25 son nulos. Se observa que en el nodo 5 no existe 
desplazamiento en X y presion, y en el nudo 28 son nulos los desplazamientos en X y 
en Y. 






Manta1-3.dat 

valen: 



que 


contiene las cargas en los nodos. Para el ejemplo estas cargas 


4 


1 .00000e+01 


0.00000e+00 


0.00000e+00 


11 


2.00000e+01 


0.00000e+00 


0.00000e+00 


16 


3.00000e+01 


0.00000e+00 


0.00000e+00 


21 


4.00000e+01 


0.00000e+00 


0.00000e+00 



Se especifica el numero del nodo, la fuerza horizontal, la vertical y la presion. 

• Manta1-4.dat que contiene las propiedades de los materiales. Para el ejemplo se tiene: 

3.06000e+04 2.50000e-01 1.00000e+00 7.29600e-01 1.14275e+06 1.01970e+06 

4.23000e-04 4.23000e-04 0.00000e+00 0.00000e+00 

Los datos que aparecen son \E, v, e, e 0 , k s , k w , k x , k y , 0 0. 

• Manta1-5.dat contiene datos de aplicacion de la carga, la primera fila de datos es: 

0.1 0.100 50 0.01 0.5 

Los datos son escalon de aplicacion de la carga; convergencia; dato numerico de 
integracion; dT el tiempo como se aplica la carga; 0 es para la integracion en el 
tiempo. Las fuerzas horizontales se van a aplicar en 1 0 intervalos; el dT significa 0.01 
de dfa. Hay 9 filas mas de datos con la misma informacion. 

8. PROGRAMA DE ELEMENTOS FINITOS 



8.1 Matriz de rigidez de la estructura sin restricciones de los nudos de borde 

Se indica a continuacion la forma como se obtiene la matriz de rigidez de toda una 
estructura a partir de la numeracion de los nodos que reporta GID. A manera de ejemplo, en la 
figura 15 se muestra a la izquierda la numeracion de los nodos, del elemento finito 4 y a la 
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derecha los grados de libertad asociados a la numeracion de los nodos del GID, los grados de 
libertad de las presiones estan antecedidos con la letra p . 

8 



,10 



11 12 13 

CONECTIVIDAD 

Figura 15 Conectividad y Grados de libertad de los nodos para elemento 4. 

Primero se halla la matriz de rigidez del elemento finito para suelo saturado utilizando 
el programa KelemQ8C0Sat, cuyo marco teorico y forma de uso se detallo en Aguiar y 
Mroginski (2010). Ahora con relacion a la figura 15, en el archivo Manta1.dat, esta la 
secuencia de la conectividad, se inicia en el nodo 11, para este nodo se debe obtener los 
grados de libertad 31 , 32 y 33. Esto se lo realiza con las siguientes ecuaciones: 

filag=(nodoi-1)*ndofn+igl . . . 

colug=(nodoj-1)*ndofn+Jgl ' ' 

Donde nodoi, nodoj, es la identificacion del numero del nodo, para el ejemplo se inicia 
en 1 1 ; ndofn es el numero de grados de libertad que se considera por nodo, en este caso 3; 
igl, jgl son contadores que varfan de 1 a ndofn; filag, colug, identifican los grados de libertad 
de nodo considerado, que va a ser ensamblado. Una ecuacion similar a la (2) se utiliza para 
encontrar la posicion de la fila y columna dentro de la matriz de rigidez del elemento finito, las 
variables utilizadas son: filae, colue. 

Antes de presentar el programa que obtiene la matriz de rigidez de la estructura 
considerando como que si todos los nodos tienen grados de libertad, es necesario indicar el 
significado de las variables utilizadas: 

nelem es el numero del elemento; nnodo es el numero de nodos que se considera en el 
elemento, en este caso es igual a 8; nnodA es el numero de nodos de las esquinas del 
elemento, en este caso es 4; npoin es el numero de nodos de la malla (estructura); nload es 
la cantidad de cargas nodales exteriores; ngaus es el numero de puntos que se consideran en 
la integracion por la cuadratura de Gauss, en cada direccion; ndime es la dimension en la que 
se trabaja el problema, para el caso piano vale 2; ndofn es el numero de grados de libertad 
que se considera por nodo, para este caso es 3; nfixe es la cantidad de grados de libertad 

prescritos; numat es el numero de materiales diferentes de la estructura; nine es el numero 
de incrementos de carga; nstre es para identificar el problema que se esta resolviendo vale 3 
para tension plana y 4 para deformacion plana; ntype es codigo que se utiliza para indicar si 
se trata de un material solido, en este caso vale 1 ; material saturado, vale 2; material no 
saturado, vale 3; en este caso ntype = 2 . 

El programa que obtiene la matriz de rigidez de la estructura, por ensamblaje directo se 
denomina EnsamblajeK y los datos de entrada son: 





‘32 p 36 34 38T 

GDL 
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coord Matriz de coordenadas nodales; conec Matriz de conectividades; matpr propiedades 
de los materials; timer propiedad de la integracion en el tiempo; despG es el vector que 
contiene a los desplazamientos nodales y presiones de poro en coordenadas globales del 
incremento anterior de carga; FGlob es el vector de cargas generalizadas contiene a las 
cargas actuantes en los nodos y a las cargas debidas a las presiones de poro; loads es el 
vector de cargas nodales en coordenadas globales; inc es el numero de incremento de carga 
(forma de aplicar la carga). A continuacion se lista el programa. 

function [KGlob, FGlob, FGele] = EnsamblajeK(coord, conec, matpr, timer, despG, ... 

FGlob, loads, inc) 

% Ensambla la matriz de Rigidez KGlob y el vector de fuerzas FGlob 

global nelem nnodo nnod4 npoin nload ngaus ndime ndofn nfixe numat nine nstre ntype 

KGlob = zeros(ndofn*npoin,ndofn*npoin); %M.R. de toda la estructura. 

FGele = zeros(ndofn*npoin,1); 
for ielem=1:nelem 

[Kelem,Felem] = KelemQ8C0Sat(ielem, matpr, timer, conec, coord, inc);%M.R. de E.F. 
for inodo=1:nnodo 

nodoi=conec(ielem,inodo);% Se identifica el nodo del E.F. (Fila) 
for igl=1 :ndofn % ndofn es el numero de grados de libetad por nudo 
filag=(nodoi-1)*ndofn+igl;% Se identifica el gdl asociado al nodoi 
filae=(inodo-1)*ndofn+igl;% Posicion en M.R. de E.F. de la fila 
for jnodo=1:nnodo 

nodoj=conec(ielem,jnodo); %Se identifica el nodo del E.F. (Columna) 
for jgl=1 :ndofn 

colug=(nodoj-1)*ndofn+jgl; % Se identifica el gdl asociado al nodoi 
colue=(jnodo-1)*ndofn+jgl; % Posicion en M.R. de E.F. de la columna 
KGIob(filag,colug)=KGIob(filag,colug)+Kelem(filae,colue); %Ensamblaje 
end 
end 

FGele(filag) = FGele(filag) + Felem(filae) ; 
end 
end 
end 

return 



8.2 Matriz de rigidez de la estructura con restricciones de los nudos de borde 

La estructura que se esta analizando tiene 28 nodos, con el programa EnsamblajeK se 
encuentra una matriz de rigidez de la estructura de (28*3,28*3), es decir una matriz de rigidez 
de 84 por 84. Pero en la figura 5 se aprecia que la estructura tiene 67 grados de libertad, luego 
la matriz de rigidez tiene que ser de ese orden. 

Para encontrar la matriz de rigidez de 67 por 67, se necesitan las restricciones de los 
nodos, que para el ejemplo se hallan en el archivo: Manta1-2.dat. En base a este archivo se 
determina los nodos que no tienen grado de libertad y se procede a eliminar la fila y la columna 
que esta asociado al supuesto grado de libertad que se obtuvo con el programa EnsamblajeK. 
Esto se lo hace con el programa CondBorde, que se lista a continuacion. 



function [KGcon,FGcon, FGlob, ifixe,ndesp] = CondBorde(KGIob, FGlob, fixed, ndesp,iiter) 

global nelem nnodo npoin nload ngaus ndime ndofn nfixe numat nine nstre 

ifixe=zeros(ndofn*npoin,2); 

for i=1:nfixe 

npos = (fixed(i,1)-1)*ndofn + 1 ; 

j=o ; 

for idof=npos:ndofn+npos-1 

j=i+i ; 

ifixe(idof,:) = fixed(i,[j+1 j+ndofn+1]);% Se coloca 1 en los nudos restringuidos 
end 
end 
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nfix=find(ifixe(:,1)==0); %Vector con los gdl de la estructura, sin restricciones 
KGcon=KGIob(nfix,nfix); %Se encuentra la matriz de rigidez para los gdl de la estructura 
fix=find(ifixe(:,1)==1); %Se identifica los gdl de los nudos restringuidos 

if iiter==1 

FGlob = FGlob - KGIob(:,fix)*ifixe(fix,2)./ninc;% nine es los incrementos de carga 
ndesp = find(FGIob(:,1)~=0);% Identifica los gdl donde actua la carga y .... 
end 

FGcon = FGIob(nfix); % Vector de cargas en los gdl reales de la estructura. 
return 



Los datos de ingreso al programa CondBorde, son: KGlob , FGlob que contiene la 
matriz de rigidez y el vector de cargas sin tomar en cuenta las restricciones de los nodos; 
fixed es la matriz que contiene las restricciones de los desplazamientos y presiones 

prescritas, para el ejemplo es el archivo Manta1-2.dat, ndesp es un vector que identifica el 

grado de libertad donde actua la carga y donde la fuerza es diferente de cero, dado que el 
problema acoplado de consolidacion en suelos saturados es no lineal, el proceso de solucion 
debe ser iterativo, liter define el numero de iteraccion que se esta realizando. 



A la izquierda de la figura 16, se muestra la numeracion de los nodos que definio el 
programa GID y a la derecha se tienen los grados de libertad con los cuales se resuelve la 
estructura. Notese que son diferentes a los utilizados en la solucion manual, indicada en la 
figura 5. 
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Figura 16 Numeracion de nodos y grados de libertad de estructura analizada. 



8.3 Vector de cargas en estructura sin restriccion de los nudos de borde 

En el programa GID se identifican los nodos donde actua la fuerza horizontal, la fuerza 
vertical y la presion. Ahora bien, en este apartado se indica como se obtiene el vector de 
cargas en la estructura sin tener en cuenta las restricciones de los nudos, para cuando solo 
actuan las fuerzas horizontal o vertical. 
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Sea inodo la identificacion del nodo cargado, con la ecuacion (3) se determina el 
primer grado de libertad de ese nodo. 

npos=(inodo-1)*ndofn+1 ( 3 ) 

Donde ndofn es el numero de grados de libertad que se considera por nodo, para el 
presente caso es 3; npos es la identificacion del grado de libertad horizontal del nodo cargado. 
El programa que calcula el vector de cargas sin considerar las restricciones de los nudos se 
denomina EnsamblajeF y se lo lista a continuacion. 

function FGlob = EnsamblajeF(loads,FGIob,Carga,nload, ndofn) 

%%% Ensambla el vector de fuerzas de las cargas exteriores %%% 

for ifix=1 :nload 

inodo=loads(ifix,1);% Identifica el nodo cargado de acuerdo a numeracion GID 
npos=(inodo-1)*ndofn+1 ;% Identifica el grado de libertad para fuerza horizontal 
FGIob(npos,1)=FGIob(npos,1)+loads(ifix,2)*Carga;% Coloca la fuerza horizontal 
npos=npos+1 ; % Identifica el grado de libertad para fuerza vertical 
FGIob(npos,1)=FGIob(npos,1)+loads(ifix,3)*Carga;% Coloca la carga vertical 
end 



8.4 Vector de cargas por la presion de poro 

Sea Q t el vector de cargas. Cuando recien se inicia el problema el vector de cargas 
esta formado exclusivamente por las cargas laterales o verticales pero con el transcurso del 
tiempo, este vector se incrementa en nuevas filas, debido a la presion de poro P PORO por la 

presencia del agua. Se destaca que se incrementan las filas asociadas a los grados de libertad 
de las presiones y se obtiene con la siguiente ecuacion. 

Qt =Q t + dT * H ww * P PORO ( 4 ) 

Donde dT es el tiempo como se aplica la carga; H ww es la matriz de rigidez del 

agua de 4 por 4; P PORO es el vector que contiene las presiones de los poros en las esquinas 

del elemento finito. El programa que calcula el vector de carga por la presion de poro se 
denomina: FGFlux y se lo lista a continuacion. 

function [FGlob] = FGFIux(FGIob, coord, conec,despG, timer, matpr,inc) 

% Calculo del vector de cargas por la presion de poro. 

global nelem nnodo nnod4 npoin nload ngaus ndime ndofn nfixe numat nine nstre ntype 
[gausp,weigp] = gauss_points(ngaus); % Valores de la cuadratura de Gauss y sus pesos 
DTime = timer(inc,4); % es el tiempo de aplicacion de la carga 
for ielem=1:nelem 

corlo = coord(conec(ielem,1 :1 :nnodo),:); %Coordenada de los nodos del E.F. 
corl4 = corlo([1 3 5 7],:); %Coordenadas de las esquinas del E.F. 
imat=conec(ielem,9); %ldentificacion del Tipo de Material 

Kx = matpr(imat,7)./10 ; Ky = matpr(imat,8)./10 ; %se divide por 10 por las unidades 
Kperm = [Kx 0.0 ; 0.0 Ky ]; %Matriz de permeabilidad 

npos = (conec(ielem,1 :2:nnodo)-1).*ndofn+ndime+1 ; %Numero de nodos esquineros para Presiones 
Pporo = despG(npos); %Presion de Poro de los nodos esquineros 
Hww = zeros(nnod4,nnod4); %Matriz de rigidez del agua 
for igaus=1:ngaus 
psi = gausp(igaus) ; 
for jgaus=1:ngaus 
eta = gausp(jgaus); 

[shap4,deri4] = ShapeQC0(psi,eta,nnod4) ; 

[detJB,DerCa,MatJB] = Jacobian_Matrix(corl4,shap4,deri4, ndime, nnod4) ; 
difV4 = detJB*weigp(igaus)*weigp(jgaus)*matpr(imat,3) ; 
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Hww = Hww + DerCa*Kperm*DerCa'.*difV4 ; 
end %jgaus 
end %igaus 

FGIob(npos,1) = FGIob(npos,1) + (DTime.*Hww*Pporo); % Incremento por la presion de poro 
end 



Este programa llama al programa gauss_points que se lo lista a continuacion. 

function [gausp,weigp] = gauss_points(ngaus) 

if ngaus==1 
gausp(1)=0.0; 
weigp(1)=2.0; 
elseif ngaus==2 

gausp(1)= -1/(sqrt(3)); %-0.5773502691 89626 

gausp(2)= 1/(sqrt(3)); %0.5773502691 89626 

weigp(1)=1.0; 

weigp(2)=1.0; 

end 

return 



8.5 Programa Principal 

Una vez que se tiene la matriz de rigidez y el vector de cargas se obtienen los 
desplazamientos horizontal y vertical y las presiones de poro, utilizando la ecuacion (1), que 
era el alcance del articulo sin embargo de ello a continuacion se presenta el programa 
completo de consolidacion desarrollado por Mroginski (2008). 

% Programa para el calculo de suelos saturados 2D con elementos finitos 
% Pre-post procesador: GID 9.0.2 
% Autor: Mg. Ing. Javier L. Mroginski 
% Fecha: 14/06/2010 

clc 
clear 
close all 

global nelem nnodo nnod4 npoin nload ngaus ndime ndofn nfixe numat nine nstre ntype 
displayf MEF 2D - elementos Q8C1') 

problem^Mantal'^/o'QS-lxISat'; % , Unidim' %'Zapata' % 'Q8-1x1Sat' % , Pamin5.3.1'; % Nombre del problema a 
ser ejecutado 

%[conec, coord, fixed, loads, matpr, timer, exp] = input(problem) ; % Lee los ficheros de entrada 

% INPUT 

conec = load([problem '.dat']) ; % Matriz de conectividades 

coord = load([problem '-1' '.dat']) ; % Coordenadas nodales 

fixed = load([problem '-2' '.dat']) ; % Restricciones/Desplazamientos prescriptos 

loads = load([problem '-3' '.dat']) ; % Cargas Nodales en coord globales 

matpr = load([problem '-4' '.dat']) ; % Propiedades Materiales 

timer = load([problem '-5' '.dat']) ; % Propiedades de la integracion en el tiempo 

nload = size(loads,1); % cantidad de cargas nodales exteriores 

nelem = size(conec,1); % cantidad de elementos 

nnodo = size(conec,2)-1 ; % numero de nodos del elemento 

nnod4 = 4 ; % nodos de esquina (elemento Q8, 2D 

ngaus = 2 ; % numero de ptos. de Gauss por direccion 

ndime = size(coord,2); % dimension del problema 

npoin = size(coord,1); % numero de nodos totales de la malla 

nfixe = size(fixed,1); % cantidad de grados de libertad prescriptos 

numat = size(matpr,1); % numero de materiales distintos 

nine = size(timer,1); % numero de incrementos de carga 

nstre = 4 ; % 3: para EPT , 4: para EPD 

ntype = 2 ; % 1 : Solidos , 2: Suelos Saturados , 3: Suelos No Sat. 

ndofn = ndime+(ntype-1); % grados de libertad por nodo 
DESP = zeros(ninc , 1) ; 
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PPOR = zeros(ninc , 1) ; 

EXP = [DESP PPOR] ; 

% 

StresG = zeros(4,ngaus*ngaus*nelem) ; StresT = zeros(4,ngaus*ngaus*nelem) ; 

StraiG = zeros(4,ngaus*ngaus*nelem) ; StraiT = zeros(4,ngaus*ngaus*nelem) ; 
despG = zeros(ndofn*npoin,1) ; despT = zeros(ndofn*npoin,1) ; 

Tload = zeros(ndofn*npoin,1) ; Eload = zeros(ndofn*npoin,1) ; 

FGlob = zeros(ndofn*npoin,1) ; 

ndesp = []; safe=[] ; DESP=[]; 

figure(l) ; 

plotter(coord, conec, nelem,nnodo,'b-') %grafica la malla indeformada con azul 

figure(2) ; hold on ; grid ; 
figure(3) ; hold on ; grid ; 

Carga=0; tol=0; maxit=0; DTime=0; Theta = 0; Time=0.0; 
for inc=1:ninc 

Carga = timer(inc,1) ; tol = timer(inc,2) ; maxit = timer(inc,3) ; 

DTime = timer(inc,4) ; Theta = timer(inc,5) ; 

fprintf Incremento de carga : %3i \t \n',inc); 

deAnt = despT ; iiter = 0 ; conv = 1.0; Time = Time+DTime; 

FGAnt = FGlob ; 

FGele = zeros(ndofn*npoin,1) ; 

FGlob = EnsamblajeF(loads, FGlob, Carga, nload,ndofn) ; 

FGlob = FGlob - FGAnt ; 

[FGlob] = FGFIux(FGIob, coord, conec, despT, timer, matpr,inc) ; 
while conv==1 

iiter = iiter + 1 ; fprintf (’> Iteraccion : %3i \t \n', iiter) ; 
if iiter>=maxit+1 

disp(’> Se supero el maximo nro de iteraciones ') 
sound(wavread('aycaramba.wav'), 23000) ; 
return 
end 

[KGlob, FGlob, FGele] = EnsamblajeK(coord, conec, matpr, timer, despG,. .. 

FGlob, loads, inc); 

[KGcon,FGcon, FGlob, ifixe, ndesp] = CondBorde(KGIob, FGlob, fixed, ndesp, iiter); 
despl=KGcon\FGcon ; 

despG(find(ifixe(:,1 )==1 ))=ifixe(find(ifixe(:,1 )==1 ),2).*(1/ninc) ; 
despG(find(ifixe(:,1 )==0))=despl ; 

[StresG, StraiG] = StressReductionEL(despG,matpr, conec, coord, StresT, StraiT) ; 

[ResEn] = Conver_Energy(despG, conec, coord, matpr, StresG, loads, ... 

FGlob, FGAnt, ifixe) ; 
if ResEn<=tol %ResPI<=tol|ResPS<=tol 
conv=0; 
else 

FGlob = FGlob - FGele ; 
end 

fprintf (’> Residuo ENERGY: %10.5e \t \n', ResEn); 
end %iiter 

StresT = StresG ; StraiT = StraiG ; 

despT = deAnt+despG; %StresT = StresT+StresG ; StraiT = StraiT+StraiG ; 
coord2 = actual_coord(coord,npoin,ndofn, despT); 

figure(l) ; plotter(coord2, conec, nelem,nnodo,'r-') %grafica la malla deformada con rojo 

inodo = 8 ; %4 ; %8 ; % 237 ; % nodo a graficar 
gl = 2 ; % 1=x, 2=y, 3=L, 4=Lx, 5=Ly, 6=Lxy 

if strcmp(problem, , Unidim2D')==1 
inodo = 4 ; gl = 2 ; 

elseif strcmp(problem,'Q8-1 xl Sat')==1 

inodo = 1 ; gl = 2; % para presion de poro (Fig 3) 
end 

[safe] = GraficoXY(inodo,gl, Time, EXP, coord, matpr, npoin,ndofn, inc, despT, StraiT, StresT, KGlob, safe) ; 

GID_result(inc, problem, despT, StraiT, StresT); %RESULTADOS PARA GID 
pause(O.I) 
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end % inc 

hold off ; %figure 2; hold off ; %figure 3 
coord = actual_coord(coord,npoin,ndofn,despG); 

figure(l) ; plotter(coord,conec,nelern,nnodo, , r-') %grafica la malla deformada con rojo 
hold off ; %figure 1 

GID_mesh(problem, coord, conec); % Geometria para el postproceso del GID 9.2 

dispC****** THE END ******’) 

9. PROGRAMA GID COMO POST PROCESADOR 

Una vez que se ejecuta el programa de elementos finitos se deben generar archivos 
con la extension FLAVIA.RES y FLAVIA.MSH. El nombre que antecede a estos archivos es 
aquel que se lo dio cuando se utilizo el programa GID como Pre Procesador. Estos dos 
archivos se crean directamente dentro de la carpeta Mantal, para el ejemplo que se esta 
desarrollando, si no es asi se deben copiarlos e insertarlos en esta carpeta, que forma parte de 
CONSAT. 

> Ahora se ingresa al programa GID, se abre la carpeta Mantal y se cambia la 
modalidad de trabajo de Pre Procesador del GID a Post Procesador seleccionando 
la tecla que esta senalada con una flecha roja, en la figura 17. (Ver parte superior 
el 9 icono). 




Figura 17 Se cambia de Pre Procesador a Post Procesador 

> Para ver resultados se sigue la secuencia indicada en la figura 18. Se selecciona el 
icono Ver resultados, luego la forma como desea ver esos resultados, en este 
caso se ha seleccionado Areas coloreadas, luego se indica que tipo de resultado 
se desea ver, en este caso Displacement y finalmente el sentido de los 
desplazamientos, en el ejemplo se ha seleccionado X Displacement. Al dar click 
se encuentran los desplazamientos mostrados en la figura 19. 

En la figura 19 se aprecia que los mayores desplazamientos se tienen en los 
puntos donde se aplican las cargas y se van atenuando hacia la derecha debido a 
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que en el extremo se tiene confinado el suelo. Las mayores cargas laterales se 
colocaron conforme se incrementa la altura en consecuencia los mayores 
desplazamientos laterales se tienen en la parte inferior. 

Los datos que se utilizaron son los correspondientes a un perfil de suelo semejante 
al suelo del Nucleo de la Presa la Esperanza y con los valores mfnimos, de tal 
manera que para comparar estos resultados con los indicados en la tabla 5 se 
deben ver los mostrados en la parte superior izquierda para el modelo de 8 nodos, 
para los grados de libertad 2, se obtuvo 0.0016 m., para el grado de libertad 3 se 
hallo 0.0012 m., para el grado de libertad 4 se encontro 0.0008 m., y para el grado 
de libertad 5 se tuvo 0.0004 m. 




f Inicioj Mantal.gid | ._J Elementos Finitos | iig] Art_03.doc -Microsoft... | \j Dibujo- Paint || t ' GiD Pro... g S3 ^ 17:00 

Figura 18 Procedimiento a seguir para ver resultados. 



En la parte inferior de la figura 10 se indica que los desplazamientos encontrados 
son para el paso 10. Se recuerda que se dieron 10 escalones de carga. 

Para ver las cantidades indicadas de los desplazamientos se ingresa al archivo 
Mantal.flavia.res, se visualiza los resultados con Word, agrandando la 
presentacion . Ahf aparece el paso de carga y los valores encontrados. 
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Figura 19 Desplazamientos horizontales de columna de suelo confinada. 



1 0. COMENT ARIOS Y CONCLUSIONES 

Se ha encontrado la respuesta estatica de una columna de suelo confinado, sometido a 
presion hidrostatica del agua, empleando dos modelos de elementos finitos, el Q4 y el Q8; en 
el primer modelo no se considera la presencia del agua, en cambio en el segundo si se lo toma 
en cuenta, la solucion se lo ha realizado aplicando toda la carga en un solo instante de tiempo. 
Los resultados obtenidos se muestran en la tabla 5 para los tres tipos de suelo y para las dos 
condiciones de valores minimos y maximos; se observa que los desplazamientos encontrados 
con los dos modelos son semejantes. 

Posteriormente, se ilustro el uso del Programa GID, en las etapas de Pre y Post 
Proceso, con los datos de la columna de suelo confinado y se presento los programas con los 
cuales se realiza el analisis estatico. Se aspira de esta manera haber contribuido al desarrollo y 
ensenanza de los elementos finitos en suelos saturados. 
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